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On the S bility of Viscous Fluid Motions 


JAMES SERRIN 


Communicated by C, TRUESDELL 


This paper deals with one aspect of the classical problem of hydrodynamical 
stability, namely, the determination of sufficient conditions for a basic (usually 
laminar) flow of an incompressible fluid to be stable under arbitrary disturbances. 
The technique to be applied is the well known method of energy, originated by 
OSBOURNE REYNOLDs and WILLIAM MCF. Orr and used since that time by 
many other writers}. In spite of this intensive study, it appears that a number 
of new results can be secured from the method, and it is to these that the paper 
is devoted. 

Our main conclusion is a Reynolds number criterion for the stability of an 
arbityary fluid motion in a bounded region. In particular, we show that a basic 
flow in a bounded region Y is stable whenever its Reynolds number Re = V d/y 
is less than 5.71; here V is the maximum speed of the basic flow, d is the diameter 
of ¥, and »v is the kinematic viscosity of the fluid. The number 5.71 in itself 
is neither especially good nor especially bad as a criterion for stability, but 
what is interesting is the fact that it is absolutely rigorous and applies in- 
dependently of the geometry of the flow region and the particular flow involved. 
For this reason the result may appropriately be called a Reynolds number for 
universal stability. We also obtain similar criteria for the stability of flows in 
unbounded regions, and applications are made to the problem of uniqueness of 
steady flow. 

In the second part of the paper (§ 4) we state a general variational problem 
connected with the stability of an arbitrary motion. The Euler-Lagrange equa- 
tions corresponding to this problem bear an interesting and remarkable resem- 
blance to the Navier-Stokes equations, but they are in general too difficult to 
solve, except for special cases. 

The last part of the paper treats a particular example, the laminar Couette 
flow between rotating coaxial cylinders. The methods of the earlier sections of 
the paper, when applied to this case, yield the stability criterion 


2.—2 
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Shake { = V 4 
ae 7 R, R, log (R2/Rj) (1) 
in which the notation is practically self-explanatory. The author knows of no 
other formula of this type which applies to arbitrary disturbances of the basic 
laminar motion. (The criteria of G. I. TAyLor and others (cf. [14], Chapter 2) 


1 Cf. references [1]—[9] at the end of the paper. 
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all refer to infinitesimal disturbances of a special type; thus, although they have 
indisputable importance if one is considering the breakdown of already established 
flow, and give quite sharp results in certain limited situations, they do not by 
any means apply to the whole realm of possible disturbances of the motion.) 
Finally, the stability of Couette motion is considered from the point of view of 
the variational method noted above. The results here do not have the finality 
of criterion (4), but they are nevertheless of interest. 

The paper begins with a derivation of a fundamental identity for the rate of 
change of energy of a perturbation motion. This formula is the basis for all that 
follows. 

1. The Reynolds-Orr energy equation 

We consider a basic fluid motion occupying a region ¥ =¥Y (t) of space and 
subject to a prescribed velocity distribution on the boundary Y of ¥. In the 
cases of greatest interest Y is bounded by material walls and the boundary 
conditions arise from the motion of these walls, as, for example, in the case of 
Couette flow. Now suppose the velocity field of the basic flow is altered at 
some initial instant t=O; it is natural to ask whether the subsequent motion, 
subject to the same boundary conditions, will alter only slightly from what it 
was, or whether it will change radically in character. To investigate this problem, 
we shall consider the energy “ of the perturbation (difference) motion. If “ 
tends to zero as too, then the basic motion is said to be stable, or, more 
precisely, stable in the mean. 

To be specific, suppose the region ¥ is bounded, and let v and wv’ denote, 
respectively, the velocity vectors of the basic and altered motions. The velocity 
u=v’—v of the perturbation motion obviously satisfies 


u=0o on F, (2) 


and its kinetic energy is given by 
H=1*[ He (3) 


(in writing integrals, we shall consistently omit the conventional volume infini- 
tesimal ; moreover, all integrals are understood to be extended over the entire 
region ¥’). The rate of change of # is governed by the fundamental formula 


Ak 
ae — f{ (u-Deu + grad u: grad u) (4) 


essentially due to REYNoLps and Orr. In this equation grad u denotes the 
matrix with components (grad u);,—=«,,; and D denotes the deformation matrix 
of the basic motion, D;,=4(v, , +0, ): 

To prove (4), we begin with the observation that both w and wv’ satisfy the 
Navier-Stokes equation 


Ow 
oar + + grad v) = ef —grad p + eV? De. 


By subtraction there arises 


du 4 
“gp 7 Us grady + v’« gradu = grad et tyVe. 
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Forming the scalar product of this equation with the vector uw, and using the 
incompressibility conditions div v =div v’=0, then leads to 


Gifs ; 
gala ©) =~ (ee Dow + v gradu: gradu) + div®, (5) 
where , 

= Pik . tla \ Sepals ote 
@ Saye u +y grad (50) 7 ee. 


On the other hand, we have obviously 


dk Ca ; 
dt =fal u) ! howe, (6) 


the latter integral being taken over the entire boundary of Y. The required 
formula (4) now follows easily from (5), (6), the divergence theorem, and the 
fact thatu=—P —( on 4.2 

This derivation does not hold when the region Y is unbounded, since both 
formula (6) and the divergence theorem are applicable only to bounded regions. 
Under suitable conditions on the asymptotic behavior of v, v’, however, the 
above formal steps can still be justified (we shall omit the details). Another 
justification of (4) for infinite regions is available whenever the flow geometry 
is such that the disturbances can be assumed spatially periodic at each instant. 
This will be the case, in particular, for the important Poiseuille and Couette 
flows. The disturbance wu being supposed periodic in the direction of the axis 
of symmetry (what Lin calls sustained oscillations), the region YW can then be 
chosen to cover exactly one period. Then the boundary integrals at either end 
of ¥,, neither of which vanishes separately, just cancel one another. Formula (4) 
may therefore be assumed to hold in these two important situations. 


2. Criteria for universal stability 

In this section we shall use the method of energy to establish certain criteria 
for the stability of arbitrary fluid motions. This method is based on the observa- 
tion that if # tends to zero, then u must likewise tend to zero almost everywhere. 
Thus a basic flow will be stable (stable in the mean) provided the energy of any 
disturbance tends to zero as ¢ increases. To apply the method one considers 
the right-hand side of (4): if it is negative for arbitrary non-vanishing vectors u 
satisfying div w=0, then obviously @.%/dt<0 and there is stability. Since the 
second term on the right of (4) is always negative, it is seen that viscosity tends 
to damp out any disturbance. On the other hand, a high rate of shear in the basic 
flow can cause the first term to be highly positive, thus fostering the growth 
of a disturbance. The relative importance of these two terms, then, determines. 
the stability of the flow (cf. also the discussion in [14], § 4.5). 

Another criterion of the same sort arises when the right-hand side of (4) is 
written in slightly different form. Indeed, since div u=0, we have 


u-D-u =div|(u-v) u| —u-graduev, 


2 Several alternate and occasionally useful forms of (4) arise from the simple 
identity { stad c= grad u— (curl e)* = 2,/ D’: D’, 
where D’ denotes the deformation matrix of the perturbation motion. 


4* 
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whence, by application of the divergence theorem, (4) can be written in the form 


O% = [ (wegraduev —y gradu: grad w). (7) 


If the first term on the right of (7) is less than the second for all admissible 
vectors u, then clearly the basic flow will be stable. But the size of the first 
term is governed by the magnitude of v, from which it follows that high speeds 
in the basic flow, as well as high rates of shear, tend to cause instability. The 
qualitative nature of these effects will be investigated in the next several para- 
graphs. 

It is important to note that the energy method cannot provide accurate 
knowledge of the limits of stability, such as can be gained from the linearized 
perturbation theory ({6], Part II; [14], Chap.1). The reason is that in the 
energy method one establishes stability relative to arbitrary disturbances, while 
in reality only those satisfying the hydrodynamical equations need be considered. 
Nevertheless, because the energy method gives insight into the physical situation, 
and because the results have the merit of simplicity and complete mathematical 
rigor, the investigations based on it are both interesting and valuable. 

With these preliminary remarks aside, we may now turn to the main result 
of the paper. 

Theorem 1. Let 7 =VY (t) be a bounded region of space, which can be included 
in some cube of edge length d. Let v be the velocity vector of any motion in V 
satisfying prescribed conditions at the boundary of ¥. Then the kinetic energy of 
an arbitrary disturbance motion u=v'—vw satisfies the inequalities 


Mi, eo re (8) 
HEH, Val) ty (9) 
Here %, is the initial energy of the disturbance, —m is a lower bound for the 


characteristic values of the deformation matrix of the basic flow over the time interval 

0 to t, V ts the maximum speed of the basic flow in the same time interval, and « 
as a pure number, 

2 a 

a = SEN3 aes 32.6. (10) 

This theorem is a generalization and improvement of certain results of 

T. Y. THomas and EBeruarD Hopr3. Before giving the proof of the theorem, 

we note two important corollaries, which in a certain sense constitute the heart 


of the paper. It is to be emphasized that these corollaries are absolutely rigorous, 
and apply to all possible motions in V. 


Universal stability criterion I. J} the dimensionless “Reynolds number” 
ma*|y of a flow is less than 32.6, then &#—>0 as t—> oo, and the motion is stable. 


3 THOMAS [8] proved that an arbitrary fluid motion is stable if its deformation 
matrix is sufficiently small, but did, not compute numerical values; the present 
analysis is an extension of Tuomas’ work, Similarly, Hopr has shown (essentially) 


that a fluid motion is stable if its maximum speed is small enough, again with no 
discussion of numerical values. 
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Universal stability criterion II. If the Reynolds number Vdlv of a flow is 
less than 5.71, then X-—>0 as t-> co, and the motion is stable. (Here V is the 
maximum speed of the flow, and d is the maximum diameter of the flow region.) 

Proof of Theorenv 1. We begin with an auxiliary computation whose goal.is 
the i lit 

Ad x ne ad*fut< f gradu: gradu, (14) 


where « is given by (10). Let h be an arbitrary continuously differentiable vector 
field in ¥. Then for any value of the constant «, 


OS (4; My, + U5, + Up, 5) (Uz hy + My, ~ + E Up, 5) 
= (1 + 8) UU bP +h, (u?) , + 28 (Uz hy + Uj, 5) Up, 3 


(For convenience we use tensor rather than vector notation in carrying out this 
and several subsequent calculations.) Some of the terms on the right-hand side 
of (12) can be transformed as follows, 


h, (wv) = (wv hi) i Zs uh; ; 


U; hy, UN oa (u; Ny Up) 5 — 4; hy, 5 Up, 


(12) 


Uj; np UR i = (u; Uni), h- 
Here we have made use of the incompressibility condition uv; ;=0. Making the 
above changes in (12), integrating the result over W, and using the divergence 
theorem, we obtain the inequality 
JS (Qi; —?) + 26%; i; , Up] SS (A + €*) J Us, Uk: (13) 


Now the particular vector field 4;=C tan C x; is differentiable in a cube of edge 
length z/C centered at the origin. If we set C=a/d and suitably locate the 
origin in Y, then this vector can be substituted into (13). An easy computation 
then yields the inequality 
(+243 +2) C fs su pmp- 
The left-hand side is maximized by choosing ¢ = }(|/13 — 3), and (41) is thereby 
proved. 
The term u-D-w in (4) satisfies the inequality 


ueDeu=—muv (14) 


during the entire time interval 0 to ¢. Therefore, from (4), (11) and (14) 


Ak oe =) ioe — —-2)\ - 
a =m avd ) fu =2(m—avd) ZX. 


Writing this in the form 
d y —2(m—avd-)t 
5, Goll oA te tes)? ve 0 


and integrating from 0 to ¢ leads at once to (8). 
The proof of (9) is similar but requires in addition the inequality 


u-graduev Ss sl grad wu: grad u +="), (15) 
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which follows at once from the identity 


A:A—2u-A-04+ v= (A—Uv):(A—Uuv) 20 


with A=y gradu. From (7), (15), and (11) we obtain 


A 1 = 
ee 2 2 7-2 
Wie a (V avd Wes 


and (9) then follows by integration as in the proof of (8). 


3. Extensions and applications of Theorem 1 
The preceding method can be applied equally well to flows which take place 
in channels and pipes, and to plane flow problems, so long as the fundamental 
formula (4) applies. To illustrate the method, consider for example a straight 
pipe whose cross section has a maximum diameter d. If the pipe is directed along 
the z-axis, then the vector 


h =C(itanCx+jtanCy), (C= melas 
can be used in (13). Choosing e=0, we get an inequality of the form (11) with 
a —=2z%. By using this result, the reader can easily obtain estimates for the rate 


of change of .#, and thus formulate analogues of Theorem 1 and the universal 
stability criteria for flow in a pipe. 


Other cases can be treated similarly and the results conveniently grouped 
in the following table. 


A. Straight channel, maximum width d: Ci rn al 
B. Straight pipe, maximum diameter d: oS: 


C. Plane flow in a bounded region, maximum diameter d, 


i) Three dimensional disturbances: = 225 
ii) Plane disturbances only: a = (1 + 2) 2? 
. 43 
D. Bounded region, maximum diameter d: oo Sai. ee 
2 


(Note: It is only in cases C, ii) and D that we have been able to make use of the 
incompressibility condition div w= 0.) 


We conclude the section with two simple applications of Theorem 1. First, 
suppose the boundaries of Y consist of rigid fixed walls, so that any motion 
initially present will presumably die out due to lack of an energy source. By 
choosing v=0 for the basic motion, we see from (8) that the energy % of an 
arbitrary motion v’ in Y must in fact tend to zero according to the law 


TATA Oe = 3216. 


Similar estimates, but with smaller values for the coefficient a, have been obtained 
by Leray (for plane flow) and Kampé pe Ferret? and BERKER (for spatial 


4 KaMPE DE FERIET obtained the 2 i 
a value «= 322, using a method somewhat simi 
to the one presented here. . ae 
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flows), and RAYLEIGH in a much earlier paper [4] proved that #—>0 as t-> ov, 
though without estimating the rate of convergence. It remains an open question 
whether the velocity itself must tend to zero as t->0oo; certainly one would 
expect this, but a strict proof seems to be a matter of more than ordinary diffi- 
culty. 

As a second application, we have the following uniqueness theorem concerning 
steady motion in a fixed bounded region . 


Theorem 2. Let v and v' be two steady flows in WV, subject to a prescribed 
velocity distribution on the boundary of VW. Let —m be a lower bound for the 
characteristic values of the deformation matrix of the motion v, let Vmax v, and 
suppose that either 

malySa or VdlyS Ya. (16) 


Then the two flows are identical. 


Proof. The kinetic energy of the difference motion v’—v must be constant, 
since the flows are steady. On the other hand, it must satisfy both (8) and (9). 
In view of (16) this can happen only if #=0, which in turn implies v =v’. 

This theorem depends strongly on the condition (16), but without some such 
assumption it is extremely unlikely that the conclusion is true. 


4. Variational techniques 

The key step in the proof of Theorem 1 lay in establishing inequality (11). 
We are concerned, however, not only with the validity of (11), but also with the 
determination of the largest possible coefficient «, for the size of « evidently 
determines the numerical values in the various stability criteria. Now the method 
of proof in Theorem 1 clearly gives no guarantee of providing the “best possible”’ 
value for «; moreover, it is quite crude in its estimate of the term us D-u. 
For these reasons, it is of great interest to consider an alternate approach which 
not only supplies a theoretical procedure for determining the best possible «, 
but in addition gives a way to avoid the estimate (14). 

Stated in precise terms, our problem is twofold. First, we must determine 
the greatest coefficient « such that the inequality 


ad fui<f gradu: gradu 
holds for arbitrary vector fields u satisfying 
div == 0,5 tj—=-0,0n 7% (17) 
Second, we must determine the least coefficient » such that the inequality 
f(u-D-u+ gradu: gradu) = 0 


holds, again for arbitrary wu satisfying (17). In this case it is clear that the basic 
flow will be stable provided simply that »>%. The two problems above can be 
consolidated into the single variational problem: 


—fu-D-u = Maximum, (18) 
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where the vector field wu must satisfy the side conditions 
fgradu:gradu=1, divw=0, w=0on Fi (19) 


The first problem occurs when D is the negative of the identity matrix. 

The variational problem (18) —(19) can be reformulated as a partial differential 
equation for the extremal function u, according to well known procedures of the 
calculus of variations. Thus, through introduction of Lagrange multipliers »* 
and A=A(x, y, z, t), the problem becomes 


Of (ueD-u + * gradu: gradu — 2A divu) =0. (20) 
The Euler-Lagrange equation corresponding to (20) is easily found to be 
ue D ——graddA+7*V?u, (21) 


and this is to be solved subject to the side conditions (19). The reader may 
observe the remarkable similarity between the eigenvalue equation (21) and the 
equations of hydrodynamics. The equations corresponding to the case D=—TI 
are of sufficient importance to be noted separately, namely 


yv*V2u+tusegraddA, divu=0. (22) 


In spite of the relative simplicity of (22), we have been unable to determine its 
solution, and in what follows can only offer some general remarks concerning the 
system (19), (21). 


Now for any solution of equation (21) we have 
—fu-D-u=f (uegradA—v*u-V2u) =r* f gradu: gradu = 2*, 


where use has been made of the divergence theorem and the conditions (19). 
On the other hand, any vector u which provides the integral (18) with its maximum 
value must be a solution of (21). It follows that the eigenvalue % associated with 
a maximizing vector is precisely the maximum value of the integral (18). Moreover, 
no eigenvalue v* can be larger than 9, for if this were the case then the cor- 


responding eigenvector u* would give to the integral (18) a value larger than >. 
This proves the following theorem. 


Theorem 3. Suppose there exists a vector u which solves the variational problem 


(18) —(19). Then the eigenvalue problem (19), (21) has a greatest eigenvalue >, and 
the basic motion will be stable provided that y> >. 


The difficulty with Theorem 3 is, of course, that there is no direct way of 
verifying its hypothesis. In circumstances such as these the best that can be 
hoped for is to determine all the eigenvalues (by Theorem 1 they constitute a 
bounded set), and to assume that the least upper bound of these eigenvalues is 
the required maximum of the integral (4 8). If the eigenvectors are complete in 
the set of admissible vectors, then this gives another method for proving that 
the least upper bound of the eigenvalues is the maximum of (18); unfortunately 
this appears at least as difficult to verify as the hypothesis of Theorem 3. . 
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5. Example. Couette flow 


We consider the well known circulatory flow between rotating concentric 
cylinders. If the inner and outer cylinders, respectively, have radii R, and R,, 
and rotate with angular speeds 2, and Q, (Q,>0), then the velocity field is given 
b 

y Ue ee) 
where 


A= R Q,—Rj Q, == (R, RP)? (Q,—2,) 
Ri-R?’ R2—R? 


One finds easily that the scalar vorticity wm has the constant value 2A, while 
in polar coordinates 


The characteristic values of D are + B/r?, whence 


B 
u-Deuz — Fle (23) 
The special form of inequality (23) suggests that the stability estimates of § 2 
can be considerably improved. In particular, in place of the inequality (11) we 
shall look for an estimate of the form 


Bf = S [gradu :gradu. (24) 
To this end, let us seek a vector field h such that 
div ih? = C2 /77 = hk <7 = Re. (25) 


A radially symmetric solution of (25) is readily found, namely 


h = = tan(C logy + D)i,, (26) 
where the constants C and D are given by 
Ce: Dat log (R, Ra) 
log (R3/Fy) ° 2 log (R,/f,) 


Substituting (26) into (13) and setting e =0 then yields an inequality of the form 


(24), with 2 It 2 
p= C= leet IR 


Combining (23) and (24) with the fundamental energy equation (4) yields 
dk e 
< S(|B|-8») | 


and from this it follows that Couette flow is stable relative to arbitrary disturbances 
whenever | een? 


a 2 
estat eG Ri) te, wath) . a 


The reader should notice that this stability is in the “‘strong’’ sense, #—0 as 
t+ o. 
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Previously, SYNGE [13] has proved stability relative to infinitesimal disturb- 
ances whenever A has the same sign as 2,, which in the present case (Q,>0) 
gives the stability conditions w=2A=0, or equivalently, 2 2 (Ro/Ry)?Q 2) 
The relative zones of stability are shown in Fig. 1. For the celebrated pair of 
radii R,=3.55, R,=4.03, inequality (27) reduces to 


om < 10.92. (28) 


In Fig. 2 we have indicated the stability zones based on (28), on SYNGE S criterion, 
and on the calculations and experiments of G. I. TAyLor. If we suppose TAYLoR’S 


- x.°(R3-R?) 
[RiRe log(Ro/R,) 1? 


4 
-¥C 0 ve Qe 


Fig. 1. Stability zones for Couette flow. The zone w20 lies below the line w =0 


experimental data to be an accurate representation of the mathematical situation, 
it is seen that the right-hand side of (28) cannot possibly be greater than about 50, 
for otherwise there is an observed secondary motion. Our value is certainly not 


é ? = Q2/y 


Fig. 2. Stability zones for Couette flow, R,= 3.55 and R= 4.03 


very near 50, but considering the difficulty of the problem and the fact that (28) 
applies to arbitrary disturbances, one does feel that 10.92 is a quite respectable 
value. 


6. Couette flow (continued) 
In the previous section we studied the stability of Couette flow using the 
methods of § 2. In this section the same problem will be treated by the variational 
techniques of § 4. The end result will be seen to have a somewhat tentative 


Per. In this connection we may recall some earlier work of RAYLEIGH and SyNGE, 
in which the same criterion is established for inviscid fluids; cf. reference [14], § 4.2. 
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character, in line with the remarks at the close of § 4, though it still has con- 
siderable interest. For Couette flow the differential equation (21) and the side 
condition div u=0 take the following form (in polar coordinates) 


BO 2 Ode grr, 2 ov uU 
72 bias {4u re Oo aE 
Bian AO Aar ts. 2 ou v 
Va fo ole ve? 
aA (29) 
QS oie wi 
0z 
(ru) ala: an ow i) 
) vy oo Oz , 
- where 
sl fo) “oh iba, wae e 
— =S | 2 " 
A y OY” (r a 7? Oy 672? 


and u,v, w denote, respectively, the 7, # and z components of the perturbation 
velocity u. Equations (29), as they stand, are too difficult to solve in full 
generality. We shall therefore seek a particular solution of the form 


u=Uu(r)coskz, v=d(r)coskz, w=w(r)sinkz. (30) 


Eliminating / and w from the system (29), we get 


Agen leler 
(L— Ra =, 
(31) 
L—R)O Bie 
( il acre es 


where L denotes the differential operator 


1 d fe d Js 
y ay ( Fl 72 — 

(The derivation of (31) is most simply carried out as follows. From the last 
equation of (29) we have 


whence uw can be written in the form u=curlY, where 


Y= — (di, —Gi,) sink. 
Substituting w—=curl W directly into (21) and taking the curl of that equation 
yields 
Bcark(ol72;uj72, Oreo curl* ¥, (32) 


where the identity curl curl =grad div —V? has been used. The first equation 
of (31) is now obtained as the #-component of equation (32). Finally, the unkown 
AZ must necessarily be independent of #, so that the second equation of (29) 
reduces immediately to the second equation of (31).) 
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Equations (31) are formally similar to the classical first order equations for 
small disturbance in Couette flow. The solution of (31) is a difficult task unless 
the cylinders are close together, in which case we can suppose (31) to be suf- 
ficiently well approximated by the simpler equations 


(33) 


where D=djdr, A=B)y*R®, and R is some appropriate mean radius, say 
R=YR,R,. The boundary conditions associated with (33) are 


G20 = 20" At = Ry Be (34) 


According to the calculations of JEFFREYS and others (cf. [14], § 2.3) the critical 
value of A is 41.2 (R,—R,)~® (that is, this is the smallest value of A for which 
a non-trivial solution of (33), (34) exists). We are thus led to the following 


stability criterion, 


Pest cieeeai Nis Vil (35) 


a _ *, 
v R, R,(R,—R) ; 


For the pair of radii 3.55, 4.03 this gives, in particular, 


et | 45.5. (36) 


(36) is certainly a real improvement over the earlier estimate (28), and, moreover, 
the right-hand side comes very close to the “‘best value’”’ 50. It must be borne 
in mind, however, that (35) and (36) are not yet rigorously established, since it 
remains to be shown that we have really found the greatest eigenvalue of the 
system (29). There is some reason to believe that this is so, since hydrodynamical 
experiments have shown indisputably that the secondary motion occurring in 
Couette flow is approximately of the form (30), but barring a proof of this fact, 
we must accept (35)—(36) as only tentatively established, in contrast with the 
absolute certainty of (27)—(28). 


Note: The above research was supported by the United States Air Force Office of 
Scientific Research under Contract No. AF 49(638)-262. 
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Rayleigh’s Problem in Hydromagnetics « 
The Impulsive Motion of a Pole-piece 


GSeSALUDFORD 


Communicated by HILDA GEIRINGER 


Introduction 


Consider the flow of an inviscid fluid past the pole-piece of a permanent 
magnet, t.e., a perfectly conducting rigid body containing a magnetic field. ii 
the fluid is non-conducting, the field of the magnet exerts no force on it, and, 
as is usual in the ordinary theory of inviscid fluids, no condition on the tangential 
component of velocity at the pole-piece need be satisfied. At the other extreme, 
when the fluid is a perfect conductor also, its velocity, v, must vanish at any point 
of the pole-piece at which the normal component of the magnetic induction B 
is non-zero*. It follows that any initial discontinuity in velocity at the boundary 
must immediately move away**, in contrast to the case of a non-conducting 
fluid, where the discontinuity can cling to the pole-piece as a vortex sheet. 


The object of this paper is to give a simple exact solution for which the 
transition from zero to infinite conductivity can be traced and the modifying 
effects of viscosity determined. Thus we consider the motion of an incompressible, 
viscous, electrically conducting fluid contained between the parallel, plane pole- 
pieces, y =0, A, of a permanent magnet which provides a uniform external field 
in the y-direction. Starting at time ¢=0, with the fluid at rest, the magnet is 
made to move uniformly with velocity wu) in the negative x-direction. 

The limiting case of zero conductivity and infinite separation (k= oo) was 
considered by RAYLEIGH [/]; the fluid velocity distribution is then given by 


(1) u = uy erf (y/2 Vv?) 


where » is the kinematic viscosity and the axes move with the magnet. Recently, 
Rossow [3] has extended the discussion to a conducting fluid, but his equations 
are inexact (see Section 5). In the present paper the solution of the general 


problem is obtained by means of the Laplace transform; in certain cases it can 
be expressed by known functions. 


Baling tangential component of the electric field E is continuous, and hence 
ae at the interface, while E+ x B is zero everywhere in a perfectly conducting 
uid. 


** For the resolution of such shear flow discontinuities in a compressible perfectly 
conducting fluid see [6]. 
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Infinite separation of the pole-pieces is considered first. The boundary-layer 
solution for vanishingly small viscosity gives, at any time ¢, a velocity profile 
which is similar to (1), the similarity factor being an exponentially decreasing 
function of ¢. This transport of vorticity out of the boundary layer becomes 
more rapid as the conductivity increases, the vorticity tending to become con- 
centrated in a diffused Alfvén wave. This is borne out by the special case of 
equal, but not necessarily small, diffusion coefficients, for which the solution is 
given by known functions, and the case of small and comparable coefficients, 
when the principal effects can be determined by the method of steepest descent. 
Finally, we give general series expansions which are particularly suited to the 
initial stages of the motion, and also discuss the character of the flow for large 
values of ¢. 

The flow with finite separation of the pole-pieces is the sum of a series of 
flows, each of which is of the kind described above; the same result is obtained 
by the method of images. It immediately follows, in particular, that each of 
the two Alfvén waves which leave the pole-pieces when the viscosity and resistivity 
are small, is repeatedly reflected, first at one and then at the other of the two 
pole-pieces. 

Whenever the diffusion coefficients are not both zero, the magnetic disturbance 
dies out as {> co. For infinite separation, however, there is a residual magnetic 
field, its value being independent of these coefficients. This is of importance in 
calculating the shearing stress on a pole-piece, since the Maxwell stress contributes. 
Thus, the assumption that / may be replaced by oo for widely spaced pole-pieces 
leads to completely wrong estimates when ¢ is sufficiently large. 


At all stages of the motion the shearing stress is larger than 
(2) T = Oto Vo/xet, 


the value obtained in RAYLEIGH’S case, although the viscous part is smaller. 


1. Equations of Motion 

The fluid is assumed to be incompressible, viscous, and electrically conducting, 
and to fill the space between the pole-pieces y =0, h of a magnet which is supposed 
to be a perfect conductor. All quantities are functions of y and ¢ only, so that 
if the velocity v and magnetic field H are given by (u,v, 0) and (H, Hp, 0), 
respectively, we find*: 

(1) from the equation of continuity: div v =0v/éy=0, so that v must be 
zero since this is its value at either pole-piece; 

(i) from div H =0H,/dy =0, it follows that Hy has at each instant its value 
at either pole-piece, which, by continuity of normal induction, is constant; 

(ii7) since the current J is given by curl H it follows from the conduction 
equation: 
(3) J=o(E+p0xH), 


that the electric field E acts in the z-direction: E =(0, 0, £). 


x Displacement currents will be neglected and the material coefficients y (viscosity), 
uw (permeability), and o (conductivity) assumed constant. 
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The equations governing the motion are therefore the x-component of 


MAXwELL’s equation curl E = —  0H/0t, v2. 
jichalen go 2 

(4) By ba 

the z-component of (3), = 

(5) anew Cane 


and the x- and y-components of the momentum equation, 


ou oH au 

(6) o ap Ht Ho ay + 00” Bra 
ow oH 

(7) O ae Oy” WH ay? 


where @ is the pressure and @, the density. The remaining components of these 
equations are satisfied identically. 

Equation (7) integrates immediately to show that @ + «H?/2 is constant at 
each instant, while on eliminating E from (4), (5), and (6) we find that H and u 


tisf 
pee @ 2) 1 yo 
(n ay? a aren Ey ‘ 


(8) a oH 


— + Holy =. ea. “-)u =0 
ay? ot : 


where Ay= |/uHé/oy is the Alfvén velocity and 7 =1/uo is the magnetic dif- 
fusivity. The initial and boundary conditions on the solution of this system are* 


(9) t=0: HO t= for alley, 
(10) SPOR Y 57 =H =0 for all ¢, 


the latter being equivalent to the vanishing of E and wu, see (5). Since the magnet 
is a perfect conductor there is no restriction on H itself at a pole-piece; any 
discontinuity in the tangential field corresponds to surface currents (in the 
z-direction). 


2. Laplace Transform of Solution 
From (8) and (9) the Laplace transforms of H and w satisfy the equations 


qd? d 
(7 dy? p)H Hy 7, =9, 


(11) ad ait 


7 + Hol? qt —P)“=—H, Uo, 


where # is the transform parameter and the same symbols are used for the trans- 
forms. The general solution of this system is 


Hi Ae) Be-"Y 1 C ety ey) 


Wo Uo p—nm [A emmy Be-™] cI of temas a he De”), 


Hym Hyn 


* As in the Introduction, the axes move with the magnet. 
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where A, B, C, and D are arbitrary constants, and y= -+m, + are the roots of 


(12) (nr? — p) wr — p) — Afr =o. 
We take 7 
m=Vatbp+Vatch, n—=Va+bp Va+ch, 
where 
(Spel ia we peel ae) a atid (Le Adley 
47 v 4n 4nv 


The constants A, B,C, D must now be chosen so that H and u satisfy the 
boundary conditions (10). Thus we find 


Ayu mn sinh m(y—th) sinh n(y—+th) 
13a H= OTe F. ce 3 ES Pours Ree rs 
A) p? m*—n2 cosh mh ie cosh nh F 

se, UA, Uo 
u x 
13) Peep nrin?) 


cosh m(y — $h) 


m2 (p — nn?) 


x [n2(6 — nm?) ee aha | 


cosh 4m h coshinh — 

or, for the case h = on, 

A = —n 
(14a) Hee So eee ne-™Y 1 me"), 

U U —m —n 
(14b) Vy a . Pin nd) [n? (pb — m*) e—™ — m2 (p — nn) ec"). 
A somewhat simpler function is the current / =—dH/dy: 
(15) p= yuo 1 a cosh n(y—th) | 

4yv V(a+0 p) (a+c p) cosh mh cosh $x h : 
or, for / infinite, 
(16) [== Hy Uo : e Verh sinh Va+cpy. 


2» V(a+b p) (atc p) 


3. Infinite Separation and the Special Case n =v 


In the next four sections we shall be concerned with infinite separation of 
the pole-pieces: 4 = oo. Even then none of the transforms H, wu, or J can, in 
general, be inverted in terms of known functions. However, when 7 =v =k (say), 
1.e. C=O, this can be done; we find 


ee ee er eee at 
3 i= 1 Hots erf (222) Le ¢ (2st 
Fieger 2VRt aki ] 
Cia greased eae) 


These functions are graphed in Figs. 1—3. 
Arch, Rational Mech. Anal., Vol. 3 2 


18 Gus. Ss, LuprORD: 


0125 


kT [Hp Uo 


aS ‘i ~/ocus of maxing 


% Y 2 J 
Agy/4k 
Fig. 1. Variation of current J with distance y for various values of T=A2 t/4k, when y=v=k 
1.0 


HAg / Ng uo 
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Agy/4k 
Fig. 2. Variation of the disturbance field H with distance y for various values of T=A3 t/4k, when v=)=k 
/=1/6% 


Aoy/4k 


Fig. 3. Variation of the velocity w with distance y for various values of T=A3t/4k, when v=yn=k 
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Clearly J is the total effect of two diffusion sources of strength + Hym/Ay 
moving with velocities +A,, the diffusion coefficient being k. In the limit 
k->0 it becomes (for y>0) Hou /Ay times the delta function 6(y— Apt), and 
then H and uw are Hyu,/A, and um, times the Heaviside unit functions H(A, ¢— y) 
and #(y—Apot) respectively. The initial discontinuity moves away from the 
pole-pieces as an Alfvén wave which is both a current sheet and vortex sheet: 
H jumps from zero to Hyu)/A, while w jumps from u, to zero. It is a mistake to 
speak of this wave as a shock, however (cf. STEWARTSON [4] who obtains a 
similar limit in another problem), since it is a valid solution of (8), for 7 =» =0, in 
which the variables have finite jumps across the characteristic y= A,?*. 

It is easy to see why the equality of 7 and y leads to simple results. The 
operator of the system (8) is then the product of the two operators 


and these are diffusion operators relative to axes moving with velocity +A. 


When )7—J/y is small, but |/7y is not, successive approximations to the 
solution can be found by regular perturbation. Thus in (16), for example, we 
expand in c: 


exp {I/atbp+Vatcply} _ exp{lVa+op + Val y} r 
V(a+0 Pp) (a+c p) Va(at op) 


where K = — (1 Vay)/2a, and hence find (to the first order) 


basics (( a hs a) layane oP | 
(+8 a) apt aad 


eee a 


2\Vakt Akt 
where 
farts 4ynv es 4nv igGxstee (Va—Vv)? l= ae |. 
ote” (n+yy?’ 2AR 240? 


The motion is therefore effectively the sum of a pair of diffused waves with 
reduced velocities --Ag and diffusion coefficient k, whose square root is the 
harmonic mean of the square roots of 7 and ». 


4. General Case with y and v Small but Comparable 


When 7 and y are small, but otherwise quite general, an indication of the 
character of the motion can be obtained by the method of steepest decent. Setting 
p =A§s/|/nv in the Bromwich contour integral defining the inverse of (16), we 
have to evaluate 

y+too 


| ds 
(18) es exp 0 jst — y{V1 +as+V1+ 8 s}/2A,| lGeeeaHaeES| 


* See the remarks of von Miszxs [5] on the similar question of contact discon- 
tinuities in ordinary hydrodynamics. 


QF 
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for large values of # = Aj/|/nv. Here . 
et ce i i ee = ate | 
6 \2 ap Ve i Diy p= \74 te , a B 


and y is any real constant > — 4/«. 
Now the stationary points of the arguments of the exponential are deter- 


mined by 
a B wot 
ee 
Vites V1+Bs y 
and with s the corresponding functions of x the largest argument is obtained 
for ~—1 with the lower sign. Hence, for x sufficiently close to unity, the dif- 
ference of the integrals (18) has the asymptotic form 


ull. B(1tas)i(1+6s) pr exper ls — (1 Fas Vi Rg 
[a2(1 +B s)’—B2(1+a5)*]? 
where s=2(1—x)/(2 +B) +O(1 — x), 


and elsewhere is negligible. | 
It follows that J has the approximate asymptotic form | 
Hyu 1 (y—A gt)? 
4 oe ee eX ieee \, 
a) J Ay V2x(y+») ¢ : 2(n+%) ¢ | 
which should be compared with the exact result (17) for 7 =v, and interpreted 
| 


in a similar fashion; the effective diffusion coefficient is the arithmetic mean 
of 7 and ». 


5. The (Viscous) Boundary Layer and the Flow Outside 


Clearly the above arguments do not apply when 7 and y are not comparable, 
since then « and f are also large. Of particular interest in this respect is the | 
case where y is small compared to 7. | 

For any fixed 7-0, the boundary-layer solution at the pole-piece is obtained 
by letting »>0 with y/J/» fixed. In the limit we find 


(20) Vo ah Vp ay tn a | 
so that VuP nia’ | 
(21) jm Hate edith ent (2) 

di 2 \» t 


while from (14) 


(22) Eh = oe erf (4, / “|, U = tye 434" erf eae i 

0 fT 2) vt 
Across the layer H is constant at any instant, but increases with ¢ from 0 to 
Hyuy/Ag. On the other hand J increases with y from 0 to (Hyuo/n) exp(—Ae/n), 
the latter value decaying from Hyuy/7 to zero as ¢ increases. The velocity prow 


at any instant is similar to RAYLEIGH’s [see (1)], the scale factor exp (—A$t/n) 
decreasing from 1 to 0. 


*~ Note that J =H,u,.se that B= 0 [see (5) ]. 
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Besides this similarity property of the velocity profiles there are two other 
interesting points to notice. First, the same result for u is obtained by assuming 
that the disturbance of the electro-magnetic field is small [H =F =0 in equa- 
tions (4) and (5)], but that the gradient of H, i.e. the current J, is important 
for the velocity field [J =— @H/éy in (6) replaced by —wH,/n from (5)]. This 
is the argument used by Rossow [3] in discussing the present problem; what 
he in fact obtains is the boundary-layer equation for uw, though under the un- 
necessary assumption that H/H, is small. However, this does not give H, which 
can only be obtained by determining the flow outside the boundary layer. His 
calculation of the total shearing stress on a pole-piece (cf. Section 9) is wrong 
since it involves using this boundary-layer approximation outside the layer. 
Secondly, the assumption that the magnetic field is not disturbed at the pole- 
piece is incorrect in O (1/7) for 7 large. 

Owing to the presence of the exponential factor in uv, the boundary layer 
effectively disappears in a time of the order of 7/A§. For a poor conductor this 
time may be quite long, and there will be a boundary layer of the same type 
as RAYLEIGH’S (7 =o) at the pole-piece. On the other hand, for small 7, the 
boundary layer will barely form, and the initial discontinuity may be expected 
to be propagated away as a diffused Alfvén wave*. This we now show by 
discussing the inviscid flow outside the boundary layer. 


This time we keep both y and y fixed as we let »->0, and obtain 


— Hy | Cul Pee Wotly kee 1 —ny 


Ay p\/1+n p/AR p 1+npj/4g 
fie Ele toi foe 
Aj 1+ p/AG’ 
where 7 is given in (20). It is easily checked that this is in fact the solution of 
the system (8), for »=0, which satisfies the boundary condition E =0 (but not 
adH/dy=0) on y=0O. Moreover the values of H, u, and J for y=0O agree with 
the corresponding values in (21) and (22) for y= oo. 

As in the case of general 7, » the inversion cannot be performed in terms of 
known functions. However, as there, the behavior of the inverses for small 7 
can be determined by the method of steepest descent. The change of variable 
in the Bromwich integral is now ~=A§s/y and the large parameter becomes 
9 —=A§j/n. We omit the details, since they are now quite similar to those in 
Section 4, and only state the approximate result: 


iy Uo 1 : oa 


which is the same as that obtained from (19) by formally setting y=0. The 
effective diffusion coefficient is $7. 


6. Small and Large Values of t. The Residual Field 
For unequal values of 7 and y the flow must be computed from series expan- 
sions, if more than its asymptotic character for small 7 and y is required. In 
obtaining these expansions we follow the ideas of GOLDSTEIN [2]. Thus for 


* When the vanishing of 7 keeps pace with that of », as in Section 3, nothing that 
can be identified as a boundary layer even forms. 
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small ¢, we treat p as large in equation (4 6) and expand in terms of P=? +a/lb: ‘ 


exp[—Vatop+Vateply 


Ba b—c) 1 a? (b—2)"2 nl ee 
= exp [ (Vo ae Ve) VP | “€Xp oe 7 : \P sb2ci =P 1 ly 


- a a(b—c (Das) ye arh 
= exp, (lo = Ve) VP y] {1 = oe re | 3 ies iM 


exp[—Ja+b (Di a Va+cply __ ©xp [= (Vo Vc) \yP y] x 
V(a+0 p) (a+ep) JbeP nes : 
a(b—¢) 9 pig) b= 9): (e(0— 8) ety Ww te wee 
x {1 a 2b Ve ye : 2bce 4b y 1) Ve + } 
Hence we find that, when 7 +», 


lig cmedeenr eri 
(23) x emer D(A) + Alan Dae) + B(2t) Ds) : a3 
ae Pal) cad A(2t)*D_, (Tea) + B (22) Des Gea p< ~f} , 
where Aby a 2A nv Az y? | 
~ nV) Wn FV)? ’ (7—»)? | 4yn» (Vii? ; 


and D, is the parabolic cylinder function of order ». By examining the con- | 
vergence of the two expansions involved we see that for computational purposes 
A \nvt|(y —)? must be small, while at the same time (7 —|) y//vt should — 
not be large. | 

The corresponding result for »=0, 7.e. for flow outside the boundary layer, 
can be obtained in the same way, or as the limit of the above (y fixed): 


Beall othe Abt yh wied (asbanae ( y 
a 4 ex 12 ed! pee ie 
i \2 n Xp n is Bat Loa a 7 n—1 Vant 


This expansion can be used for computation whenever 4jy /¢/7! is small. 


Each of the functions H — Hyu)/Ayp, u, and J has an infinity of order one 
half at = —a/b in the p-plane, and no other singularity (the one at p= —ajc 
is removable). Hence each behaves like ¢~? exp [— AG ¢/(J/7 + Vv)?] as t oo with 
y fixed. More particularly we find 


(24a) H mw gts 4 — Wane os F Agy 
o VrAy Vn» (Vn+|») 


ae l- aap 


a ey Oe eee ~~ 
(24b) “° Vn(Vnt)>) Ho %o 


, Ayy | 1 Agt 

Cee ee ID r =f | § exp |= =| 
Vz4y Vn» Vnv( Vntyr) ! \t (Vn +)? 

The most important thing to notice is that there is a residual field H =H o%/Ao 


and hence a residual shearing stress on the pole-piece. We shall return to this 
point in Sections 8 and 9. 
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7. Finite Separation of the Pole-pieces 


If in the previous analysis the pole-piece is located at y= vy, instead of at 
y =0, then y must be everywhere replaced by y—y,. If in addition the fluid 
lies below the pole-piece (y< yg) instead of above, then y must be replaced by 
Yo — y and the sign of H changed. We shall now show that, with finite separation 
of the pole-pieces, the flow is the superposition of a series of flows, each of which 
is of one of these two kinds. 


Consider, for example, J as given by equation (15). Since 


cosh m (y—+h) em(y—3h) EE ey —3 h) ey —h) ae git 
coshimh gimh 4 —hmh 7 1ten™ 
[o-e) 
= (eae ert) (4 yh th 
j=0 
CO 


co 
=) (=A Per? isis le») (Sid ener i) 
j=0 i=0 


and similarly for cosh (y—$h)/cosh $nh, it follows that 


(25) p=¥ (—1)) £-GAh) a (<1) $*(—Gh); 


here Y* (yo) and Y~ (yo) are the current functions of the two flows referred to 
above, respectively. Similar results are found for H and uw. 

The same conclusions would have been reached by the method of images; 
JZ (7 + 2h) is the reflection, in y =h, of — f*(—7h), which is itself the reflection 
nee? (72) in 4 = 0. 

The results obtained for infinite separation therefore yield immediate descrip- 
tions of the flow for finite separation. Thus, when 7 =v (see Section 3), J is 
total effect in 0< y<h of diffusion sources + H)u%/A y moving with velocity A, 
away from the points y=0, +2h, +4h,... and y=+h, +3h, ..., respectively, 
together with their negatives moving with velocity — A, away from the same 
points. For 7 and y small and comparable (Section 4) the motion consists essen- 
tially of two diffused Alfvén waves, one starting at y =O and suffering successive 
reflections at y=h and y=0, and the other originating at y= and moving 
exactly oppositely to the first. In the limit of zero 7 and » the velocity is always 
zero in the regions between the pole-pieces and the wave fronts, and alternately 
+u, between the fronts. On the other hand, H is always zero between the 
fronts, and is alternately +H )u)/A) and -— Hyu,/A, in the lower and upper of 
the aforementioned regions, respectively. The wave fronts themselves are moving 
current and vortex sheets. 

As y>0 with 7 fixed, boundary layers form at the pole-pieces. The values 
of J and w in the lower boundary layer are again given by equations (21) and 
(22) of Section 5, since the correction terms for the presence of the second pole- 
piece are constructed so as to give zero contributions to J and w at the first. 
However, the correction terms for H do not sum to zero and consequently, 
although 4 is still constant throughout the boundary layer, its dependence on 
time is not given by (22) but by the function whose transform is 


H = —2o%o— tanh{ De i 


pVnp +48 2\np +48 
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For small ¢ we obtain the H of (22), but for t-> co this H tends to zero instead 
of Hyu/Ay. The corresponding formulas for the upper layer are obtained by 
replacing y with #—y and changing the sign of H. Outside the boundary layers 
the motion can be described, at least for small 7, as in the general case above. 


8. The Initial and Final Stages of the Motion 


Up to values of ¢ comparable with the smallest of h?/7, h?/v, and h/Ay, only 
the first few terms in each of the series in (25) will contribute to J for O<y<A. 
In fact, for O0< y<#h, the first term of the second series will predominate, while, 
for 4h<y<h, the same term in the first series will. Any term retained may then 
be evaluated by means of (23) if 7=-y, and (17) if 7=yv. Thus the initial stage 
of the motion is computed by correction of the solution for infinite separation. 

For large values of ¢ we apply the partial fraction rule. First note that 
p =—alb is now a removable singularity as well as =—a/c. For if in equa- 
tion (15), J is considered to be a function of Va+bp and a+ cf, it is easily 
seen that it is an even function of each. Thus the only singularities are poles 
at the points determined by the vanishing of cosh $mh and cosh $nh, 1.e. by 


m == (27 +1) a and n= (2 +4)", 
where 7 is an integer. According to (12) these points are 
; / PED) 
muiihis otk patna ec treet a is 4ARh? |. 
(26) iE P; 2h2 Ui ipa ee | (7 y) (27+1)2 m2 |? 


clearly only non-negative values of 7 need be considered. 
In general these poles are simple, the residue at #7 being 
Ayu | 1 _, Cos [(2j+1) a(y—sh)/h] 
4nv  V(a+b pF) (a+c p;) t(—1) [5h (dm/dP) ps) 
= — 4H, uy sin [(2j-+1) « y/h] 
V (25-41)? 2 (n—»)?— 4.45 1? 

and that at #; the negative of this value. Hence, on combining terms in pairs 
and assuming h<2|7—v|/2A,, we find* 


(27) J=8H uU S sin [(2j-+4) w y/h] e“i* sinh B.¢ 
h y pa V(2i+1)? 22 (n—v)2—4 AB hy 
where 
Sey ays 2 m2 (n+ 27 —— = 
a= (2) +1) B= OE" Veazey 0 — 4Aee, 


When h2a|n—v|/2Ay, f; is pure imaginary or zero for each value of 7 such 
that 27 +15 2A ,h/x|n—v| and the corresponding sinh-term in (27) changes into 
a sine-term or a ¢-term. Similarly, or by use of (8), we find 


0% 


8h 
Tie al colar lage] -e~%' sinh B,t, 


JB (27-+1) VQF+1)? 2 —v)2—4 AB 


u—4u S sin[(27 +1) ay/h Laating (4—») sinh B;¢ , Cosh Bit | ast 
od, i Wieeeiaeiccceesi | Gears ; 


* Of course, separation of variables leads to the same result. 
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However, when =», the dominant term in these series for large ¢ is not 
necessarily the first. For consider p*, as defined in (26), to be a function of 
the continuous variable 
pera) 2S 

h? : 
Then we find that as x decreases, p; increases through real values to — a(1/b-+1/c), 
when x=a(b—c)?/bc, and then becomes imaginary with real part increasing 
further to zero when x=0. But #7; increases to a maximum of —a/c for 
x =a(b—c)/c before decreasing to —a(1/b+1/c) at x =a(b —c)?/bc, after which 
it takes values conjugate to #7. So long as this maximum point cannot be 
exceeded for non-negative integer 7, 7.e. so long as 


— Oe OOn 


—r Ay ’ 
the term 7 =0 dominates in (27). When 
2Va9Vn=Vol 4 < 2ln-9l 
72 Pe Aad” 


however, one or other of the terms whose 7’s straddle the value 


EEE i} 

2 Un Vn>|Vn—|>| 

is the most important. As / increases further this term will remain dominant 
at first, and then at some value before the range 


(28) po» eVatelva—VeL 
V2 4, 
it will give way to the first term once more (which is now oscillatory). 
In particular we see that there is no residual field H, as there was with infinite 
separation (Section 6). In fact, assuming (28) holds, we have 


ke 8h up 
a 4A h?— 2? (n—»)? 


= a 5 Ay Uo : cos ( Tt = sin (220 t| 
It Ala h h Fj 


H~ 4H, 


RY e~ (nt) 20 gin (2 4.43 h?— 2? (n—0)? é| 


ae ( h 2h? 


when / is very large. This result should be compared with (24a). 


9. Shearing Stress at a Pole-piece 


From the right-hand side of (6) or directly from its definition as the sum of 
viscous and Maxwell stresses, the total shearing stress is seen to be 


tT =UHyoH + oor Ouldy. 
Hence, according to equations (14), the transform of its value on y =0 is given by 


A = a : ae [Ap +»(p +nmn)| 
0 


— Mo , : ido (1 \/z)- heaton 
2\nv platbp 2 qn) Va+bp’ 
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of which the inverse is a 


Ss ae ert( a a \+ V2 [- qeimel: 


0 


For small values of ¢ this gives 


vu AG(2/7+/) o(t 
=, + (ny)? rage 


(29) 2o 


af 


70 ; v/s 
= 
S 
S 
S 
LY 
S 
locus of minima ~ 
as} = 
| Rayleighs case 
| 
| | co 
| 
| 
| 
| 
i. Q 7 a J 4 
Aptly 


Fig. 4. The total stress t on a pole-piece as a function of time, for various values of 7 


while for large ¢ we find 

Ait cee ue 
(30) =~ toy — exp | — pee | [to a + OH], 
in agreement with equations (24). The first term in (29) is the value (2) found by 
RAYLEIGH [1] for the purely viscous case (7 = oo); the first term in (30) is the 
residual Maxwell stress arising from the residual magnetic field (see end of 
Section 6). Hence t decreases from infinity more slowly* than in RAYLEIGH’S 
case, reaches a minimum [at ¢=|/»(]/7 +)/»)?/2A4$ /y], and then tends asymp- 
totically to 09% Ao, see Fig. 4. 


When the pole-pieces have finite separation the corresponding stress (at each 
of them) is given by 


T mn U 
=o x 
20 a ae 


x [— 2 {46 + »(p — 7 m?)} tanh 4 gmh + m{Ap+(p — yn n*)} tanh dnhj. 


* Note, however, that the viscous part of the stress is always less than it is for 
n = 00 by a factor exp [— 46 ¢/(/m +)r)?]. 
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This cannot be inverted in terms of known functions. However for sufficiently 
small values of ¢, depending on h, the inverse will differ little from (29), see 
Section 8. On the other hand, for large values of ¢ we apply the partial fraction 
rule as in the preceding section and obtain 


T | 4Uq | 2h? AS+ 2? v(n—») ye ee ie ‘ ’ 
Qo ah V4 ARh2— 22 (n— y)2 2h? . 


He i 4A oh? — 7 (n weit a ewe (n+y) t/2h? 
2h? 


+ a11005( 


= es Uy Ay sin (3 Q tle Diao ead Ute 
LG 


where it has been assumed that / is very large. Thus we see that, however large 
h is, t is eventually oscillatory, the period tending to infinity and the damping 
to zero as h tends to infinity*. 

Thus, for / large, the assumption = oo leads to incorrect estimates for both 
H and 1 for sufficiently large values of ¢t. In other words the limit operations 


are not interchangeable: 


lim lim + lim lm, 
t—>oo h->oo h—>oo too 


for any values of 7 and ». 
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Feblerabschatzung bei linearen Gleichungssystemen 
mit dem Brouwerschen Fix punktsatz 


JOHANN SCHRODER 


Vorgelegt von L. COLLATZ 


Mit Hilfe des Brouwerschen Fixpunktsatzes ({J], [7]) werden Fehlerab- 
schatzungen fiir lineare Gleichungssysteme Gu =r hergeleitet. Alle vorkommen- 
den GroBen seien reell. Die gegebene Gleichung wird in der Form u= Tu 
geschrieben und eine Menge ® ermittelt, welche der Operator 7 in sich abbildet. 
Diese Menge enthalt dann nach dem Brouwerschen Fixpunktsatz einen Fixpunkt 
u*—Tu*, d.h. eine Lésung u* des gegebenen Gleichungssystems. 2D besteht 
hier aus allen Vektoren wu, fiir welche x u<y bei geeignet gewahlten Vektoren 
ey gilt. 

Die sich auf diese einfache Weise ergebende allgemeine EinschheBungsaussage 
fiir die Lésung w* wird fiir mehrere Sonderfalle formuliert und auBerdem in 
verschiedener Weise umgeformt. Dabei stellt es sich heraus, daB eine Reihe von 
bekannten, bisher getrennt bewiesenen (z.T. aber noch nicht speziell auf Glei- 
chungssysteme angewendeten) Fehlerabschatzungen in dem allgemeinen Ergebnis 
enthalten sind. Dazu gehéren Fehlerabschatzungen, welche sich mit dem soge- 
nannten Fixpunktsatz fiir kontrahierende Abbildungen ({14], [4] S. 35 ff.) be- 
weisen lassen, und solche, die man aus [13] folgern kann, sowie die in [10] ange- 
gebenen Abschatzungen; ferner die bei monoton wachsendem oder fallendem 
Operator médgliche EinschlieBung der Lésung zwischen die Naherungen eines 
Iterationsverfahrens ([6], [2], [8], [11], [5]) und die EinschlieBungsaussagen bei 
Gleichungssystemen monotoner Art [3]. Verschiedene dieser Abschatzungen 
werden verbessert und verallgemeinert, die Klassen von Gleichungssystemen, fiir 
welche solche Abschatzungen gelten, erweitert. Es miiBte noch untersucht 
werden, ob man auf diese Weise wesentlich neue Typen von Gleichungssystemen 
behandeln kann. — Die Ergebnisse sind in § 7 zusammengefaBt. 

Die Resultate kann man ohne groBe Schwierigkeiten auf lineare Gleichungen 
in einem Banachschen Raum iibertragen, indem man den Schauderschen Fix- 
punktsatz benutzt. Damit sind dann z.B. Anwendungen auf Differentialglei- 
chungen méglich. Derartige Anwendungen sollen jedoch zusammen mit einer 


funktionalanalytischen Behandlung nichtlinearer Probleme erst in einer folgenden 
Arbeit besprochen werden. 


§ 1. Bezeichnungen 


A, B,... bedeuten m-reihige quadratische Matrizen mit den Elementen Ajh, 
Oar und 7,...,u4,v,... m-dimensionale Vektoren mit den Komponenten 


4 4 
iota Ory Cee OO OL 
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In der Menge dieser Vektoren wird eine Beziehung u<v in folgender Weise 
definiert : 


u<v bedeutet Yh;,u*< dh, v0" CES wise (ied) 
k=1 k=1 


Dabei ist H =(h;,) eine gegebene regulare Matrix. Praktisch wichtige Beispiele 
sind: 


1. u<v bedeutet w'< v’' fiir alle 1 =1,2,...,m: (122) 


2. u<v bedeutet u'<v' fiir bestimmte Nummern 7 =14,,%5,..., 4, (1.3) 
ui = vy’ sonst. 


(Wenn die Aufteilung (2.2) so gewahlt wird, daB die Ergebnisse von der Nume- 
rierung der Unbekannten unabhangig sind, kann man im Falle (1.3) etwa 


1; =] (je AeRe et) (1.4) 
annehmen.) 
Jeder Vektor w laBt sich dann in eine Differenz 


u=u*— u- (4.5) 


von Vektoren u* und u~ mit folgenden Eigenschaften aufspalten: 
a) Es ist oS ut und uSu* (dh. uw Zo). (1.6) 
b) Aus oSv, uSv folgt u*Sv. (17) 
|w| sei der Vektor 
|u| =ut+tu. (1.8) 


Zum Beispiel erhalt man im Falle 


(1.2): || = (|e), 
(1.3): w= | elie qtr reas) ea =e? 91t) 


—|u'| sonst, 
(14); || =H (|3 hie). 


Von dieser Ordnungsdefinition fiir Vektoren ausgehend, definieren wir eine 
solche fiir Matrizen: 


‘ 


A<B bedeutet AuxBu fiir jeden Vektor uw=Zo. 
Der Ungleichung A<B entspricht bei der Ordnungsdefinition 


Dye ia Ss Us, fiir alle 7, k; 
(1.3): 4,35, fir 1=1,, R=i, und t+ 1,, R-%, (x, 8 =1,2,...,2) 
A; = Os, sonst ; 


k= 
(14): &,<b;, fir alle i,k 
bel 


~ oy 


A=(a;,) =HAH>, B=(b;,) =HBH2. (1.9) 
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Formal genau so wie u*, w,|w| werden dann 4°, A- und | A| erklart. Man 
hat in (4.5) bis (1.8) nur 0, w, v durch O, A, B zu ersetzen. Zum Beispiel erhalt 


man im Falle 


(1.2): afg=4(a;,| +42), az, =4 (| 42] — 4%) fiir alle 7,k; 
(1.3): “ate =# ((ae,| + ain), ain =4 (| 4:4] — 4) 
fir i=%, kot, und it, bee, GP S120 
Gip—s\— |e a; x), Gi, = (—| Fel d;,) sonst; 


(44): AtT=HAAMH, A-=HAAOH 
mit : 
ajy re a (|a;.| + G;,) , asp -s 3 (|; | ae a; ») fiir alle 7, 
und (4.9). 
§2. Allgemeine EinschlieBungsaussagen 
und Fehlerabschatzungen fiir Iterationsverfahren 
2.1. Herleitung der allgemeinen EinschlieBungsaussagen. Die Matrix G des 
gegebenen Gleichungssystems 
Gu=r (e229) 
werde aufgespalten in eine Differenz 
G=A—B (2.2) 
mit einer regularen Matrix A. Die Gleichung w= Tu mit dem durch 
Tu=Mu+ts, M=A1B, s=Aty 
definierten Operator ist dann dem Ausgangsproblem (2.1) aquivalent. Wir 
suchen also einen Fixpunkt der durch T vermittelten Abbildung. 
® sei die Menge der Vektoren wu, fiir welche 
*Susy 
gilt, wobei x und y zwei feste Vektoren mit x<y bedeuten. Fiir w€®D schatzt 
man ab: 
Tu=Mu+s=M*tu—M-u+tsSMty—M-x+4+s (2.3) 
und entsprechend 
Tu=M*x—M-y+s. (2.4) 
Unter der Voraussetzung 
4—M*x+M-ySssSy—Mty+ M-x 


gilt also xSTusy fir w€®, d.h. der stetige Operator T bildet die konvexe, 
abgeschlossene und beschrankte Menge ® in sich ab. Diese Menge enthalt des- 
halb nach dem Brouwerschen Fixpunktsatz ((Z], [7]) einen Fixpunkt u* = Tu* 
der Abbildung 7, d.h. eine Lésung w* des Gleichungssystems (2.4). Wir fassen 
das Ergebnis zusammen: 

Ergebnis: 

V oraussetzung : 


*—M*x+M-ySssSy—Mty+M-x, OS 


/\ 
— 
N 
ee 
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oder 

Pal ai My SU) lye Me (ney, lh iesly (2.5’) 
oder 

«sTy—M*(y—4), Tx+Mt(y—x)Sy, GES: (2.5”) 

Fehlerabschatzung : 

x= M*x—My4+ssu*¥SiMty—M-x4+sSy (2.6) 
oder 

xSTx—M (y—x)Su*¥=Ty+M-(y—x)Sy (2.6’) 
oder 

= Ty—M* ys) aS Tx My —x) Sy. (2.6) 


Die genaueren der jeweils angegebenen Schranken folgen mit u*=Tu* aus 
(2.3), (2.4). Gilt A420, so kann man M~ und M* durch die gréBeren Matrizen 
A+B bzw. A B* ersetzen, wie man an den Formulierungen (2.5’), (2.6’) bzw. 
2.5”); (2:6) erkennt. 


2.2. Spezialfall des monoton wachsenden Operators. Es sei M~=O, also 
M=O. Dann wachst der Operator T monoton (im schwachen Sinne), d.h. aus 
usv folgt Tu<Tv. Definiert man 


Xo = %, Mi SHSE Kyi = TX, Vui=Ty, (K=O, 1, 2,..+)), (2.7) 


so geht (2.5’) in (2.10) und (2.6’) in (2.11) tiber. Unter der Voraussetzung (2.10) 
gilt ferner (2.12), denn sei 


Koss tye Sh, STS YS SS Vo (2.8) 


fur n =p (24) — fir p~=—1 ist dies nach (2.14) richtig —, so folgt 


al Be ee gee) UN.) Yee, 1 Vy sd Vo, 


d.h. (2.8) fiir n =p +1 und damit alle n. 


Ergebnis: 
Voraussetzung : 
Ma== (0 (d.h. MZO) (2.9) 
und 
%ySVo, MSM, WSo- (2.40) 
Fehlerabschatzung : 
KySMySUSUYSVo (2.41) 
und 
Rs ge es SSS Sp. (2.12) 


Iterationsverfahren mit monoton wachsendem Operator in halbgeordneten 
Raumen behandelten L. V. KANToRovitTcH [6] und D. MoRGENSTERN [8]. Prak- 
tische Anwendungen des monotonen Verhaltens von Iterationsfolgen bei linearen 
Gleichungssystemen zur numerischen EinschlieBung der Lésung findet man bei 
L. Cortiatz [2]. In der obigen Formulierung stehen die Ergebnisse fiir Verfahren 
in halbgeordneten Raéumen auch in [//] und [5]. 
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2.3. Spezialfall des monoton fallenden Operators. Ist Vivo also MSO, 
so fallt der Operator T monoton: aus uSv folgt Tu = Tv; Mit der Bezeichnung 
(2.7) geht (2.5) in (2.14) sowie (2.6) in (2.15) iiber. Entsprechend wie bei 
monoton wachsendem Operator beweist man hier auBerdem (2.16) fiir die durch 
(2.7) erklarten Folgen x, und y,. 


Ergebnis: 
V oraussetzung : 
M* = 0 (d.h. MSO) (2:43) 
und 
XeSVo, MSM, USYo- (2.14) 
Fehlerabschatzung : 
XpSVySUSUMSV (2.15) 
und 
Xo YS Kes Vrs 500 ES HIRES G0 SS 45 Vo vor (2.16) 


Im Spezialfall yy =x, lautet die Voraussetzung (2.14) 
an Gee (2.17) 


und die Fehlerabschatzung (2.16) : 


Nga toa SSS Sees Ae (2.18) 


Ein entsprechendes Ergebnis kann man fiir den Fall x)= vy, formulieren. 

Die speziellen Aussagen (2.17), (2.18) fiir das abstrakte Iterationsverfahren 
mit monoton fallendem Operator wurden in [17] formuliert. Dort wurde auch 
ihre Anwendung auf lineare Gleichungssysteme kurz angedeutet. Die allgemeinen 
Aussagen (2.14), (2.16) findet man in [5] funktionalanalytisch bewiesen. 


2.4. Fehlerabschatzungen fiir eine Naherungslésung uy=T uy. Es sei up 
eine Naherungslésung fiir w* und u,=Tu,. Wir setzen 


X=U+v, y=utw 


mit Vektoren v, w, fiir welche v<w gilt. Dann geht das Ergebnis von Nr. 4 
uber in 


Ergebnis: 
Voraussetzung : 
v—Mtv+Mwsu—u4Ss0—Mtw+ Mo, vw. (2.19) 
Fehlerabschatzung : 
vSu* —uySw (2.20) 
und 
Mtv — M-wSu*—u<Mtw— Mov. (2.21) 
Im Spezialfall v = — w erhalt man als 
Ergebnis: 
Voraussetzung : 


| -- %|<w—|M|w, wo. (2°22) 
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Fehlerabschatzung : 
|u* — u.|<w (2.23) 
und 
|u*—u,|<|M|w. (2.24) 


Im Falle A+2O kann man in diesen Ergebnissen M~, M* und |M| durch 
die gré8eren Matrizen A1B-, A1B* bzw. A|B| ersetzen. 

Praktisch interessant diirften auch noch die Falle v=o und w=o sein, 
insbesondere dann, wenn M~=O bzw. M+=O, also z.B. A120 und B-=O 
bzw... 5b* =O ist, 

Das zweite der hier formulierten Ergebnisse (2.22), (2.23), (2.24) kann man 
auch aus [13], §1 folgern. Jedoch wurde dies bisher noch nirgends ausfiihrlich 
dargestellt. Die Satze aus [13] liefern die zusidtzliche Aussage, daB die durch 
U,14=Tu, (n=0, 1, 2,...) definierte Folge u, unter der Voraussetzung (2.22) 
gegen u* konvergiert, und sie sind auch auf komplexwertige Gleichungssysteme 
anwendbar. — Die mit (2.22) aus (2.24) folgende gegeniiber (2.24) grébere Ab- 
schatzung 
|u* — u,| Sw —|u,— u| (2.25) 
wurde fiir die Ordnungsdefinition (1.2) und Aufspaltungen der in den §§ 4 und 5 
behandelten Arten in [10] angegeben, allerdings unter einer zusatzlichen Kon- 
vergenzvoraussetzung. Mit der Schranke o =w—|u,— | erhalt man im Falle 
A1=0 die fiir (2.22) hinreichende Bedingung 


Ao.=|B|(o+ |%,— %l)- (2.26) 


Das Ergebnis (2.22), (2.23), (2.24) stellt gleichzeitig auch eine Verbesserung von 
Aussagen dar, welche man mit dem Fixpunktsatz fiir kontrahierende Abbildun- 
gen ((14], [4] S.35ff.) bei der Norm ||«||=Max bal (w'>0) beweisen kann. 
Vgl. hierzu die Nrn. 4.3 und 5.3. : 

Wahrend im Spezialfall v =— w der Betrag | w*—u,| abgeschatzt wird, besteht 
bei dem ersten der in dieser Nr. formulierten Ergebnisse die Méglichkeit, u*— u, 
nach oben und unten in Schranken verschiedener Art einzuschlieBen, z.B. also 
auch Aussagen wie “*— u, =o zu erhalten. 


§ 3. Vom Defekt ausgehende Fehlerabschatzungen 
3.1. Allgemeine EinschlieBungsaussagen. Es sei A120. In diesem Falle 
kann man zur Fehlerabschatzung von hinreichenden Bedingungen fiir die in den 
Ergebnissen des § 2 genannten Voraussetzungen ausgehen, mu dann aber auf 
die scharferen der jeweils angegebenen Abschatzungen verzichten. So folgt aus 
(3.2) wegen (3.1) die Voraussetzung (2.5) mit A+ B* statt M* und A7*B statt M. 


Ergebnis: 

V oraussetzung : A1>0 (3.1) 
und (A — BY) x+ Boy SrS(A— B*) y+ Bx, HSY (3.2) 
oder (statt (3.2)) 

Gx+B-(y—*)SrSGy—B(y—n), (3.3) 
LiSYys (3.4) 
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Fehlerabschdtzung : xsutsy. (3.5) 


(3.3) kann man als Bedingung fiir die Defekte d,=—Gx*x+rundd,=—Gy-+r 
deuten. 

L. Corratz [3] definierte Aufgaben Gu =r monotoner Art in halbgeordneten 
Raumen so, daB Gx<Gy bei sonst beliebigen Elementen x, y immer *Sy zur 
Folge haben soll. Diese Forderung ist bei linearen Gleichungssystemen mit 
G0 identisch. Gilt G1=0, so folgt wegen ry =Gu* aus GxSrSGy (d.h. (3.3) 
mit B-=0O) die Abschatzung (3.5). Diese Méglichkeit der Fehlerabschatzung 
bei Systemen monotoner Art ist in dem obigen Ergebnis enthalten: Man setze 
A=G und B=O. Dann ist auch B~-=O und (3.4) folgt aus (3.3) wegen G?= 
A1=0. — Im Ergebnis (3.3), (3.4), (3.5) braucht das gegebene Gleichungssystem 
auch fiir den Fall B-=O nicht von monotoner Art zu sein. 


3.2. Fehlerabschatzung fiir eine Naherungslésung uy. Es sei uv) eine Nahe- 
rungslésung, deren Fehler abgeschatzt werden soll und d der zu uy gehdrende 


Defekt 
d=—Gu+r 


der Gleichung (2.1). Dann folgt im Falle (3.1) aus (3.6) durch Multiplikation 
mit A die Voraussetzung (2.19) mit 41B* und A+B statt M* und M-. 


Ergebnis: 
Voraussetzung : 
An2O 
und 
(A— Bt)v+ Bwsds(A—Bt)w+ Bo, vsw (3.6) 
oder (statt (3.6)) 
Gu+ B (w—v)<d<Gw— B(w—v), vsw. (3.6’) 
Fehlerabschatzung : 
vSu*¥—uwSw. (3.7) 
Im wichtigen Spezialfall v=—w hat man das 
Ergebnis: 
Voraussetzung : 
Antes 0 8 
Pe (3-8) 
|4|<(A—|B|)w, wo. (3.9) 
Fehlerabschatzung : 
| u* — u| Sw. (3.10) 


§4. Aufspaltung wie beim Gesamtschrittverfahren 
4.1. Aufspaltung von G. Die Hauptdiagonalelemente der Matrix seien positiv : 


(0) (1 =1,2,...,m). (4.1) 


ve ee die Matrix G auf in die Summe einer Diagonalmatrix D, einer oberen 
reieckmatrix —C, und einer el fib 
1 unteren Dreieckmatrix —C, (C, und C, sollen 


f 
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auch in der Hauptdiagonalen nur Nullen enthalten): 


—C, 
Ge D Be nite aC aC 


=+C, 
und setzen dann 
G—A SB mit—A=]D und B=—G. 


Man erhalt Tu—D141Cu+D"s, das Iterationsverfahren u,.,=Tu, (n=0,1, 
2, ...) ist das sogenannte Verfahren in Gesamtschritten. 


4.2. Monotonie beim Gesamtschrittverfahren. w<v sei durch (1.2) oder 
(4.3), (1.4) definiert. Dann ist A742O wegen (4.1). Wir geben vier Falle von 
Vorzeichenverteilungen bei den g;, an (allgemein war (4.1) vorausgesetzt), fiir 
aie -5s—O oder 6'—O und damit dann auch M~ =O bzw. M*t=0O ist, der 


Operator T also monoton wachst bzw. fallt: 
SY 


12) ¢5= 0, fircalle z2,k mit 728: G 
Oe e250 ttt alle 2, ke G=[ + ], 


fee = O- (Ul k,n.) Und 2k 1-1... Mm Mt tsk 


£;r—= 0 sonst: 


AN = Om iit = 452). b and t,R=1-+4,...,m 


he: 
Poe 0 ssonst- ee are 
aie 
4 es 0) (1.2) 
Im Falle ; gilt Facer bei der Definition one (1.4) 
4 B= M0 (1.3), (1.4) 


Das gewohnliche Differenzenverfahren fiir Randwertaufgaben mit der ebenen 
oder raiumlichen Potentialgleichung fiihrt oft auf Gleichungssysteme, bei denen 
gleichzeitig die Bedingungen der Falle 1 und 4 erfiillt sind, wenn man die Unbe- 
kannten in geeigneter Weise numeriert. Die Vorzeichenverteilung 2 erhalt man 
z.B. bei den Clapeyronschen Gleichungen. 


4.3. Fehlerabschatzungen fiir Naherungslésungen uy und u,. Es sei u, eine 
Naherungslésung und 1, die daraus nach dem Gesamtschrittverfahren berechnete 


1 Die Vorzeichen +, — bedeuten, daB die an den betreffenden Stellen stehenden 
Elemente >0 bzw. <O sind; nur fiir die Hauptdiagonalelemente ist g;;>0 vor- 
ausgesetzt. 


3* 
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erste Naherung. w<v sei durch (1.2) definiert. Aus (2.22), (2.23), (2.24) und | 
(3.9), (3.10) erhalt man dann das 


Ergebnis: 
V oraussetzung : 
™m i : 
| u, — up| < w* — : YY eal @, w=0 as ee eres) 
Sit fy 
ki 
oder 
, m k ‘4 se k F ’ tei ; i 
Ci — > ginmtr = g,,0° — » lgi,| 2", wZz0 (7 = » Ay erry . 
= k=1 
k=1 i 
Fehlerabschatzung : 
| u*¥*— uy| Sw" (i =A, 2S). mM) 
und Me 
|ju**—_ui|< d Dae Lib w (C12, ne 
Sit SS 
hi 


wie (4, Dee on) (4.2) 
die folgenden Aussagen. 
Ergebnis: 
Voraussetzung: 
m 
MS mit pa yare (G4) 2) as ene 
Bi 
hi 
und ae 
= Wileb ons 
ae} ; ae 4. 
[ad — a Sree ines S (45) 
i-1_*' = 0 gesetzt fiir |w, — | =0 und y;=1). 
1— oy 
Fehlerabschatzung : 
|u**— u|<Se (¢ 32 wat, 000) 
we |u**— wil <u, ¢ Cee vine Ace cre (4.4) 


Ersetzt man in (4.3) und (4.4) alle 4; durch # =Max mu; — dabei ist dann 


f#<1 vorauszusetzen —, so ergeben sich Fehlerschranken fiir die |u**— ui], 
welche auch mit dem Fixpunktsatz fiir kontrahierende Abbildungen zu beweisen 
sind. — Hier erhalt man im Falle {<1 im allgemeinen scharfere Schranken, 


und man kommt auch im Falle # =1 zu Abschatzungen, wenn man | ui — uo| =0 
fiir w4;=1 erreichen kann. 


§5. Aufspaltung wie beim Einzelschrittverfahren 
5.1. Aufspaltung von G. Wir fordern wieder (4.1) und setzen jetzt 


GA “oR ‘mitt ep Ope pea e, 
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und den Bezeichnungen der Nr. 4.1. Dann wird Tu =(D—C,)1C,u+(D—C,)4, 
und das Iterationsverfahren w,,.,= Tu, (n=O, 1, 2, ...) ist das sogenannte Ver- 
fahren in Einzelschritten fiir das System (2.1). 


5.2. Monotonie beim Einzelschrittverfahren. w<v sei durch (4.2) oder (4.3) 

mit ; : ! 

1,=27—1 (FAL eet) (54) 

—erklart. Wir nennen Voraussetzungen, aus denen A1+=O folgt, und solche, unter 

denen M~=O oder M*+=O ist. Die Elemente der Matrix A1+=(f;,) lassen sich 
in geschlossener Form angeben. Es sei 


. Bian 8a, a2 Sard ++: Sapk ( Sik os 
En ebint in boop onan poses a ee PO 
: ( Anges P ) 811 S01 01 See 02 +++ Ekk SiiSkr p 

un 


[i Rlp= D  Gyey,eg,.-- 0,2), 


I> >>: >k 

wobei tiber alle ganzzahligen «; mit 1>0,>%,>---> summiert wird. Dann ist 

hi => aE eae kj, fir 1>k, Ms ae fi,x=O fir 1<k. 
A1=O gilt daher unter folgenden Bedingungen: Bei der Definition 

(22): Wipe fiir alle 7,k, 

Caran (4) "7 = 0 e fiir’ alle 2, k. 
Dafiir reicht hin bei der Definition 

Chea pO Gill tees, 

(1:3), 10o.1)2 (= Aye 0m fin > k. 

Fir folgende vier Vorzeichenverteilungen wachst oder fallt J monoton im 
Sinne einer der Definitionen (1.2) und (1.3), (5.1)?: 


wes + ofS 


aN 
D 
I 


1 
Im Falle | : gilt A420 und bei der Definition 


3 B= M*=0 


5.3. Fehlerabschatzungen fiir Naherungslésungen uy und wy. uy sei eine 
Naherungslésung und wu, die daraus mit dem Einzelschrittverfahren berechnete 
erste Naherung. «Sv sei durch (1.2) oder (1.3), (5.1) erklart. Von (2.22), (2.23), 
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(2.24) und (3.9), (3.10) ausgehend erhalt man das 


Ergebnis: 
Voraussetzung - 
| —u%|S(E—|(D—C,)4C,|)w, wo. (5.2) 
Hinreichende Bedingung (fiir (5.2)): 
(D —|C,|) |m — %| = (D—|C}) #, w =o (5.3) 
oder 
(DC, Oy ld | pre es Uy [Onl) xe eee (5.4) 
Fehlerabschatzung - 
|u* — u| Sw 
und 


|u* — u,| <|(D —C,)1C,| w. (5.5) 


In (5.2) und (5.5) kann man |(D—C,)+C,| durch die gréBeren Matrizen 
|(D —C,)+||C,| oder (D—|C,|)4|C,| ersetzen. DaB die Bedingung (5.3) hin- 
reicht, ergibt sich durch Multiplikation mit (D—|C,|) +20. (5.4) ist die Vor- 
aussetzung (3.9), welche wegen A-t= (D — C,) 1&0 (2.22), d.h. (5.2) zur Folge hat. 

Der einfache Ansatz (4.2) fiihrt hier bei der Definition (1.2) auf folgendes 


Ergebnis: 
Voraussetzung : 
V1 (CES ee ree) 
mit 
y=) =(D—|C,|)*|G\e, e=(e), o==1, 9 1,2, 9m) NS .) 
und 
SEEMS | wi —Uo| Saas 
ay eae (5.7) 
| ws —Uo| ae ae a sles, e> 
a5, =O gesetzt fiir |w — u|=0 und m5 =1). 
Fehlerabschatzung : 
|u**— vil <e (¢ =1,2,...,m) (5.8) 
und 
| u**— ui| <9; € (Ghee 4 ane (5.9) 


Statt des Vektors » in (5.6) kann man auch die kleineren Vektoren 
y=|(D—C,)4||C\le oder »=|(D— Cae Cale, 


im Falle (D — C2) 720 also z.B. y=(D—C,)+|C,|e verwenden. — Das ge- 
nannte Ergebnis gilt auch, wenn man die Zahlen »; durch die gréBeren (in 4.3 
definierten) jv; ersetzt. 

Benutzt man den Vektor y aus (5.6) und in (5.7) und (5.9) statt »; iiberall 
y = Max »; — wobei »<1 vorauszusetzen ist — so erhalt man entsprechend wie 
beim Gesamtschrittverfahren in (5.9) eine Fehlerabschatzung, welche sich auch 
mit dem Fixpunktsatz fiir kontrahierende Abbildungen beweisen laBt (s. auch [9]). 
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§ 6. Allgemeines tiber die Aufteilung der Matrix G 
6.1. Notwendige Bedingungen. In der Voraussetzung (2.5) ist die Forderung 
: x—Mtx+M-ysSy—Mty+M-x, xy, 
ies) ASB le sit eai= Vi 
_enthalten. Ebenso folgert man aus (3.2) 

Ave |B on Wei Ve 1.210% 

Die gleichen Bedingungen mit (z.T. anderen) Vektoren zo fordern die Voraus- 
setzungen aller in § 2 bzw. §3 genannten Ergebnisse. Man hat also als 


Notwendige Bedingung: Damuit die Voraussetzungen der in § 2 bzw. § 3 
genannten Ergebnisse erfiillbar sind, muB es einen Vektor z mit folgenden Eigen- 


_schaften geben: eo 20, z=|A+Blz, (6.1) 
in § 3: Re0r | Ave | Bi 35 (6.2) 
Besitzt z.B. die Eigenwertaufgabe 
(Acad? weemiAce. Abewin |B pA 


einen Eigenvektor z2o0 zu einem Eigenwert 4 mit O<AS<1, so geniigt dieser 
der Forderung (6.1) bzw. (6.2). 


6.2. Aufteilung fiir Abschatzungen nach §2. In den Abschatzungen des 
§2 treten Vektoren der Form Tu=A™(Bu-7) auf. Will man diese Ergebnisse 
benutzen, sollten also Gleichungssysteme mit der Matrix A einfach zu lésen sein. 
Das ist z.B. méglich, wenn A eine Diagonalmatrix ist wie in § 4 oder eine Dreieck- 
matrix wie in § 5 oder aber von solchen Formen nicht sehr stark abweicht, etwa 
bei Aufteilungen der Art (siehe [15}): 


SchlieBlich kann man auch Gleichungssysteme mit Matrizen der Gestalt? 


c —1 
— CG ==4 
24 c = 1 
— 4 eo =F é c = const 
—1 el 
— Cc 


oft recht einfach lésen, insbesondere dann, wenn m eine Potenz von 2 ist (s. [12]). 


2 Nullen als Elemente sind nicht hingeschrieben. 


40 JOHANN SCHRODER: 


6.3. Aufteilung bei Abschatzungen nach § 3. Bei den Abschatzungen des 
§3 darf A kompliziertere Gestalt haben. Hier bestehen u.U. Méglichkeiten, 
Fehlerschranken fiir wesentlich neue Typen von Gleichungssystemen zu gewinnen, 
fiir welche man bisher noch keine Fehlerschranken kennt. Dies mtiBte noch aus- 
fiihrlicher untersucht werden. sare 

In §3 wird A1=O verlangt. Diese Forderung ist bei der Definition (1.2) 
z.B. erfiillt, wenn 


Po Se m =O fiir alle z 
a es fr =F Pe S| 


=O tur 2=ek k=1 > 0 fiir mindestens ein 7 


gilt und die Matrix A nicht zerfallt [3]. Durch Produktbildung kann man weitere ° 
Matrizen mit A440 erzeugen, z. B.? 


| lt 5 —4 1 
cet 2 —1 —4 6 —4 | 
A= —1 2 —1 = 4 —4 6 —4 1 [= *-(6.3) 
asl 2,-—1 4 —4 6 — 
—4 | 4 —4 a 


Bei der Ordnungsdefinition (1.3), (5.4) gilt A420, wie man aus dem eben 
genannten Ergebnis fiir die Definition (1.2) folgert, z.B. fiir die Matrix? 


lee 
1.2, 4 

2 
Ae Ay FY 
re | 
{ig 


Die notwendige Bedingung (6.2) kann man auch in der Form 
Bie 0 iG mes 


schreiben. Sie ist also um so leichter erfiillbar, je kleiner B~ ist, am giinstigsten 
ware in diesem Sinne eine Aufteilung (2.2) mit B-=O. 

Bei der Ordnungsdefinition (1.2) bedeutet dies also, daB man méglichst viele 
der positiven g,, zu A schlagen soll. Zum Beispiel enthalt die Matrix (6.3) positive 
Elemente auch auBerhalb der Hauptdiagonalen. Dreieckmatrizen A mit A+=>O 
k6nnen ebenfalls positive Elemente g;, mit i=:k aufweisen (s. Beispiel 2). 


§ 7. Zusammenfassung 


Die in den §§ 2 und 3 formulierten Ergebnisse enthalten u.a. Fehlerabschatzun- 
gen folgender Art: 


a) Abschatzung des Fehlers w*— uw, einer Naherung u,= Tuy, ausgehend von 


Schranken fiir die Anderung u,— 1, also Fehlerabschatzungen fiir das Itera- 
tionsverfahren #,,,—= Tu, (m =0,1, 2,...) (Nr. 2A). 


b) EinschlieBung der Lésung zwischen Naherungsfolgen bei Iterationsver- 
fahren mit monoton wachsendem oder fallendem Operator, ausgehend von einer 
Schachtelung der ersten Naherungen (Nrn. 2.2 und 2.3) 
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c) Abschatzung des Fehlers w*— wu, einer Naherungslésung u,, ausgehend von 
Schranken fiir den Defekt d= —Gu,y+r (Nr. 3.2). 


Zur praktischen Anwendung der Ergebnisse auf ein gegebenes Gleichungs- 
system Gu =r hat man zunichst G in G=A — B mit regularer Matrix A aufzu- 
spalten und die Beziehung wv zu definieren. Es ist dann Tu =A1Bu-+ Ay. 
Dabei ist auf folgendes zu achten: 


Die Matrix B soll méglichst ,,klein“ gegeniiber A sein (als notwendige Be- 
dingung ist in allen Voraussetzungen enthalten, daB es einen Vektor z mit (6.1) 
bzw. (6.2) gibt). 

Bei Abschatzungen der Arten a) und b) soll A so einfach gebaut sein, daB 
man ein Gleichungssystem mit der Matrix A leicht lésen kann. 

Bei Abschatzungen der Art b) wird A+ B=O bzw. A1BSO gefordert. (Fiir 
die beim Gesamt- und Einzelschrittverfahren iiblichen Aufteilungen sind in den 
Nrn. 4.2 und 5.2 verschiedene Falle genannt, bei denen eine dieser Ungleichungen 
erfiillt ist.) 

Bei Abschatzungen der Art c) mu8 A1+=ZO gelten, A kann dabei jedoch von 
komplizierterer Bauart sein. Hier sollte B=2O, jedenfalls aber B~ méglichst 
klein sein. 

In den Voraussetzungen der Ergebnisse treten Vektoren auf, welche bestimmte 
Ungleichungen erfiillen sollen. Es seien zwei (allerdings eng zusammenhangende) 
Méglichkeiten genannt, solche Vektoren zu berechnen: 


1. Man benutzt die Relaxationsrechnung, ermittelt die gesuchten Vektoren 
also durch systematisches ,,Probieren“. — Bei b) geht man, um Vektoren %9, Vo 
der gewiinschten Art zu bekommen, etwa von einer Naherungslésung u%» aus, 
bringt Anderungen an und beobachtet deren Wirkung. Bei a) und c) kann man 
z.B. mit v =w =o beginnen. 

2. Man macht fiir die gesuchten Vektoren einen Parameter enthaltenden 
Ansatz und bestimmt die Parameter so, daB die jeweiligen Forderungen erfillt 
sind. Bei b) kann der Ansatz z.B. %» =u )— OV, Yo=%+ew lauten, mit einer 
Naherungslésung u,, festen Vektoren v,@ und Parametern 0, ¢. Bei a) und c) 
setzt man etwa an: v=—06v, w=ew. Zum Beispiel kommt man oft schon mit 
den einfachen Vektoren 0’ =wW = (e’), e=1 ES ONE AL NSS EPPS Ba (pani 
diesen Vektoren folgenden Ergebnisse fiir Aufteilungen wie beim Gesamt- oder 
Einzelschrittverfahren sind in den Nrn. 4.3 und 5.3 formuliert.) 


§ 8. Beispiele 
Beispiele fiir das EinschlieBen der Lésung mit der Relaxationsmethode bei 
monoton wachsendem Operator findet man in [2] und [4] (S. 341ff.). In [10] 
werden Fehlerabschatzungen fiir Iterationsverfahren mit Hilfe von (2.26), (2.25) 


Tabelle 1 
G | A A | B 
pvp i-gau | 
1,82 —2 (0) 9 Die i=) (0) | Oro O 0 
—1 1,84 —1 9 =| Pg aeicraghs hime tall, 0,16 O 
©) —1 1,9 9 =—1 29° Tio ©) 0,1 


42 JOHANN SCHRODER: 


durchgefithrt. Hier ermitteln wir Fehlerschranken fiir ein Beispiel durch Ein- 
schlieBen des Defektes und fiir ein zweites Beispiel mit Hilfe der Formel (5.9). 


Tabelle 2 Tabelle 3 


100, 53 _. | 0,087 95 
86,98 0,0068 | 0,5:/3 | 0,09349 | 0,083 
50,52 — 0,0080 0,5 | 0,08397 | 0,048 


8.1. Beispiel 1. Bei der Randwertaufgabe 
¢@) +1,62(1—#) pO +81—0,. o(—)— 9) —0 


fiihrt das gewohnliche Differenzenverfahren zur Schrittweite h =+ auf das 
Gleichungssystem Gu=r mit G=A—B5, wobei G,v, A und B in Tabelle 1 
aufgefithrt sind. Die Naherungslésung “» in Spalte 1 der Tabelle 2 ergibt den 
Defekt d = —Guy+r in Spalte 2. 

u<v sei durch (1.2) definiert. Dann ist A +2O und B=O. Gilt also (vgl. 


(3.6'), (3.7) 
—GwsdxGw fiir einen Vektor w Zo, (8.1) 


so hat man die Fehlerabschatzung 
—wsu*—uySw. 


Wir machen den Ansatz w=ew mit dem Vektor @ in Spalte 3, fiir welchen 
IU 


Aw =(2—Y3)@ ist [w’=cos + 3). Man erhalt als Gw den Vektor in Spalte 4, 
und die Forderung (8.1) ist erfiillt mit «0,096. Als Fehlervektor ergibt sich 
also w=ew in Spalte 5. 

Die Naherungslosung uw, stimmt mit der exakten Loésung u* in den ange- 
gebenen Stellen tiberein. Das kommt in der Fehlerabschatzung nicht zum Aus- 
druck. Wir kénnen jedoch bei diesem verhaltnismaBig ,,instabilen‘‘ Gleichungs- 
system auch nicht erwarten, da sich vom Defekt ausgehend wesentlich genauere 
Schranken herleiten lassen, denn man kann uw, in den angegebenen Stellen 


andern, ohne da sich die Komponenten des Defektes dem Betrage nach ver- 
groBern. 


Zum Beispiel wiirde. der Vektor w%»=uy—z mit z in Spalte 1 der Tabelle 3 
den Defekt d in Spalte 2 ergeben. Vom Defekt ausgehend kénnte man also 
fiir eine bessere Naherungslésung als wy halten. Fithrt man fiir % eine Fehler- 
abschatzung wie oben durch, so ergibt sich €=0,0693 und der Vektor ® =%@ 


in Spalte 3 der Tabelle 3. Der Fehler von ii, — er ist az — ist tats&chlich von 
der GroBenordnung der Schranke w. 


Man sieht an diesem Beispiel, da8 man aus dem Defekt nicht ohne weiteres 


auf die Giite einer Naherungslésung schlieBen kann, sondern dazu eine Fehler- 
abschatzung braucht. 
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8.2. Beispiel 2. L.CoLiatz berechnet in [4] (S. 367) folgende Randwert- 
aufgabe mit dem Differenzenverfahren: 


AAg=const in dem trapezformigen Gebiet der Abbildung, 


= a =0 auf der langen Trapezseite (« = innere Normale), 
CAG aa) auf den anderen Seiten des Trapezes. 


Bei dem in die Abbildung eingezeichneten Dreieckgitter erhalt man das Glei- 
chungssystem Gu =r mit der Matrix G in Spalte 1 der Tabelle 4 und dem Vektor 7 
in Spalte 2. i A 


2A | 
Fig. 1. Eingespannte Trapezplatte 


uv sei wieder durch (1.2) definiert. Zur Fehlerabschatzung benutzen wir 
das in Nr. 5.3 fiir den einfachen Ansatz (4.2) formulierte Ergebnis. Von der 
Naherungslésung “) in Spalte 3 ausgehend, berechnet man mit dem Einzel- 
schrittverfahren die erste Naherung u, in Spalte 4 und die Anderung u,— Mp» in 
Spalte 5. Der Vektor » aus (5.6) steht in Spalte 6. Wir erhalten 
| wi — uo] 
ey 


é= = 0,001 
und die Fehlerabschatzung 
|u**—_ uwil<n Ceo ae (8.2) 


Der verhaltnismaBig groBe Wert »,; hat keinen EinfluB auf ¢, da |w{—wé| sehr 
klein ist. 
Man kann die Fehlerabschatzung etwas verbessern, indem man den aus 


12-2 OF 0 7 
[58 o)s-(3 
2-12 11 0 


zu berechnenden Vektor » =(D —C,)1|C,| statt » benutzt. Dies ist méglich, da 
(D —C,)41=(f;,) 2O ist. Das einzige positive Element der Matrix D—C, auBer- 


Tabelle 4 
?(S) 
12 —6 1 1 0,241 0,240 583 — 0,000 417 0,584 0,584 
= 3 3 == 3 1 0,394 0,394 094 0,000 094 0,594 0,594 
DB —— GO) {15 1 0,477 0,477 087 0,000 087 0,754 0,542 


halb der Hauptdiagonalen, g,,=2, hat nur Einflu8 auf f,;, und eine kurze Rech- 
nung liefert /,,=0. Mit dem Vektor » in Spalte 7 erhalt man wieder € = ¢ =0,001 
und die Fehlerabschatzung (8.2) mit 9; statt »;. 
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On the Inversion of Integral Transforms 


M. RIBARIG 


Communicated by A. ERDELYI 


1. Introduction 


This article deals with two methods for finding the inversion F(x) of the 


integral transform 
b 


G(y) = J K(y, x) F(x) dx (1.1) 
where G(y) is a regular analytic function, defined in an open domain & of the 
complex plane y. 

These methods are the construction of the function F(x) (7) from the values 
of the derivatives G"(y)) of the function G(y) at a point y»€S, and (ii) from 
the values G(y,;) of the function G(y) at the points y;€ S. Both methods have 
been described in the literature many times, usually in connection with the 
construction of F(x) from a complete set of orthonormal functions [/, 2, 3]. It 
is the aim of this note to show that under certain conditions F(x) can be con- 
structed without any regard to the completeness of the system by which /(x) 
is approximated. 

The first method is based on the following consideration: Let the kernel 
K(y, x) and an inversion F(x) of the integral transform (1.1) be such that, for 
any vo€ S, the order of integration with respect to x and differentiation with 
respect to y can be interchanged. Then, the value of the derivative G"(y 9) is 


b 
equal to the integral [ K(y,, x) F(x) dx. This integral can be interpreted as a 


scalar product of the functions F(x) and q(x) K" (yo, x), where g(x) is a suitably 
selected weight function. Since the analytic function G(y) is uniquely determined 
by the values of its derivatives G(y9), equation (1.1) may be replaced by the 
values of the scalar products of F(x) with the q(x) K"(y 9, x) and by the require- 


b 
ments that the transform { K(y, x) F(x) dx be regular and analytic for ye S 


and that the order of integration and of differentiation be interchangeable. Thus, 
it should be possible, with suitable choice of the weight function q(x), to ap- 
proximate in mean the function F(x) by the set {g(x) K" (yo, x)}. 

For the second method we can say: Since G(y) is a regular analytic function, 
it is determined by the values it takes at an infinite set of points {y,} C S having 
an accumulation point in G. The scalar product of the functions F(x) and 
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q(x) K(y;, x) is equal to G(y;); hence, equation (1.1) may be replaced by the 
scalar products of F(x) with the g(x )K(y,, *) and by the requirement that the 


b 
transform f{ K(y, x) F(x) dx be regular analytic for y€ ©. Thus we expect that 


we can approximate F(x) by a set {q(x) K(y;, x)}- 


2. Construction of F(a) from derivatives of G(y) 


A Hilbert space L%(a, b) is formed by all measurable complex-valued func- 
tions f(x) having finite norms ||f(x) ||, where 


wa 
IIF(~) leaf lite x)? qu" (x) dx (2.1) 


and the weight function g(x) is measurable and positive almost everywhere. 
The scalar product (f, g) of two functions /, gel; (a, b) is defined by 


(f(x), 8 (*)) =S f(a) B(x) a7 (x) dx. (2.2) 


A function F(x) satisfying (1.1) can be constructed from the derivatives 
Gy») of the function G(y) if we can find a weight function g(x) such that the | 
space Li (a, b) satisfies the following conditions: 

41) At least one of the solutions of (1.1), say A(x), belongs to L7(a, b). 

2) For some yy€ S 


q(x) Kyo, x) E Lj (a,0), 7 =0,1,2 (2.3) 
3) For every f(x) € Li (a, 6) 
b 
= | K(y, x) f(x) dx (2.4) 
is a regular analytic function of y in the domain ©. 
4) For every /(x) € iiss b) 
= J KO (Vora) (ala en meee Ota reas (2.5) 


If the conditions 1) and 4) are satisfied, we can express the scalar products 
of F(x) with the functions g(x) K"(y», x) by the derivatives of Gy). 25, 


G (Yo) = (F(x), ¢(*) K (yo, x)). (2.6) 


Now, from the set {q¢(x) K" (yy, x x)‘, an orthonormal system of functions {p, (x) } 
is formed with 


Px (*) = 9 (x) di Ce, KO (yg, x), (2.7) 
where a due to the possibility of linear dependence among the function 


q(x) K (yg, x). The Fourier coefficients a, of the function (x) can be expressed 
by the derivatives of G(y), 


a, = (F(x), o a) Dies, GO (yo). (2.8) 
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It follows from BEssEL’s inequality that the series 


ae es > a, pu (2) (2.9) 


converges in mean to-some function F,,(x) € Li (a,b). By condition 3), 
b 
Gaply) = J K(y, *) Fip(x) dx (2.10) 


isa regular analytic function in the domain &. The Fourier coefficients of F (x) 
and F,,(*) being equal, the derivatives of the functions G(y) and G,,,(y) are equal 
at the point y). Since both functions are regular analytic, we have for all ye S 


G(y) = Gap(y), (2.41) 
and the function F,,(x) is a solution of (1.1). 
The difference of the two solutions (x) and F,,(x) of (1.1), 


Fy (x) = F(a) — Fp (2), (2.12) 


is one of the solutions of the homogeneous equation 
b 
J Ky, *) Fox) dx =0. (2.13) 


Therefore we can conclude: The set {¢(x) K" (v9, x)} is closed in L3 (a,b) if, 
and only if, none of the non-trivial solutions of the homogeneous equation (2.13) 
belong to Li (a, 0). 


3. Construction of f(x) from the values of the function G(y) 
Let the Hilbert space L7 (a, b), (2.1) and (2.2), satisfy the following conditions: 
1) At least one solution, A(x), of (1.1) belongs to L%(a, b). 
2) It is possible to choose an infinite set of distinct points {y;} CS having an 
accumulation point in S so that {¢(x) K (y,, x)} CL (a, d). 
3) For every function f(x) € L(a, b), 


g(y) =f K(y, x) f(x) dx (3.1) 


is a regular analytic function of y in the domain ©. 
From the set {¢(«) K(y;, x)}, the system of orthonormal functions {,(x)} is 
formed with 


Pa(*) = (4) Dieu Ky, a), nek. (3.2) 
The Fourier coefficients a, of the function /,(x) can be expressed as 
a, = (F(x), (2) = DIF, Uy). 3.3) 
It follows from BESSEL’s inequality that the series 


B, (2) = yy, an x (2) 3.4) 
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converges in mean to some function F,,(%) € L? (a,b). By condition 3), 


Gap(y) = J K(y, *) Fap(*) ax (3.5) 


is a regular analytic function in the domain G. The Fourier coefficients of F, (x) 
and F,,(x) being equal, the regular analytic functions G(y) and G,,(y) are equal 
in an infinite set {y,} of distinct points. Hence, for all ye S 


Gap(¥) = G(Y), (3.6) 
and F,,(x) is also a solution of (1.1). 

It is evident from the proof that it is not necessary that the condition 3) 
be satisfied by every {(x)C Li (a,b). Strictly speaking, that it be satisfied at 
least by A(x) and F,,(x) suffices. The same can be said for the third and the 
fourth conditions of the previous paragraph. 

Besides, it follows that the set {q(x) K(y,;,x)} is closed in Li(a,6) if, and 
only if, none of the non-trivial solutions of the homogeneous equation (2.13) 
belong to Li (a, 0). 

4. Examples 


a) The kernel K(y, x) ts a regular analytic function of y € & for every al xb, 
and a uniformly bounded function of x€la,b| and ye S. The difficulty in the 
application of both methods lies in the construction of the weight function g(x) — 
satisfying the conditions of Sections 2 and 3. The problem is simplified if the 
kernel K(y, x) is a regular analytic function of v€ S for x €{a, b] anda uniformly . 


bounded function of x€[a, b] and y€G, and if an inversion /,(x) exists which | 
is integrable, 


NEiAtten (4.4) 


In this case we can always select g(x) so that all conditions are satisfied for the 
application of both methods. We choose an almost everywhere positive meas- | 
urable weight function so that 


A (~)|at (4) <M (4.2) 
and 
b 


J a(x) dx< My, (4.3) 


a 


where M, and M, are arbitrary. Let us verify that such a weight function satisfies 
all conditions for the inversion of the integral transform by both methods. 
Clearly F(x) € L?(a, b) since 


A) [P< My< 09. (4.4). 


The second condition of both methods is satisfied since, for every yC GS, 
lla (*) KO (y, «|| << sup |KM(y, x) | VE Se gi = Oa aaa (4.5). 
asxsb 


for the Ky, x) are also bounded functions of x €[a,b] and y€G, as one can. 
see expressing K")(y, x) by K(y, x) in the form of the Cauchy integral. 
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The third condition of both methods requires the function g(y) to be regular 
and analytic in the domain © for every /(«) € L3(a, b). We shall show that ¢(y) 
can be expanded in a power series arround every point yy€ S. Now, 


Ky, x) = yy 290" KH (yg, x) + 74(Y, o> 2), (4.6) 


0 


and it follows from the expression of 7,,(y, v9, *) as a contour integral involving 
Ky, x) that, for each y)€ ©, there exists an N that 


l7n(¥s Yoo X)| < £0 
for every 


m>N, |y—¥|Se% and x€[a,d], 


where 9<1 and ®% is the distance between yy and the boundary G. Inserting 
the series (4.6) into the integral transform (2.4) we get 


g(y) =i 2%! (F(x), g(x) K (99, 2)) + Ry lv, M0), 
where 0 (4.7) 
R,,(y, Vo) == (f(x), q(x) TVs Vo. x)) , 
and 


|R(Y, Yo)| < VM, ||7(*)]| Es (Y, Yo, *)|<e (4.8) 


for every n>WN and |y—yo|Se%). Thus g(y) is indeed a regular analytic func- 
tion of y€ & for every f(x) € L? (a, b). 
It follows from cquation (4.7), with ~=oo and R,,=0, that the fourth 
condition of Section 2 is also satisfied. 

In the following two cases, the method will be applied to two kernels which 
are not analytic functions of y€ © for every x € [a,b], but are uniformly bounded 
function of x€ [a,b] and yEG. 


b) Inversion of Laplace transforms. Both methods can be applied to the in- 
version of Laplace transforms 


gly) = J exp(— vx) f(a) ax (4.9) 
provided y is in the domain 
Soke ve ke. yea 0, hoM., 
where M and & are arbitrary, and the weight function g(x) is chosen so that 


f(x) €L7(0, 2) 
and "i (4.10) 
gix)exp(— #4) CL, (0,7). 
‘Consequently, the first two conditions of both methods are satisfied. By Schwarz’ 
inequality we get from (4.10) 


fexp(— kx lf @)ldz<o. (4.11) 


Thus, all other conditions are fulfilled too, as follows from the theory of Laplace 


transforms [4]. 
Arch. Rational Mech. Anal., Vol. 3 4 


50 M. Ripari¢: Inversion of Integral Transforms 


For an arbitrary solution of the integral equation (4.9) it is not always possible 
to choose a weight function g(x) which satisfies (4.10), since solutions are known 
which do not satisfy the inequality (4.11) at any value of & [4]. 

If we put g(x) =1 and choose a positive yo, the method of Section 2 becomes 
identical with Tricomi’s method for the inversion of Laplace transforms by 
means of Laguerre polynomials [/]. 

With g(x) =1, the method of Section 3 becomes a particular case of ‘the 
method of A. ERDELYI [3]. 


c) Inversion of Fourier transforms. The inversion of the Fourier transform 


co 


g(y) = J exp(¢y x) f(x) dx (4.12) 


—oco 


is possible if we can choose the domain 
Sg; |Imy|<k, |Rey|<M k,M>0 
and the weight function ¢(x) so that 


f(x) € Li(— 0, 0) 
and (4.13) 
q(x) exp (k|x|) € L3(— oo, oo). 
That such a weight function satisfies all conditions of Section 2 and Section 3 
can be shown in a similar way as in the case of the Laplace transforms. 

Some types of integral equations cannot be solved by the described methods, 
e.g. integral equations of the second kind. However, the methods may be ap- 
plicable to an equivalent integral equation obtained from the original one by a 
proper integral transform. 
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Vaguely Normal Operators on a Banach Space 


G. L. KRABBE 


Communicated by W. MAGNUS 


1. Introduction 


This is the first of two articles concerned with extending to Banach space 
operators certain properties of normal (Hilbert space) operators. 


If B is the Haar measure of a locally compact abelian group G, then the relation 
(Cala i= ae )-y(A)-B(da) for all (x, y) in L?(G) xL*(G) 


defines (when # = 2) the inner product in the space L?(G). If #(L?(G)) is the set 
of all operators T on L?(G) UL? (G) such that 


ey ea lel 9) eetoralle(x oy) in’ e(G) Ee (G)., 


then (L?(G)) will denote the set of all operators of the form 7, +77, where 
T, and 7, are both in #(L*(G)). An operator will be called “vaguely normal”’ 
if it is a bounded ae on L?(G) that has an extension in W(L?(G)). I here 
take G ={0, +1, +2, ..} and consider a certain algebra Ug, of convolution 
operators on J, =L?(G). 2 aie out that %~, consists of vaguely normal opera- 
tors. Generalizing previous results [5,6], I show that Ws, is included in the 
weak closure of the algebra generated by the Hilbert transformation *—H * = 
{{H x],},c¢g, where 


[Hx], 7 » ip ee ( ie n) : (1) 


To each T in Ys, corresponds on operator-valued measure E T and a function / 
such that the relation 


(P= reice = S10 ). B41 (ad) (2) 


holds in a certain sense (see 5.5.). When p=2, this implies that 7 —/(H) (using 
the functional notation customary in Hilbert space). If f 22 then the resolvent 
set of T is void (Theorem 5.4); this extends another property of normal operators 
(15, p. 305]. M. Riesz and E. C. Tircumarsu [/4] have studied certain out- 
standing members of the algebra gq, (see the examples in 5.7). 

4* 


52 G. L. KRABBE: 


2. Three measure-theoretic lemmas. 

Since this paragraph deals with immediate consequences of well known results 
in [1], the explanations will be very sketchy. It will be convenient to say that 
(E, S,v) isa (TF) complex measure space if y= do, %, (where n=1,...,4 and 
w, =exp(inn/2)) and if (EZ, S,»,) isa totally finite measure space [J, p. 73 and 
p. 31] for all »=1,...,4. Let w be Lebesgue measure on xX =[0,4),-whilese 
denotes the o-ring of all -measurable subsets of X. All functions are complex- 
valued. 

2.1 Lemma. Suppose that A€L1(u, X) and let » be the function E +» (E) 
fined on X b 
defined on y Repeat Sih 

E 
then (X, X,v) is a (TF) complex measure space, v<m, and dv/du =f. 

Note. By writing h=)>\w,h,, (n=1,..-.,4) with h,20 and h,€L*(u, £), 

we can define »,(E) = h,, du; accordingly y= >), »,. From [1, p. 127 (1)] fol- 
E 

lows that »,<j for all m (this is our definition of <j). From [J, p. 133] follows 

that dv,/du=h,, whence h=))o, dv,/du=dr/du; the last equality embodies 

our definition of dy/dyu. 


2.2 Lemma. Suppose that (X, X,v) is a (JF) complex measure space and 
v<. If gis a bounded Borel function in L1(y, X), then 


fedvafe- du. (3) 


2:8 Note. “Here 14 (9%) = 0{L) 0,2) Ay ..4 4) Letegs andig, beathe 
real and imaginary parts of g; the g, are then bounded Borel functions and 
g,CcLi(v, X). By [1, p. 134, Theorem B], the equality (3) holds for g=g,, v=», . 
A suitable linear combination of these results yields the conclusion, if we define 


1 4 
fg -dy me XP On J 8 : AV. 


2.4 Notation. If g and A are functions, then (go A) denotes the function 
%—>g(A(x)) (called gd by Harmos [J, p. 162]). The set L®(w, X) of all «-measur- 
able, essentially bounded functions will also be written L®(X). 


2.5 Lemma. Suppose that d€L™(u, X), and let w(A) denote the closure of 
the range of A. If ¥) denotes the class of all Borel sets of @(A), then the set 
C= 4k) -F EY} is a o-subring of X [7, p. 162]. If gis a bounded Borel func- 
tion defined on w(A) then (go A) EL1(y, S) and 


J (go A) -dy=Jg-d(voA). 

a w (A) 
Note that A is a measurable transformation (in the sense of [7, p. 162]). With 
the notation of 2.3, it follows from [/, p. 162, Theorem B] that all (g,0A) are 
S-measurable ; but, since any bounded G-measurable function belongs to L1(6, S) 
whenever (X, S, 8) is a totally finite measure space, it results that (g,0A) EL*(v,, S). 
Therefore (go A)€L4(v, S). The rest of the proof follows from [J, p. 163] by means 
of the indications given in 2.3. 
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3. Laurent operators 


Suppose that peEL*(X), X=[0, 1], and let Y m denote the sequence ‘ile Dnt CG 
of Fourier coefficients of ; 


[Ye], ra enn (8) u(a9) =~ (mG). (4) 


i {x ),c6; I shall denote by y® x the convolution [Y g]* x; in other words 
Y exr= {lp ® CAM €G> where 


[p@x],= 2% LY Puy: (5) 


veG 
First studied by Torpiitz in 1907, the operators pa={b Cl,> ob} were called — 
“Laurent operators’”’ by F. Riesz [12, pp. 171—176, 152—155]; they are repre- 
sented by Toeplitz matrices (HARTMAN & WINTNER [4]). It is assumed through- 
out that 1<p<oo. I shall write p’=#/(p—1) and* 


(x|¥) = 20 Tn ((x, y) El, X ly). (6) 


If (b,c) €l,xl,, there exists two functions B and C in L?(X) such that b= 
[YB] and c=[YC]. Therefore [Yp]*b=[Yq]*[YB]=[Y(p-B)] and ac- 
cordingly 


(p ® bc) eC ROTC eA Sts OTs (7) 


The first equality comes from the PLANCHEREL theorem [7, 36D and 38B], the 
purpose of the second equality is to correlate with the notation (due to HALMos) 
which was used in § 2. It is easily seen that mg_ is a bounded operator on J,. 

3.1 Definition. Let y(E) be the characteristic function of a set E; then M4 
will denote the mapping o->M4 (oc) =7(A4 (c))a_ of the Borel sets of the plane 
into the space of bounded operators on 1,. On the other hand, M(E)=y(E)e2 for 
all E in & (see § 2). 

3.2 Theorem. Suppose that ACL™(X). If g is a bounded Borel function 
defined on @(A) (see 2.5), then the relation 


((g0A) © 6]<) =F g(A(B)) (L(A) Ble) = fg) (MM (42) |e) 


holds when (b, c) Ely X ly. 
Proof. Fix b and c, and write »(E) =(M(E) b|c). A successive application 
of 3.1 and (7) yields 
v(E) = (x(E) @b|c) =fz(E)-B-Cdu=J(B-C)- du. 


E 
Since (B-C)€L'(X), we can infer from 2.1 that (X, X,v) is a (JF) complex 
measure space, »<ju and 


“* —B-C  (u-almost everywhere). (8) 
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But, by 2.5, f=(goA) is a bounded function in L*(», 6), and 


1 


frawafi dua [p BT du=(f00l9. (9) 


The first equality comes from 2.2 (3), the second equality from (8), and the third 
comes from (7). On the other hand, 


ff-dy =f (0) (M (dd) b| c) =f g(a)v(A4@A)). (10) 


The first equality results from the definition of y, while the second comes from 
2.5. Since » (A4(a)) = (x (A 1(0)) bc) = (4% (0) 6] ¢), a comparison of (9) and (10) 
concludes the proof. 

3.3 Remarks. If A €L®(X) it is easily checked that M4 is a compact complex 
spectral measure in the sense described in Hatmos’ book on Hilbert space [2). 
Incidentally, the conventions used in [2] permit the more concise notation 


1 
Ag,=fA(O)-M(a9) =f a-M4 (dd) 
0 
in the case g(A) =A. 
3.4 Corollary. If AC L™(X), b€l, and m€G, then 


1 
LAO by =A aN ard Nan. 
0 
Proof. Substitute c= {4,,,},<¢¢ (Kronecker delta) in 3.2 with g(A) =A and set 
N,, (E) = (M(E) 6|c) =[M(E) 5. 


4. The property S, 

Let ©, be the class of all bounded operators which map /, into itself. If Q 
is an operator defined on /,Ul,, then Q, will denote its restriction to /, (ie., 
Q,x =Qx for all x in /,); the class €,~ &,, consists of all such operators Q which 
satisfy 0,€ ©, and Q,€ Gy. 

If AcL™(X), define A’={0— A(1 —8)} and consider the operations A > A’ 
and A->A. From now on, mw-equivalent functions are identified. I shall denote 
by S, the family of all subalgebras § of L®(X) which are closed under these two 


operations and which are such that A—Ag is an algebraic isomorphism of § 
into a subalgebra ¥~ of E,Gy. 


4.1 Remark. It is easily checked that L®(X)€S,. 
4.2 Notation. If Q=Ag I write Q'=(A’)g and Q=(A)g. If (x,y) Cl, Xlyy 
then <x, vy) =(x|¥) (see (6) and recall that ¥= {Vz hn co): 


; D seth Suppose that §€S, and QE Bq. If (x, y) El, xl, then (Ox, y> = 
BOs ys 


Proof. Suppose that Q = Ag and let Jy denote the set of all sequences x= 


Y 
woe EG 
such that x,—0 whenever | 


is sufficiently large. If a and 0 are in J), then 


22 bn a [VAn 4 =Z 34, bn [VA], —n- (11) 
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This follows from the relation [YA],,=[YA']_,, (which in turn comes directly 
from the definition (4)); in view of 4.2, (6), (5) and (11) we therefore have 


<Qa, b> =<A ®a,b> =a, A’ @ b> = (a, 0'bY (12) 


when (a, b) El, XJ. Since J, is dense in any space /,, there exist sequences {a}, 
and {b)}, with (a, b)€1,xJ, such that 


im IIx a™||,=0, lim [ly = 2 || —=0. (13) 


But Q,€ G and (Q’),€ &,, (since Q’ = Bg and B=A’€ &€S,). From (43) there- 
fore follows that 


<%,(Q')y y> = lim lim <a, (Q’) » B® (n >co, k+ov). 
Since Q,x =Qx and (Q’),,y=Q’y, a reference to (12) completes the proof. 
4.4 Theorem. Suppose that § CS, and Q€ Be. If (x,y) El, xl, then (Qx|b) = 
(x|Q0). 
Proof. Note first that 


LO 05 = 2 AB. Ole aes (ee Les 


where for typographical reasons [A denotes /. But from the definition (4) 
results easily that E( RYLA 1) =[YAl],,, and it can therefore be concluded from 
4.2 and 4.3 that 


(Qx|b) =<Qx, b> = <x, Q'dy =X4,1\2 b,[YA],—n} = (x|Qb). 


5. The main results 


Denote by (#R) the set of all real-valued functions in §. If $¢€S, and 
Q€(SR)e, then Q=@Q and from 4.4 follows that 


(Qx|b) =(x|Q) for all (x,b) in 1,Xxly. (14) 


Thus (SR) qc A(l,) (see § 1). If QE Sq, then it is easily seen that Q=Bg +1Cg, 
where both B and C are in (§R)g. Thus FgcWV(i,) (see § 1). 

Accordingly, if T is in the set Jg,={4@,:AE B}, then T=Ag,y for some 
A in %, and T has the extension Q=Ag€%q which belongs to W(l,). We 
thus obtain 

5.1 Theorem. Jf §ES,, then Sep consists of vaguely normal operators. 

5.2 Remarks. lf § € Sg, then the members of (4X) @, are selfadjoint operators, 
and %» consists of normal operators. By 4.1, this holds true in the particular 
casepey = Ly (Xx): 

5.3 Notation. Let & be the set of all functions of finite variation on [0,1], 
and let X be the set of all functions which are equal almost-everywhere to some 
member of %. It has been shown in [6, 6] that &¢S,; accordingly, all operators 
in Ug, are vaguely normal. 
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If TE &,, the notion of adjoint operator 7* of 7 is a classical one [15, p. 169]. 
The usual identifications yield 


CPG renee when (x, y)€l,xly. 
If TEAg,, then T=Ag, for some A in XY, and it follows from 4.3 that 
T*=(Aey. (15) 


Let I7(A) denote the set of all complex 4 such that mu ({#:4 (9) =/}) 20. TAR 
MAN [3] has pointed out that, when 1<P~<2, then 


Po(Bg,) HIB) Po(Bo yp) (if BE A) (16) 

where Po(T) denotes the point-spectrum of 7. 
5.4 Theorem. If 2<q<oo and TEUg,, then the resolvent set Rs(T) is void. 
Proof. By hypothesis there exists a member A of %& such that T=Ag,. 


Observe that Rs(T) ¢Po(A,) CHA’) = IA) <Po(Ag,): 


the first inclusion follows from (15) and the fact that Rs(T) <Pa(T*) (see [/J, 
p- 304, Theorem 5]), the second and third inclusion is obtained from (16) by 
setting p=q’, while the equality is easily established*. We have thus obtained 
Rs(T)<Po(T), whence the conclusion Rs(T) =0 (since Rs and Po are always 
disjoint). 

5.5. From 5.2 follows that Ag, is a normal operator (in the usual sense) 
whenever A €Y, and all the operators encountered in § 3 are normal operators. 
Suppose that TE%g,; T=Ag,, ACU. Let E T be the resolution of the identity 
for Ags. From 3.3 and [2, p. 71, Theorem 1] therefore follows that 


fa-M4 (da) =fa-ET (aA). 


Since both M4 and E*? are compact complex spectral measures, it follows that 
M4=E? (see [2, p. 65]). If A is now the function J(9)=8 and J=Ig,, then 
the same reasoning (again using 3.3) yields M=E/. A small computation will 
verify that — 2a7J is the operator H defined by (1); consequently E? = —2aiM. 

Now (8,c)€l,0(1,x1,) implies (A @b|c)=(Tb|c) and 3.2 now yields the 
relation 

1 
a: T \ ae (ee) 
(Td|c) =fa-(E (dA) bc) =| fe) (E# (49) | c) 
0 

for all pairs (0, c) in J, (1, xl,-). The equation (2) in the introduction is intended 
as a concise rendering of the above. 


5.6 Corollary. Suppose that T€ Ag, and let o(T) denote the spectrum of T. 
HarTMAN [3] has shown that there exists a function B in & such that T=Bg, 
and o(Ba,.) =o0(T) =a(B) (see 2.5). As in 5.5, let E? be the resolution of the 
identity for Bg,. Accordingly, if g is a bounded Borel function defined on o(T) 
with (go B) EU, then (go B)g, is an operator V which (in view of 3.2 and Sey) 


* Note that {8:A’(9) =/A}={1—a:A(a)=/} and “M(E) +0 if and only if 
M(t1 —aiae Eh) +0. 
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satisfies the relation 


(Vb|c) rie -(E* (dA) b|c) (17) 
for all pairs (b,c) in /,>(/,xJ,). In particular, if 1<p<2 and c={6,; hi c@ (see 
3.4), then 

ee ee (18) 


holds for all 6 in 1, and all » in G. 


5.7 Examples. Consider the translation operator x-+1x defined on /,, where 
TX={Xn41}neG- It is easily checked that tr=Ag,, where A(#) =exp(— 22198). 
It was proved in [6] that, if A(#) =>) /,,0” is an entire function, then 4g,= 
Di An(e >)”. Accordingly, we may write 

co 
T= 2a (— 2211 g,)"=expH 
since H=—2aiIg,. Let loga be the branch of the logarithm such that 
log exp(tA)=7A for —2ma<ASO; then I=(t/2m)(logoA), whence Ig,= 
(2/22) (logo A)@,, and from (18) it follows that the curious relation 


[Hx], => 4 = Plog a-[E() bn 


ISG |a|=1 
von 


holds for all 2 in G and all x in J, (when 1<f£S2). 
If we take A™ (9) =exp(—22ia0+271«), then (when |«| is not an integer) ; 
[A @ x], = ee Si ty ———— (n€G). 
p=—oOO 
Both operators H and (A), have been studied in a famous article by M. Riesz 
[17]. If «=1, then A™=-—A, so that 


6. The weak closure theorem 

In case # =2 it is customary to denote by g(7) the operator V defined by 
(17) for all pairs (b,c) in J,x/,; therefore (go A)g.—g(A@_). In particular, to 
any T in %g~_ corresponds a bounded Borel function g such that T =g(Ig.) 
(this is obtained by taking A =J and T=gg,). In the general case (6=2), 
the result t =exp(— 271Ig,) of 5.7 illustrates the fact that, if & is the class of 
absolutely continuous functions, then §@, is in the uniform closure of the algebra 
generated by Ig, (this was proved in [5]). A weaker result relating to the whole 
algebra 2, will now be established. 


6.1. The algebra U generated by Jg, is the (complex) linear span of the set 
{Ue 5)": #=0, ies it A member T of ©, is in the weak closure of U if there 
exists a sequence {7}, in U such that (Tx|y)=lim(Tx|y), too when 
(x, vy) €L,X ly. 

6.2. Two more remarks will be in order. From the definition (5) follows that 
Ag,=Bg, if A =B w-almost-everywhere; accordingly, Ug, = Be» (see 5.3). It 
is clear from STECKIN’S paper [13] that the isomorphism A—Ag, of & into ©, 
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ss continuous in the following sense: there exists a number / such that 
I|Acollp S Aol All. (19) 
where ||A||, denotes the total variation of A. 

6.3 Theorem. If TEUg,, then T is im the weak closure Ul of the algebra 
generated by Igy. 

Proof. lf TEX», then TE Veo» (see 6.2) and we can write T=fop top 
with fg, and ge,» in (VR)a» (see §5). It is easily seen that any complex linear 
combination of members of Ul belongs to U; since T is a combination of mem- 
bers of (VR) gp, it follows that, in order to prove T€U, it will suffice to establish 
(VR)g»_CU. To that effect, write T= (Kap, where P, is a real-valued function 
in &. Since P, is continuous j-almost-everywhere, we have 

FP, = lim wy, (u-almost-everywhere) (20) 


where P is the ec Bernstein polynomial of P, [8, p. 4]. But, from Aye (VR) 


Uj 


follows [8, p. 25 or 9, Satz 6] that 
sup {||Pl|p: «== 0} =My <0. (21) 


Take any 0 in the set J, of all sequences 0 such that b,,=0 when || is sufficiently 
large. Now |P|<||P||,S% for s-+0; take any m in G, then the measure N,,, 
defined in 3.4 enables us to write 


{P,-aN,, =lim [ P-dN, 


from LEBESGUE’S dominated-convergence theorem (recall that N,,<wy, so that 
(20) holds N,,-almost-everywhere). From 3.4 therefore: 

[P, 0], = lim [P@)b], (bEl,, mEG). (22) 
The second part of the proof begins by taking an arbitrary x in /, and observing 
a IRI) Shelby to-M EO) (23) 


(from (19) and (21)). Fix m for the time being and consider the linear functionals 
~.={«El,>[P ® x],,}; from (23) follows that sup {||q,||:44+=0}<k)My. Since 
moreover @%(b) =lim y,() (from (22)) for all 6 in the set J) which is dense in ly, 
it results that @»(x) =lim y,(x) (too) when x €l, [10, p. 131]. In other words; 
for any x in /,, the property [F, ® x],,=lim[P ® x],, holds for all m in G. In 
view of (23) and [J0, p. 139, Satz 6], this in turn implies that 


(Py ® x|y) =lim (P ® x|y) when (x,y) €l,xly. 

t—> oo 
Now FP, is a polynomial >) /, 1" (6+=0) and therefore P®x—=Tx, where T= 
(Aap Ade »)” (since the mapping A >A @, is a homomorphism, as pointed 
out in 5.3). But then the 7) are in the linear span U, and since Tx =P, @ x: 

(I «|y) =lim (T, «|y) when (x, y) €l,x ly, 

t—> Oo 

which by 6.1 concludes the proof. 
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Sulle Funzioni di pin Variabili a Variaxione Limitata 


LuciANo DE VITO 


Memoria presentata da G. FICHERA 


Il presente lavoro é originato dalle lezioni tenute negli anni accademici 
1956—57 e 1957—58 dal Prof. FicHERA presso I’Istituto Nazionale di Alta 
Matematica. Esso ha per oggetto la dimostrazione della equivalenza di diverse 
definizioni di funzione di pit variabili a variazione limitata. 

Prima di enunciare il tipo di problemi considerati in questo lavoro € opportuno 
introdurre alcune convenzioni che ci permetteranno di rendere pit spedito il 
discorso. 

Indicato con S” lo spazio euclideo ad 7 dimensioni, denoteremo con P il 
punto generico di S” e con x,,..., x, le coordinate di P rispetto ad un fissato 
sistema cartesiano ortogonale. Sia h un intero tale che 1/7; indichiamo 
con «” una combinazione di classe # dei numeri 1, 2, ..., 7i cui elementi a, ..., a 
pensiamo disposti in ordine crescente: a4,<---<a,. Se non é h=rs, indicheremo 
con },,..., 0,-,1'(y —A)-upla dei numeri interi compresi tra 1 ed 7, non contenuti 
nella combinazione «) e tali che b,<---<b,_,. Indicheremo inoltre con é e 
rispettivamente una h-upla ed una (y —h)-upla di numeri reali e precisamente: 


(1a) é = (%,,» brat KXq,)> 1) ss (%>,» eas Xb,-n) i 


Supposto 1<h<y, se u é una funzione definita in un insieme di S”, conver- 
remo di indicare u(P) anche con il simbolo w(&, 7). 

Sia ¢ una variabile, prescelta fra le x,,, ..., ,, che costituiscono la j-upla &. 
Se € h>1, indicheremo con y 1’(h—1)-upla costituita dalle h—1 rimanenti 
variabili (sempre disposte in ordine di indice crescente). Nel seguito, dovendo 
esplicitamente indicare la dipendenza della funzione w dalla variabile ¢ ado- 
pereremo, per denotare la w(P), anche il simbolo w(t, y, 7). Se h=1, si converra 
porre: u(P) =u(t, yn) =u(&, ); se h=r invece: u(P) =u(t, y) =u(&). Inoltre, se 
I é un intervallo di S1, indicheremo con J , la proiezione di J sull’asse S®™ della 
variabile ¢, cioé l’intervallo descritto da ¢ in S!), quando P varia in J. Analoga- 
mente indicheremo con J; la proiezione di J sulla varieta lineare S$” di equazioni 


Xp, ns | Om aae ys con I, ie ee la proiezione di J sulla varieta lineare 
n) . . 

Si ahs equazioni %,,=0,..., X%,,=0 ed infine, se h>1, con I, la proiezione 

di I sulla varieta lineare ‘ste =i) di equazioni /=0, %,,=0,...,4%,,=0. Se 


1 
Usando il termine intervallo intenderemo sempre riferirci, salvo esplicito avviso 
in contrario, ad intervalli superiormente aperti; pertanto, con la locuzione intervallo 


di S® intendiamo un insieme dello spazio euclideo ad vy dimensioni definito da 
condizioni del tipo c¢;< 4;<d; con d;>c; (i=1,..., 7). 
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1<h<y, i punti di S™ verranno indicati, secondo che pit torni comodo, con 
uno dei due simboli, (¢, y, 7) e (€, 7), i punti di S{” con € oppure con (¢, y), i punti 
dell’iperpiano Seat. d’equazione ¢=0, con (y, 7). 

Se h=1, per i punti di S” verranno usati i simboli (¢, y) oppure (€, 7), per i 
punti di Sf? i simboli to & Se h=z, per i punti di S” adotteremo i simboli & 
oppure (t, y). In tutti i casi, per i punti di S® useremo la lettera ¢. 

Se 1<h<y, indicato con J un intervallo di S”, con I, xJ, denoteremo 
Vintervallo di SY;", descritto da (y,) al variare di y in ford yy in: T. 

Quando, nel seguito, in corrispondenza ad una funzione u(P)=u(x,,..., %), 
considereremo una coppia di interi # e k con 1<hfS<r, 1<kS<y, converremo 
sempre di fare le posizioni (1,1) e di assumere t=~,. Indicheremo inoltre con 
Pea NUplan (its: apa, Lean sak): 

Cio posto, una funzione u(P), definita nell’involucro dell’intervallo J di 
S 2 ed ivi continua, dicesi a variazione limitata secondo TONELLI rispetto ad X,, 
oppure, come sempre, diremo nel seguito, [CVL],a) , se: 


1°) fissato il punto 7 quasi ovunque in J,,? la funzione di ¢: u(t, 4) é VL in J,; 

2°) la funzione V(y) che coincide con la variazione totale dell’anzidetta 
funzione nei punti 7 di J, in corrispondenza ai quali u(t, 7) ¢ VL e che assume 
il valore zero nei restanti punti di J,, € sommabile nell’intervallo J,,'. 


Il Prof. FicHERA ha dato la seguente definizione di funzione a variazione 
limitata®: 

Una funzione u(P), definita nell’involucro dell’intervallo I di S” ed ivi 
continua, dicesi a variazione limitata rispetto a x,=t se, fatte le posizioni (1,4) 
e posto, in corrispondenza ad ogni intervallo-dimensionale J <I: F(J) =fu(&)y,do 

Jil 


(ove Y J indica la frontiera di J, v, la componente secondo l’asse ¢ del versore 
normale esterno ad YJ, dw é l’elemento di misura ipersuperficiale su Y J), la 
funzione d’intervallo F(J) é VL in I®. 


2 Per involucro di un insieme J s’intende l’insieme dei punti che hanno distanza 
nulla da J. L’involucro di J sara indicato con I. 

3 Le locuzioni ,,quasi ovunque“ ,,funzione misurabile“ ,,funzioni equivalenti sono 
qui sempre da intendersi nel senso di Lebesgue. 

4La funzione V(7) € misurabile in J,; cfr. [6] p. 430. 

5 Cir. [6] p. 428. 

6 Sia F(J) una funzione definita nella famiglia {J} degli intervalli contenuti in J. 
La funzione F(J) dicesi a vaviazione limitata in I se, considerata un’arbitraria de- 


n? 


composizione 6 di J in intervalli di S™: IT=f,+-:-+J e posto o(6)=))|F(fj) 
j=1 


la variabile o(6) riesce superiormente limitata al variare comunque di 0. L’estremo 
superiore di o(d) prende il nome di variazione totale e sara indicata con Vp(J). E 
evidente che, se F é VL in J, essa é tale anche in ogni J/¢J. Resta cosi definita in 
{J} la funzione Vp(J), variazione totale di F(/) in J, la quale € non negativa ed 
additiva se F é additiva. Una funzione F(j), additiva e VL, é derivabile in quasi 
tutti i punti di J e la sua derivata F’ é misurabile e sommabile in J. Nel seguito, 
usando il simbolo F’ o parlando di derivata di una funzione d’intervallo F, intenderemo 
sempre riferirci alla derivata calcolata nel modo seguente: se P € un punto di J, 
indicato con I» , lintersezione di J con l’intervallo quadrato avente come estremo 
inferiore P ed il lato di lunghezza /, porremo: F’(P) = im nF (Lp, 7) se questo limite 
esiste finito e F’(P) =0 nell’altro caso. Analoga convenzione faremo per le derivate 
delle funzioni di una variabile. 
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Diremo che una funzione wu, verificante tale condizione, ¢ [CVL],,, m J. 


Il Prof. FicHeRA ha dimostrato, nel caso 7=2, lequivalenza fra le due defi- 
nizioni di funzione a variazione limitata ora ricordate’. Questa dimostrazione € stata 
estesa al caso dir qualsiasi in un mio recente lavoro’. I] Prof. FICHERA, nel predetti 
corsi di Analisi tenuti all’Istituto Nazionale di Alta Matematica, ha generalizzato 
le definizioni [CVL]qw,,, [CVL Jam, z considerando funzioni soltanto sommabili 
anziché funzioni continue e precisamente ha dato le seguenti definizion1: 

A. Una funzione u(P), definita nell’involucro dell’intervallo J di S ” la quale 
sia ivi misurabile e sommabile, dicesi [VL],,, se, tra le funzioni ad essa equi- 
valenti® ve n’é almeno una — che, per semplicita di notazione, seguiteremo ad 
indicare con u — tale che risultino soddisfatte le seguenti condizioni: 

a) per ogni 7 €/, la funzione di ¢: w(t, 4) € VL in I,; 

b) la funzione V(y), che rappresenta la variazione totale dell’anzidetta fun- 
zione, € maggiorata da una funzione misurabile e sommabile in J,1°. 

Si dimostra (cfr. teor. XII) che le derivate parziali rispetto a ¢ delle varie 
funzioni equivalenti ad u e soddisfacenti le a) e b) sono misurabili ed equivalenti 
tra loro in J. Allora si assume, per definizione, come derivata di u rispetto a t 


la derivata parziale rispetto a ¢ di una qualsiasi delle funzioni equivalenti ad u 
Ou 
Ox, 


e verificanti le a) e b). Questa derivata sara indicata con il simbolo: 


a), k& 

B. Una funzione u(P), definita nell’involucro dell’intervallo I di S”, la quale 
sia ivi misurabile e sommabile, é [VL], ;, se, tra le funzioni ad essa equivalenti, 
ve n’é almeno una — che per semplicita seguiteremo ad indicare con u — tale 
che, posto per ogni intervallo y-dimensionale JCI: ®y»(J) =f u(»,dw, la 
funzione ®,»(J) risulta VL in TI. ft 

Si dimostra (cfr. teor. XIII) che, se w ed w sono due funzioni equivalenti e 
verificanti la detta condizione, le corrispondenti derivate B; e P. sono equi- 
valenti in J. Allora si assume, per definizione come derivata di u rispetto a t, 
la derivata ®,, relativa ad una qualsiasi delle funzioni equivalenti ad w e veri- 
ficanti la condizione in discorso. Tale derivata verra indicata con il simbolo 
ou 12 
hee 

Il Prof. FicHERA ha poi anche proposto la seguente generale definizione di 
funzione a variazione limitata!’, che include quella di funzione [VL]am,, € quella 
di funzione [VL], ,: 


Ci ON 50, HSA, 

8 Cfr. [5] p. 100. 

® Due funzioni si diranno equivalenti se sono quasi ovunque eguali. 

Questa definizione @ equivalente a quella di funzione a variazione limitata 
generalizzata secondo CEsarI-ToneELLI; cfr. [4] p. 299, [18] 18k, Selo 

‘1 Si tenga presente che, per il noto teorema di Fubini, tra le funzioni equivalenti 
ad una funzione misurabile e sommabile in 7 ve n’é almeno una che risulta misurabile 
e sommabile su tutte le sezioni coordinate di J cioé sulle intersezioni di J con iperpiani 
parallel ad iperpiani coordinati. 

Si dimostra facilmente che ogni funzione w la quale sia [CVL]q,% (che sia 
[CVL]qin, x) € anche [VL ]q,% (€ anche [VL jam,%) €, viceversa, ogni funzione continua 
re cran (che sia [VL], %) € anche [C VL Jae, a (€ anche [(CVL]q, ~); cfr. [7]. 

Yad 
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fissati gli interi h e k, con1<AS<7,1<kS<re la combinazione «™” contenente 
Vintero k, una funzione u(P), definita nell’involucro dell’intervallo I di S”, la 
quale sia ivi misurabile e sommabile, dicesi [VL], ,, se tra le funzioni ad essa 
equivalenti ve n’é almeno una — che, per semplicita di notazione, seguiteremo 
ad indicare con u — tale che siano soddisfatte le condizioni: 


1,0) posto, per ogni 7€J, e per ogni intervallo J, contenuto in J,: 


Dm (Jen) = f u(&,n) do, 
cap 
ove # J. & la frontiera di J; su SM) la ®.(Je, 4), come funzione dell’intervallo 
Jz, per ogni fissato 4 in I, é VL in I;; 
2,) la variazione totale Vp (Zz, 7) &€ maggiorata da una funzione g,.(%) 
misurabile e sommabile in J,,. 


Le definizioni A e B, prima ricordate, rientrano in questa generale, ora data, 
come casi particolari relativi ai valori h=1 ed h=ry rispettivamente; si deve 
solo tener presente che, nel caso di # = 1, indicati con « e f gli estremi del generico 
intervallo J;, si ha ®,(J:,7)=(B, 7) —u(%, 7), mentre Vo(Zz,7) rap- 
presenta la variazione totale di w in J;. Inoltre, per h=r, la condizione 2,) 
non é da imporre perché perde senso?®. 

Si dimostrera (cfr. teor. XI) che, se w ed u% sono due funzioni equivalenti e 
verificanti la condizione 1,~m), le corrispondenti derivate iq), Oia) sono equi- 
valenti in J. Allora si assumera per definizione, come derivata di uw rispetto a 


1=x, la derivata ®j), relativa ad una qualsiasi delle funzioni equivalenti ad 


u e verificanti la 1.) € 2, ); tale derivata verra indicata con il simbolo Ee 


ath), ke 

L’equivalenza tra le diverse definizioni che si ottengono, facendo variare h 
da 1 ad 7 e, per ogni h, la combinazione «™) in tutti i modi possibili, pud dimo- 
strarsi, nel caso delle funzioni continue, estendendo la dimostrazione dell’equi- 
valenza tra [CVL], , ¢ [CVL],m,, Cul prima si é accennato!™. I] Prof. FICHERA 
ha posto il problema di dimostrare tale equivalenza nel caso generale (cioé senza 
fare l’ulteriore ipotesi della continuita) e di far vedere che tutte le derivate 

Ou 

Sr ie k 
dell’intero f e della combinazione «, differiscono tra loro solo nei punti di un 
insieme di misura nulla. 


che si ottengono, per un fissato k, in corrispondenza alle diverse scelte 


14 Si tenga presente che, in virtu del ben noto teorema di Fusini, tra le funzioni 
equivalenti ad una funzione misurabile e sommabile in / ve n’é almeno una la quale 
é misurabile e (h — 1)-dimensionalmente sommabile su ogni intersezione di I con le 
varieta (h — 1)-dimensionali di equazione: %,,=cost., ..., %4,_,= cost. t= cost. 

15 Si osservi che ®,w(Je,) non € altro che il primo membro della formula di 


GREEN: fue y) v,a@ =| oe dé. Poiché, per y=1, la formula di GREEN diventa: 
car iy) Je é. 
u(B, 4) —u(a, 7) -(3 dt, appare naturale assumere per ®,q)(Je, 7) l’incremento 
Je 


della w relativo all’intervallo unidimensionale Jz. 
16 Si tenga presente che, per h=yv, la ®yw(Je, 4) non dipende pit da 7; in tal 


caso pertanto useremo il simbolo ®, (Js). 
PONCE, MONE We 
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Nelle righe che seguono si daranno alcuni teoremi che risolvono la questione 
posta. sia 
Poiché le funzioni [VL],m,, sono tutte e sole le funzioni dotate di derivata 
debole rispetto ad «18, i risultati che verranno ora dimostrati permetteranno 
anche di stabilire l’equivalenza tra la definizione di funzione dotata di derivata 


debole rispetto ad x, e quella di funzione [VL], ; (HA, nasty FA). Dnelere 
si avra che la derivata debole locale di w1® coincide quasi ovunque con ls oT a’ 
In particolare quindi sara dimostrata l’equivalenza tra la definizione di funzione 
dotata di derivata debole rispetto ad x, e quella di funzione a variazione limitata 
generalizzata secondo CESARI-TONELLI (cfr. nota’). Quest’ultimo risultato € 
stato di recente dimostrato direttamente da K. KRICKEBERG®®. Dobbiamo anche 
osservare che l’equivalenza delle definizioni [VL], , e [VL],m,, poteva gia 
considerarsi acquisita combinando il risultato di FICHERA citato in#® con quella 
ora ricordato di KRICKEBERG. 

I. Sia u(P)=u(t,y) una funzione definita nell’involucro dell’intervallo I di 
S”, la quale sia ivi misurabile, sommabile, e verificht le 14), 240). Suppontamo 
inoltre che per ogni nEI,, la funzione di t: u(t, ny) sta continua ner due punti estremi 
di I, e sia continua a sinistra in ognt punto interno ad I,. Allora la u gode delle 
seguentt proprieta: 

a) per ogni tE I, la u(t, y) & funzione di n misurabile e sommabile in I,, 


b) indicata con V(n) la variazione totale di u(t,n) in I, e posto of=a+s ere 


(ove ae b sono gli estrema di I,) riesce: 


n 
V(n) =lim >) | a (al”, n) — wu (a1, m)|. 
Pe hell 
Come conseguenza di tal proprieta st ha che V(y) & misurabile e sommabile in I,. 
Dato che u(t, 7) € misurabile e sommabile in J, esiste un insieme N, contenuto 
in ‘ae di misura lineare nulla, tale che, per fc LW la funzione di 7: u(t, 7) é 
misurabile e sommabile in J,. Sia ¢ un punto di N—a; indichiamo con {t,,} 


una successione di punti contenuti in J,—WN tali che riesca ¢,,<?¢ e lim #,,=—1. 
m—-> CO 


Allora si ha: u(t,7) = lim u(t,,,7) e pertanto w(t, 7) risulta misurabile in 1, peu 
ie m—> OOo 

ogni ¢¢ J,;—a, perché limite di una successione di funzioni ivi misurabili. Inoltre, 

se ’EI,—N, riesce: | u(t,,, 4) —u(t’, n)| Sg4m(n). 


18 Cfr. [7] ove & dimostrato quanto segue: siano 4A un campo di S®, w una funzione 


localmente sommabile in 4, V la classe delle funzioni lipschitziane ed aventi supporto 
contenuto in A; se la u e dotata di derivata debole, cioé se per ogni v€V si ha: 


OM os . f : alee : uh 
u aan dx=—|vda, la funzione u é [J L)qin,% IM Ogni intervallo contenuto con 
. k 

A 


il suo involucro in A e riesce «(C) = [dO per ogni insieme chiuso e limitato C 

contenuto in A (ove con il simbolo eau si € inteso Vintegrale di RimMANN-STIEL- 
6 

TJES rispetto a Dy). Se, viceversa, la wu ¢ [VL]gm,4, esiste una @ equivalente ad u 

la cui corrispondente ®,) € una misura la quale verifica la relazione fuses 


= — [vd® per ogni v€V. Cfy. anche [138]. 
4 


** Per derivata debole locale di uv si intende la derivata della misura « (cfr. nota (17)). 
20 Cr. (10). 
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Ne viene: | «(t,,,7)| S| (¢’, m)| + g.m(m) e quindi la funzione u (t, 7) = Jim 4 (6,77) 
risulta sommabile in J,. In modo analogo si dimostra che (a, 7) é funzione di 4 
misurabile e sommabile in J,. 

E cosi provata la proprieta a). 

Indichiamo con 6 una arbitraria decomposizione di J; in intervalli?! e con 
6, la decomposizione di I,, ottenuta mediante i punti «{” (s=1,...,2—1). 
Se %p=A, %, -.., %,=, con %<---<u,,, Sono i punti mediante i quali é ottenuta 


la decomposizione 6, porremo a, (0) = >) |w(a,,) —u(%.—1,7)|. Per dimostrare 
s=1 


la proprieta b) ci basta far vedere che, assegnati ad arbitrio il numero positivo 
é e la decomposizione 6, esiste un intero ,() tale che o, (6) So, (dn_)) +e. 
Sia @ la minima ampiezza degli intervalli relativi alla decomposizione 0. 


Assumiamo ”,(7) in modo tale che riesca OEE < ed inoltre, per ogni# verificante 
cee & b—=@ 5 eA é ~ 
la condizione «,— tar <t<a,, si abbia |u(«,, 7) —u(t, n)| a E allora 


facile vedere che, con tali assunzioni, si ha: o, (6) So, (O,,(,)) + é- 

La misurabilita e sommabilita di V(yj) é immediata conseguenza delle pro- 
prieta a) e b) e del fatto che o,,(d,,) € non superiore alla funzione sommabile g,«) (7). 

II. Sia u(P) una funzione definita nell’involucro dell’intervallo I di S”, la 
quale sia 1vt misurabile e verificht la proprieta 1,~). Allora esiste una funzione 
u(t, ), definita nell’involucro di I, equivalente ad u, la quale, per ogni HEI, é 
funzione di t continua nei due punts estrenu di I, e continua a sinistra in ogni 
punto interno ad I,. Inoltre, indicata con V(n) (con V (n)) la variazione totale di 
u(t,) (di a(t, n)) in I,, st ha V(n)<V(y). 

Indichiamo con a, 0 gli estremi di J,, con J, un intervallo contenuto in J, 
avente per estremi i punti a, ¢, con v(t, 7) la variazione totale di u(t, 7), rispetto 
a ¢, nell’intervallo J, e poniamo v(a, 7) =0. 

Definiamo le due funzioni p(é, 7) e q(t, 7) al modo seguente: 

p(t,n) = 3 [v(t) + u(t) — 4 (4,9) 
q(t,n) = 3 (v(t) — «(t,n) +4 (4,m)]- 

Le funzioni f(t,”) e q(t,7) rappresentano rispettivamente la variazione 
positiva e la variazione negativa di w(t, 7) rispetto a ¢ nell’intervallo J,. Per 
ogni EI, ¢E1,, si ha: w(t, n) =u(a,n)+p(t,n) —@ (4), v(t 7) =P (6 0) +4 (6 0). 
Ne viene V (yj) =v(b,y) =£(2, y) +¢(0, 7). Le funzioni p(t, 7) e g(t, 4), per ogni 
7 €I,, sono non negative e monotone non decrescenti rispetto a ¢, come subito 
si verifica. Ognuna di esse, allora, presenta solo discontinuita di prima specie 
rispetto a ¢. Per ogni fissato 7, l’insieme dei punti di discontinuita di p(¢, 7) é€ 
finito o, al piu, numerabile. 

Consideriamo la funzione: 

ee a<t<b 
PET) ) gh Oe 
= lim p (7) pees 
21 Usando la locuzione ,,decomposizione di un intervallo in intervalli‘ intenderemo 


sempre riferirci ad una decomposizione di un intervallo nella somma di un numero 
finito d’intervalli a due a due privi di punti in comure. 
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E immediato constatare che A(t, 7), per 7 €1,, & funzione di ¢ continua nei punti 
t—aet—b, écontinua a sinistra in ogni punto ¢ interno ad J, ed € monotona non 
decrescente in I,. Riesce inoltre: £(b,) — p(4,) SP (8, 7). Indichiamo con 
g(t, 7) la funzione dedotta da q(é,7) con lo stesso procedimento con cul sl e 
dedotta # da p. Allora, la funzione u(t, n) =p(t,n) —G (ty), per ogni nEL,, € 
continua rispetto a ¢ nei puntit=ae tb ed é continua a sinistra in ogni punto 
interno ad J,. Inoltre, la a(t, 7), per ogni 7€J,, ¢ funzione di ¢ a variazione 
limitata in J, ed indicata con V(n) la sua variazione totale in J;, si ha: 


V (n) <p (6, y) —p(a,n) +4 (6, 0) —9 (a, 0) SP (0, 9) +4 (6,7) =V (n). 


Infine, avendosi: % (t, 7) = lim w (‘— =, n); perynel,, a<ts b, ed essendo x (t, 7) 


una funzione misurabile in J, si ha che tale risulta anche la #(¢, 7). D’altra parte, 
per ogni 7€J,, le funzioni u(¢,y) e %(t, 7) sono quasi ovunque eguali in J; e 
pertanto u(t,) e a(t, 7) risultano equivalenti in J. 

In tal modo il teorema é completamente dimostrato. 


III. Siano I un intervallo di S, {]} la famiglia degli intervalli di S” contenuti 
in I, F(J) una funzione additiva e VL in {J}. Esiste allora una ed una sola 
misura (J), definita in {J}, la quale gode della seguente proprieta : im corrispondenza 
au st puo determinare un insieme numerabile N; contenuto nell’asse x; (t=1, ..., 7) 
in modo tale che, per ogni intervallo J €{J} definito dalle condiziont: a; <x;<; 

Ny, Dp Ne C= 1,25. 7), Pesca. NJ) —= Lf). 

Pensiamo la F definita su tutti gli intervalli di S”, assumendola nulla in 
quegli intervalli che non appartangono a {/}. Sia w la misura generata dalla 
funzione f??. Indichiamo con M; (1 =1,...,7) V’insieme numerabile contenuto 
nell’asse x; e costituito dai numeri c; tali che l’iperpiano x;=c; sia di disconti- 
nuita per la misura 4“. Denotiamo inoltre con Mj (i =1, ..., 7) ’insieme numerabile 
costituito dai punti d; dell’asse x;, in corrispondenza ad ognuno dei quali é 
verificata la seguente circostanza: se (a,,...,4,) € (b,, ..., b,) sono i punti estremi 
dil, {/’}; ¢ la famiglia degli intervalli cosi definiti: a,,, = te< b, (mad, add = 
t+1,...,7), %€(a;,4,) allora, per ogni ¢ (¢=1,...,7), le funzioni di ¢;: pp(J’) 
e dr(J’) (con J’E{J'},) sono discontinue nel punto t;=d,. 

Assumiamo N,;=M;+M;. Sia J un intervallo contenuto in {J}, di punti 
estremi (a, ..., %,), (B,,...,B,), tale chea;¢N;, B;4.N;. Se J, ¢Vintervallo definito 
da: «;— “ S%,<6;+ = ( (1=1,...,7), riesce, com’é facile constatare; lim F(J,,) = 

n—->Cco 
F(]). Si he inoltre: u(J) =u(J). Poiché, d’altra parte, riesce: lim F(J,) =(J), 
si ha: F(J) =p(J). ore 


» IV. Stano u(P) una funzione definita nell’involucro dell’intervallo I di S ed 
wut [VL Jam, zp» A un intero positive minore dir ed «” una combinazione contenente k. 
Indichiamo con Jn} la famigha degli intervalli di S’—" contenuti in I, e con J 
Vintervallo di S” definito dalle condizioni: Lely Vel ye Jan eee shbore 
1a Es) =Vo.m(J) & assolutamente continua in WA 


8 Dicendo che w é la misura generata da F(I) intendiamo che, per ogni insieme 
C chiuso e limitato di S, riesce: u(C) = J dF, Vintegrale a secondo membro essendo 


considerato nel senso di RrEMANN- Sater 
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Supponiamo che la funzione wu stessa verifichi la condizione 1,4). Indichiamo 
con 6 una decomposizione di J, mediante i punti a, 0%, ..., %,, CON by <0 <<a. 
La famiglia di tali decomposizioni sara indicata con {6}. 


Fissata una decomposizione 6€ {6}, per ogni J, €{J,} poniamo: 
F.(Jn) ae dn J [u(%s, 9.7) — 4(a1, 90m) ) dy, 
y 


n 
e inoltre: 


Van) = S| |u(a;, 7,4) — w(%.1, ¥,n)| dy] an. 


Indichiamo con /, ; Vintervallo di S definito dalle condizioni: «,_,<t<«,, 


Used, sanvdnay: 
(2,) Vr, VJ,,) = Von ls, J) x 


Infatti riesce: Oy (I, j,) =4s(J;,); allora, indicata con A una decomposizione 
di T,,3, in intervalli di Si): I, ,=0+-+-+1™ e posto o(4) =¥ |®a0() si 
ha che: Vo,,,(Is,,) € Vestremo superiore di o(4) al variare comunque di A, 
mentre V;(J,,) coincide con l’estremo superiore di o(4) al variare di A soltanto 
nell’insieme delle decomposizioni di Is. jp del tipo ee ee con 
T= (#51, %). 

Cio prova I’asserto. 

Sommando membro a membro le (2,) per s=1,...,, si ottiene 


n 


Modi dle 2 Ve @yn (Ls, Jp) = Vogm J) =¥% Ur) 
e quindi: 


@) S| Lees. 9.) — ea, 900)] dy] an S 1h) 


Si Jy Ty 


Il primo membro di tale diseguaglianza, per ogni fissato J, €{J,}, ¢ funzione 
di 6. Sia Y(J,) il suo estremo superiore su {6}. Facciamo vedere che riesce: 


(4) Yn) =V Uh) 
Dalla (3) si ha intanto: V(J,)<V(J,). Inoltre, fissata una goon 


A di J in intervalli di S”: J=M4+--.-+1™ e posto: o( = |Dan(L* 
sempre determinare una 5 tiny Sa 6 di {6} tale ie ae 


n 
ASX S| flues, yn) —u(aaa, vo dy] dni 
cat, dy 
as per questo considerare la totalita degli iperpiani d’equazione t=«, 
(= , n) ognuno dei quali contiene almeno un punto di JI ‘in corrispondenza 
a a tie i (¢=1,...,m), ed assumere come 6 la decomposizione di J, relativa 
ai punti a%,...,%,,- In tal modo é provata la (4). 


5* 
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Sia m) un punto di J,, ed indichiamo con I,,,1 bintervallo quadrato contenuto 
in I, che ha 7) come estremo inferiore ed il lato di lunghezza /. Dalla (3) si trae: 
y 


~ 1 
65) Soe [| [le e.9 0) —m (ea mI i] dn SV ln): 
s=1 Ines 
Indichiamo con N, l’insieme contenuto in J, e di misura (7 — h)-dimensionale 
nulla, tale che, per 7 € J,,—Nj, Tlesce 


: 1 
jim oy [| f belo. 2o0) Wea enldy 
Ty 


I 


dy 


Nor? 
= f ls. 910) =U (45-1, 9, Mo) | @y- 
Ty 
Supponiamo anche che, per 7€J,—Nj, esista la derivata Y'(y) di Y(J,). 
Passando al limite nella (5), per +0, in corrispondenza ad ogni punto 7 CI, —N,, 
si ha: 


> f|w(as, 99) — 4(a%-1,¥,9) [49 SV'"(n), 


s=1 I, 


donde, per ogni J, € {J,,}: 


Sf [ Slee. 90) — 4-2. nallay] ans [1 dy. 


s=1 Jn Ty Jn 
Da questa diseguaglianza si trae: Y(J,)< [V'(n) dy. Poiché Y(J,) =V(J,), 
si ha l’assoluta continuita di W(J,). Jn 


V. Sia u(P)=u(t, vy) una funzione definita nell’involucro dell’intervallo I di 
SM ¢ [VL], , 1m I. Esiste allora una ed una sola misura «a(J), definita nella 
famiglia {J+ degli intervalli di S® contenuti in I, in corrispondenza alla quale é 
possibile determinare una funzione u, equivalente ad u in I, tale che: 


Pan (J) = findo =a(J). 


IJ 
Supponiamo che la w verifichi la 1,). Indichiamo con «(/J) la misura 
generata in {J} dalla funzione ®,»)(J) = [w»,d@. Per il teor. III, esiste, per 


ogni 7, un insieme numerabile N,, homens nell’asse x;, tale che per ogni inter- 
vallo J €{J} definito dalle condizioni a;<x,<b;, con a;€1,,—N,, ;€1,,—M,, 
riesce: ®y(J) =a(J). E immediato constatare che, in questo caso, per ogni 
74=k, Vinsieme N, € vuoto. Siano a un numero contenuto in J ,—N,, R un inter- 
vallo contenuto in J,, avente un estremo nel punto a e l’altro in un punto 3, 
{H} la famiglia degli intervalli di Sk contenuti in J,. La funzione di H: 
a(kxH) é assolutamente continua in {H}. Infatti, dato che, per il teorema 
precedente, la funzione Voyn(2t XH) & assolutamente continua in {H}, indicata 
con g(y) la derivata della funzione di H: Voyn(lixH) e con {b,} una successione 
di numeri tali che: jim b, =, 6,€1,—N,, 6, <0, si ha, per ogni H€{H}: 

H 


J [4 (On. ¥) — u(a, y)| dy =e) dy. 
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Poiché: 


’ 


lim flu@n.») — ula, »)] dy = |a(k x H) 


n—> CO 


si ha: 
|a(RxH)| SJg(y)dy, 
donde l’asserto. i 
Esiste allora una funzione di y: /(b, y), misurabile e sommabile in J y> tale che: 


= J f(b,y) dy b>a 
AH 


a(R x H) 
aie Oa. 


Poniamo, per (¢, y) EI: 
=fGy)+u@y) telon,. 
Si ha allora, per ogni intervallo J € {/}: 
fivdo=a(J). 
Soap 


L’unicita della misura che gode della proprieta ora dimostrata per « segue 
immediatamente dal teor. III. 


VI. Sia u(P) una funzione definita nell’involucro dell’intervallo I di S. Se 
uw é [VL }ym , essa & anche [VL], 4. 

La funzione w, per ipotesi, € equivalente in J ad una funzione verificante la 
{,«)). Per semplificare le notazioni, supporremo che la wu stessa verifichi tale 
condizione. Siano ae b gli estremi di J,. Indichiamo con J, l’intervallo di S 
che ha come estremo inferiore il punto a e come estremo superiore il punto ¢ 
(ove a<tS<b) e con {J,} la famiglia di tutti questi intervalli; denotiamo inoltre 
con {J,} la famiglia degli intervalli di SY~” contenuti in J, e con {J} la famiglia 
degli intervalli J =J,x J,, ove J,E{ fi} e LE {IS}. 

Diciamo fo (J) © Yo,)(/) vispettivamente la variazione positiva e nega- 
tiva di Di in oe . 


Poni ; 
ee wy, py [= Pol) se a<tsd 
| == 0 se t =a 
= Veent)} se a< tb 
V(t gee 
OB See — a 
aatpat (),) se a<tSb 
t (A) 
pt} =? Sse. t= a 


J) Seat -at SD 
= () set=a. 


23 Dato che ®y@ é una funzione additiva e VL, le sue variazioni positiva e negativa 
sone definite dalle relazioni: 


POyn(J) == = (Dyn(J) air Voyn(J)] , I@4r(J) = 2 (Vom) oe Dyn(J)] 
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Si ha, per ogni edie 
(6) p(t, i) =4(O6 5) +VEF)), hh) =2 [Vie f) — Bt, J)]- 


La funzione di J,: P(t, J,) & assolutamente continua in { Jy per ogni fissato ¢ 
in I,, come segue dalla sua stessa definizione ; VG dye funzione assolutamente 
continua di J,, come risulta dai teoremi precedenti. Dalle (6) si ne allows ae 
per ogni ¢€1;, le p(t, Jj), g(t, f) sono funzioni assolutamente continue in {/J,}. 

Dato che fo,»(J) & additiva e non negativa, per ogni LeU}; pe I) 
funzione monotona non decrescente di f in J;. 

D’altra parte esiste un insieme NCJ,, di misura (ry — 1)-dimensionale nulla, 
tale che, indicato con J, ; Vintervallo quadrato di cS al contenuto in I,, avente 
il punto y come estremo inferiore ed il lato di lunghezza J, per ogni Ve 
riesce : 

lim n’~" p (0, Iy,1) << + 0. 


n—> Co ‘e 


Si ha allora, per Te Vie ave 
max lim n’~1 9 (¢, Iy,1) S lim n’™* (8, Ty,1) < + 0. 
n—> CO n n—> co 


Per ¢€ I, poniamo: 
{ = max lim n’* p(t, I,,3) se yEl,—N 


n—-> CO 


Bb; ‘ 
Neg se yEN. 


La A(t, y), per ogni yEI,, € funzione di ¢ monotona non decrescente in ie 
infatti, se yCI,—N e >, si ha: p(t, Iy,1)Sp(", I,,1), donde: 


1 
"nN 


max lim n’~* p(t, Iy,1) S max lim dilins 2a Used bo 


noo 


Inoltre A(t, y) € non negativa dato che tale ¢ p(t, Jy). Poiché p(, J,) € funzione 
assolutamente continua di J,, per ogni ¢€J;, si ha: 


ten 
eek a 


Facciamo vedere che A(t, y) &é misurabile in J. Cominciamo, per questo, con 
losservare che la funzione di t: #(¢, Jy), essendo monotona non decrescente in J,, 
puo essere discontinua soltanto in un numero finito 0, al piu, in un’infinita 
numerabile di punti di J,. Indichiamo con N, Vinsieme di questi punti. E evidente 
che anche la funzione di t: p(t, J,), fissato comunque Jy in {J,}, pud essere dis- 
continua soltanto nei punti di N. Poniamo A(t, y) =~" pit, I,1). La Alt, y) 


é funzione di y continua in I, uniformemente rispetto a t, come segue subito 
dalla diseguaglianza : 


(7) b(t, Jy) =I Pit, y) dy 


PU h)SSPby)dy tel, 

Jy 
Ne viene che, indicato con N’ V’insieme contenuto in I e costituito dai punti 
(t, y) per ognuno dei quali si ha: f€ MN, la funzione F,(¢, y) @ continua in J —N’. 
Poiché N’ ha misura r-dimensionale nulla, si ha che F(t, y) & quasi continua 
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e quindi misurabile in J. Poiché riesce: 


= max lim #1 (?, y) yel,—N 
Pty9){ ee oN 


resta provata la misurabilita di Fit, y) in I. Ripetendo per Jo,m(J) le con- 
siderazioni ora svolte per bo,(J), Si perviene a dimostrare l’esistenza di una 
funzione 2(t, y), definita per ¢€J, e y€I,, misurabile in J, la quale, per ogni 
ycéJI,, € funzione monotona non decrescente di ¢ in i é nulla in a, é funzione 
di y misurabile e sommabile in I, per ogni t€ I, ed é tale che: 


~ tel, 
(8) q(t, Jy) =p eres (eles 
Per ogni J €{/J} si ha: 
PD. (J) = Po (I) — doun(J) =P 4 J) — 2b J), 
eloe, per (7) ed (8): 


Dyn (J = J la( t, vy) —u(a,y)|dy red [Pt,v) —2hy)] dy 


y 


tel, 
HEADS: 

Ne viene (eastais di una funzione w(é, y) definita in I, ivi misurabile, equi- 
valente ad w in J, tale che, per ¢€/;, yEJ, si abbia: 


w(t,y) —u(a,y) =PAlt,y) — 2it,y). 


Pertanto la funzione w(t, y), per ogni yCl,, € VL rispetto a ¢ in J, e per la 
sua variazione totale V(y) si ha V(y)SA(d, y) + 2(0, y). La funzione A(d, y) + 
2(b, y) &é misurabile e sommabile in [,,. Ne viene che la funzione w, equivalente 
ad u in I, soddisfa le 1,«), 2,«) e pertanto u é [VL], ,- 

Cosi il teorema € completamente dimostrato. 

VII. Szano h un intero tale che 1Sh<r, k un wntero tale che 1SkSyr, u(P) 
una funzione definita nell’involucro dell’intervallo I di S”. Se u é una funzione 
[VL] um, , essa & anche [VL] x, ¢- 

Considereremo il caso h>1. Per h=1 la dimostrazione é perfettamente 
analoga a quella che esporremo ora. 

Supponiamo, per semplicita, che la ~ stessa verifichi le 1,~@), 2.m). 

Sia N un insieme contenuto in J,, di misura lineare nulla, tale che, per 
t€I,—N, la funzione di (y, 7):u(¢, y, 7) sia misurabile ‘ sommabile in I, <T,. 

Indichiamo con {J} la famiglia degli intervalli J di S” contenuti in J, ognuno 
dei quali é tale che J, non ha gli estremi contenuti in N. Per ogni J €{/} risulta 
allora definita la funzione Oy»(J) = f u(t, y,) 4d. Mostriamo che ®, (J) 

df 


é VL in {J}. Sia J un intervallo contenuto in {J}. Diciamo 6 una arbitraria 

decomposizione di J: J =J1+---+J? e poniamo a(0) =) | Dyn (J). 
Consideriamo inoltre una decomposizione di /.: = Iz+---+If e una 

decomposizione di Ip=fit-  +Tq. Riesce: j=. iy MPa Indichiamo con 


6’ una tale decomposizione di J. Fissata arbitrariamente una decomposizione 6 
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di J, puod ovviamente scegliersene una, 6’, tale che o(6)<o(0’). Si ha d’altra 
parte: 


o(0)=>¥ YiPolExMl=L DV) Sub sn) do 
i=1 s=1 i=1 s=1 |! Y(IExI}) 
Dag ps if Pam! (It, ) dn|< Jf > |G, (Le, n)| dn SJ 8m (7) dn. 
ood Jj 41 E 


E cosi provato a Dn} ) eV Lan J}. Allora, con ragionamenti perfetta- 
mente analoghi a quelli dei teor. III e V si prova l’esistenza di una funzione 7, 
equivalente ad wu in J e verificante la 1,1). 


VIII. Se u é una funzione definita nell’involucro ee I di S”, ivi 
[VL] a, ,, ndicata con {]} la famigha degl intervalli di S contenuti in I, esiste 
una junzione w equivalente ad u in I la quale, scelti comunque l’intero h e la com- 
binazione «™, verifica le 1,m), 2m) ed é tale che, per ogni fissato n in I, la 


Pm (Jen) = fondo 
SL JE 
é una misura; inoltre la corrispondente variazione totale: Vo y,(Lz,7) & musurabile 
e sommabile in I, € riesce quast ovunque: 


J Vog(les 90) 49 = Vogoy(le.m) =¥"(H), 
ove con V'(n) st é indicata la derivata della funzione dell’intervallo J, : 


V(S;) = Von (I. so) : 


Per i teoremi VI e VII possiamo supporre che la w stessa verifichi le 1,«)), 2,«). 

Per i teoremi I e II esiste una funzione u“, equivalente ad wu in J, la quale, per 
ogni fissato t in J,, € continua rispetto ad x,=¢ nei punti estremi di J, ed é 
continua a sinistra in ogni punto interno di J;. Inoltre la sua variazione totale 
V(x) & misurabile e sommabile in J,. 

Siano a e b gli estremi di J,. 


E facile vedere l’esistenza di un insieme N, contenuto in J,, di misura (7 — 1)- 
dimensionale nulla verificante le seguenti condizioni: 


1°) indicata con {J"} la totalita delle varieta lineari ognuna delle quali é 


definita da equazioni del tipo x,, =cost., ..., x;,_,=cost. (ove b,,...,b,_, & un’ar- 
bitraria combinazione dei numeri 4, 2,...,k—1,k+4,..., vr) e posto: 
V(r) = V(t) tel, —N, 
=0 tEN, 


la V(x) risulta misurabile e sommabile su ogni intersezione di J’ con l’insieme 
definito da: t=a, t€J, al variare comunque di J’ in {T}, 
2°) la funzione cosi definita in T: 
alt, 7) == (0) tel, —N, 
=) TEN, 


é misurabile e sommabile su ogni intersezione [J al variare comunque di I’ 
in {I}. 


i i 
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Sia f un intero compreso tra 1 ed 7 ed «”) una combinazione di classe h dei 
numeri 1, 2,...,k—1,k+4,...,7. E facile constatare che la u(t, t) verifica la 
1,q@). Siano 4 un punto di J,, 6 una decomposizione di J, relativa ai punti 
Dar latin chy Ot Ol ha: 


n? 


Va ies 7) =estr. (sup. & Sf laa, yn) —@(0,-4, yn) dy, 


s=l1 ly, 
La dimostrazione di questa relazione € analoga a quella con cui si é provata la (4). 
Poiché si ha 


> |# (5, 9.9) — #(a,-1,9,0)| S Vy, 7) 


s=1 


n 
estr, sup. DAAC Y,1) — 4 (a1, ¥,)| = Vo (¥ 0)» 


riesce: 


J Vo y(¥> n) dy = Vole, n) ; 
y 


Resta cosi anche provato che Vg,,(Z;,77) € funzione di 7 misurabile in /,. 
Inoltre, da (4), si trae: 


estr. soup. 2 i |S | (a,,¥,) —#(a,-1,¥,)|dy|dn =V(J,). 
7) y 
Poiché: 


cotta sup. > f la CO i} = u(x As_4, V, 1] )|dy = VE yams 1 1)» 


s=l1ly 


tenuto conto della assoluta continuita di W(J,), si ha: 


I Voyy(Le.1) dn =VU,) = [V"(n) an. 

UY] UI 
Sia N, ae contenuto in J,, di misura (y —/)-dimensionale nulla, tale che, 
per 7 é dy cau, phlesce:s 


V5 ales é n) =V¥'(n). 


Indichiamo con N linsieme, di misura y-dimensionale nulla, costituito da tutti 
e soli i punti che es della seguente proprieta: se «”) é una combinazione di 
classe i dei numeri 1, 2, ..., 7, per ogni punto (&, 7), contenuto in N riesce 7 CN. 
La funzione w(P) cosi definita in I: 
w (P) = (P) se PEIN 
=0 se PEN 


gode di tutte le proprieta ora dimostrate per # ed inoltre € tale che, fissato 
comunque «”), riesce quasi ovunque: 


Voxu(te>) rob 


Resta soltanto da provare che, ve ogni a”), la Daa J: m4) € una misura nella 
famiglia {J,} degli intervalli di S{” che sono “contenuti ine gs 
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Dai teoremi IV e VI discende che, per ogni a), esistono una misura [Mair (Je ; ”), 
definita, per ogni 7, nella famiglia { J-} e un insieme N, 2(n) verificanti le seguenti 
proprieta: 

1) N2(n) & contenuto in J, ed ha misura lineare nulla; 

h 
2) per ogni intervallo f,: yEJ,, aSt<b con a€l,— N2(n), 6€1,—N2(n), 
riesce: 44» (Jz) = Pym (Je, 7)- 

Sia ora J, un intervallo di { Ji} } definito dalle condizioni: y¢ J, eye Je ae 
intervallo di Sa 1) contenuto in I,) e a’St<b’, con @€l,— NZ“(n), pel 
Se {¢,} & una aah: di punti di I, tale che@ =<4,<0, £6l.—Dige 
im ¢, =b' ese Ie” é l'intervallo definito da: ye J,, a’ Sy<t, si ha: 


lim pia (Jé”, 4) =Ham(Jen) € lim ®,m Je”, n) = Peo (Je,)- 
N—> CO n—> 
Ne viene la coincidenza di 0) con Om, su tutti gliintervalli di S ” contenuti in J,. 


In tal modo il teorema é completamente dimostrato. Sia A un campo di 
S. Diremo che la funzione u(P), misurabile e localmente sommabile in A 
(cioé sommabile in ogni J tale che IcAé [VL], , nel campo A se esiste una 
funzione equivalente ad u in detto campo la quale verifichi le condizioni 1,«), 
2,0) in corrispondenza ad ogni I tale che <A. Fissata una combinazione «", 
indicheremo con A, l’insieme di tutti i punti ¢ di S) in corrispondenza ad 
ognuno dei quali l’iperpiano d’equazione t =f, ha intersezione non vuota con A. 
In modo analogo definiremo gli insiemi 4,, Az, A,, Ay, etc. Denoteremo 
inoltre con , (t)) !intersezione di A con l’iperpiano d’equazione tf). In modo 
analogo si definiranno gli insiemi -¥, (7), SF, (€), L4(y). Indicheremo infine con 
{I}, la famiglia degli intervalli di S contenuti con il proprio involucro in A, 
con {J,}, quella degli intervalli di Ss contenuti con il proprio eyes in 
se, (n) per 7€A,. Analogamente si definiscono le famiglie {Jz}, {Ii}, € 


Dimostriamo il seguente teorema: 


IX. Sia u(P) una funzone definita nel campo A, ivt misurabile e localmente 
sommabile. Condizione necessaria e sufficiente perché u sia [VL], in A é che 
essa sia [VL] am, , 1 ogni intervallo I< {I},. 

La necessita della condizione ¢ ovvia. Dimostriamone la sufficienza. Sup- 
poniamo 1<h<vy. Indichiamo con {J} una famiglia costituita da un’infinita 
numerabile di intervalli di {7}, a due a due privi di punti in comune, tali inoltre 
che la loro riunione coincida con A e che ogni intervallo J di {I}4 abbia inter- 
sezione non vuota soltanto con un numero finito di intervalli di LT “Const 
teoremi VI, VII, e VIII si é fatto vedere che, se u é [VL], , in un intervallo I 
di S™, si pud costruire una funzione a(t, y,) equivalente ad u in J, la quale 
soddisfa le condizioni 1,@) € 2,q)), ed inoltre, fissato ¢ in I,, risulta misurabile e 
sommabile in I, xJ,,. Applicando in ogni intervallo I di {I } il detto procedimento, 
si perviene a tesuaties una funzione w equivalente ad wu in A, la quale gode delle 
seguenti proprieta: scelto comunque un intervallo I di {7}, esiste una funzione w, 
definita in J, coincidente con @ in I, verificante le condizioni 14) © 2ya) rela- 
tivamente ad J e tale che, fissato tin I,, w(t, y, 7) risulti misurabile e sommabile 


in I, x1,; inoltre, per ogni ¢€ I,, la w(t, y, y) & essa stessa misurabile e sommabile 
in L, x I, e verifica la 14m) per ogni I €{I}4. 
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Facciamo vedere che la funzione @ verifica le condizioni 1,0), 2a) in ogni 
intervallo I di {I}. Sia I€ {I}. Indichiamo con ae 0 gli estremi di J, e con {Je} 
la famiglia degli intervalli di S” ognuno dei quali é contenuto in J ~- Poniamo 
Dm (Je, 7) fe J w(&,)%d@ ed indichiamo con g,m(y) una funzione misurabile 


e sommabile in I, tale che Voz, 7) S8.m(). Poniamo, per ogni 7 €I,,: 


P43) VJe,n) =fulé,n)y,do. 
J Je 


Se J; é un intervallo appartenente a {J,} lacui frontiera non ha alcun punto 

nell’iperpiano d’equazione t=b, riesce: O,w( Je.) =P ym (Je, 4). Indichiamo 

con 6 un’arbitraria decomposizione J, in intervalli mediante i punti a=a<a,<-:: 

<a%,=b ed indichiamo con ff l’intervallo di S{ definito dalle condizioni: 

a, ,~St<a,, y€l,. Poniamo 

o(0) = 2 J Hs, 9, 7) —# (052, 95 0)| dy 

Sina : 
oO) 22 Jlwles, 9m) — wa 9.0) dy +S |, 9,9) — #(b.9,m)| dy. 

Gaul, Ihr y 


Da tale diseguaglianza si trae: 
V5,m(Le>) S Voym(Ze,1) + fle, vn) — 4b, y,n)| dy 
Ty 


e quindi: 
Vole, n) = So.(%) (7) si if | w (0, y; 0) os u (0b, Ve n)| dy. 
Ty 
Poiché il secondo membro di tale diseguaglianza é una funzione misurabile e 
sommabile in J,,, resta provato che # verifica le 1,~), 2,0) relativamente all’inter- 
vallo I di {I}. E allora evidente che # verifica tali condizioni anche in relazione 
a qualsiasi intervallo I € {7},. 


X. Se u é una funzione misurabile e localmente sommabile nel campo A ed tvi 
[VL am, x, eseste una funzione w equivalente ad u in A la quale, fissati comunque 
Vintero h' e la combinazione a), verifica le yw) € 2ya) im corrispondenza ad 
ognt IE{I},, inoltre é tale che, per ogni fissato yn in A,, la Oy wy) (Le,) & una 
misura in {I}, e, fissato comunque I in {I}4, quasi ovunque in I, riesce: 


4 Voy (Les V0) 9 =V"'(n) = Vo, ay(Ze, 0)» 
y 


ove con V'(y) si é indicata la derivata della funzione di J,: Vo.,LeXJ,) nella 
famiglia degli intervalli J, in corrispondenza ad ognuno det quali riesce: Ie x J, EC {T} 4. 

Infatti, per i teorr. VII, IX, la uw é [VL], in A. Possiamo allora supporre 
che la wu verifichi le 1,«)) € 2,«) in corrispondenza ad ogni J € {I}4. Indichiamo 
con {I ‘)} una successione di intervalli di {J}4, a due a due privi di punti in 


comune, tale che: A = Ul ‘) e tale che ogni intervallo I€{I}4 abbia prodotto 


non vuoto soltanto con un numero finito di insiemi di {J}. Per ogni s indichiamo 
con uw) la funzione definita in J) e dedotta da u con il procedimento indicato 
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nella dimostrazione del teor.II. Siano (a, saiitaS ine Oe —_ i punti 
estremi di I), Consideriamo la funzione # la quale coincide con uw) per ogni 
Pinus, ey) tale che a’)<x«x,<b*). La funzione w, per ogni (y, 1) EA, 9» 
riesce continua a sinistra rispetto a ¢ in ogni punto di Y4(y,%). Sia we la 
funzione definita in I) e dedotta da % con il procedimento indicato nella dimo- 
strazione del teor. VIII. La funzione w che coincide con w in I Cr eae res 
soddisfa le condizioni richieste dal teorema, come é immediato verificare. 


Vogliamo ora occuparci della definizione di derivata rispetto ad x, per una 


funzione che sia [VL], , € dell’equivalenza delle derivate al variare 


A c Ox, ah), k 
in tutti 1 modi possibili della combinazione op). 


Facciamo innanzitutto vedere che: 


XI. Se u(P) é@ una funzione definta nell’involucro dell’intervallo I di So), ivi 
misurabile e sommabile, e verificante la condizione 1,m), st ha che la derivata 
P.w(P) della corrispondente funzione d’intervallo Bw (Iz, 4) & maisurabile in I. 
Inoltre, se u e un'altra funzione verificante le stesse vpotesr fatte per u e tale che 
u=U quasi ovunque in I, st ha che la derivata Oe della corrispondente funzione 
d intervallo D 4% (Iz, 7) & equivalente a Din (P) in I. 

Siano a e b gli estremi di J,: Indichiamo con J, Vintersezione di J, con 
Vintervallo quadrato di Stee! che ha l’estremo inferiore nel punto y di J, ed il 
lato di lunghezza 7. La funzione y,(P)=q,(t, v, 7) =J u(t, y’,9) dy’ & misu- 


. . . y,t Ce 
rabile in J. Basta evidentemente limitarsi a dimostrare cid nel caso che u(P) 


sia non negativa. In tale ipotesi esiste, com’é noto, una successione non decre- 
scente {u, (t, y, 4)} di funzioni misurabili e costanti a tratti, convergente ad w(P) 
in I. Siano L®,..., L°” gli insiemi in ognuno dei quali u, @ costante. Diciamo 
Wi(P) (i=1,..., m,) la funzione nulla in I—L ed eguale ad 1 in L™; se al é 


n ? n 


The oe : 
il valore assunto da u, nei punti di L, si ha: u,,(P) =» a) YW (P). La funzione 
1 


J BO, y’,n) dy’ & misurabile in I; ne viene che tale é anche la funzione di P: 


wl 
Jun (t,y',n) dy’. Poiché f u(t, y’,n)dy'=lim fu, (t,',y) dy’, si ha Vasserto. 


nN—>OO 7” 
Ty,1 Iy,1 Ty, 


Da cio si trae che la funzione di P: y,(P) ¢ misurabile in J. Sia E l’insieme dei 


punti di J nei quali non esiste finito lim nq (¢ sari) 1] — g(t, t)|. L’insieme E 
n—> co n } n 


é misurabile. Ne viene che la funzione ®j)(P), la quale coincide con tal limite 
nei punti di J — EF e vale zero nei punti di E, é misurabile in J. Sia N un insieme 
contenuto in J,, di misura (vy —h)-dimensionale nulla, tale che, per ele 
Vinsieme dei punti di 7, in ognuno dei quali w(&, 7) «(€, 7) abbia misura h- 
dimensionale nulla. Se 7 é un punto di J ,— N, indicheremo con N(v) un insieme 
contenuto in J;, di misura lineare nulla, tale che, per ¢€1,— N(y), Vinsieme dei 
punti di /, nei quali w(t, y, n) == “(t, y, 7) abbia misura (h -— 1)-dimensionale nulla. 
Sia 7 €I,—N; indichiamo con { Je}, la famiglia degli intervalli J. di SY) ognuno 
dei quali é definito dalle condizioni yé Jy, ove J, & un intervallo di S“~ tale 
ches/pOi Fl), "ast=6 con eel (N(m) +a), BEL,— (N(y) +0). E imme- 
diato constatare che la famiglia {Je}, costituisce un semianello tale che la famiglia 
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dinsiemi totalmente additiva minima che lo contiene é quella dei boreliani 
aventi il proprio involucro in J,—YI,. 
Si ha, ovviamente, per ogni 7 ¢J,—WN e per ogni Ten Jade 


(9) ; ©, Je.) = Pam (Je) - 
Consideriamo la decomposizione di LEBESGUE di ®,) relativa alla famiglia { Jia 
eo) P,0(Je1) = Oy Jes) +f Bey (E,n) a8, 

Je 


ove Dxw(Je,7) & Ia componente singolare di ®,a(J:,7). Da (9) e (10) e dalla 
analoga relazione scritta per Di, si trae, in particolare, per ogni 4€J,—N, 
IE wiry = 

a Dy (E, y) d& =i Drm (E, 7) dé. 


Ne viene che, fissato 7 in I,—N, per ogni boreliano B; di S contenuto, con il 
proprio involucro, in J,, deve essere: 


J Pane, n) dé = J Pam| &,7) dé. 


Se ne conclude, che, per 7¢J,—N, quasi ovunque in I; riesce: Dyn (E, 4) = 
Dj (€, n). Da cid segue l’asserto. Con dimostrazione analoga a questa si provano 
1 teorem1: 


XII. Se u(P) é una funzione definita nell’involucro dell’intervallo I di S™, ivi 
misurabule, sommabile e verificante la condizione 1, «), allora la derivata ae. é 


misurabile in I. Se u(P) & una funzione verificante le stesse tpotest fatte per u, 


ed equivalente ad u in I, st ha che Fide i ee quasi ovunque in I. 
OX, OX, 


XIII. Se u(P) ed u(P) sono due funzioni definite nell’involucro dell’intervallo 
I di S, ivi misurabili, sommabili e verificanti la condizione 14m), e se u(P) ed 
a(P) sono equivalenti in I, si ha ®yn(P) = Pyn(P) quasi ovnnque in I. 


Questi teoremi, come gia si é detto, consentono di definire la derivata rispetto 


ad Xp: he 


aa di una funzione w che sia [VL], , in un intervallo J di so 
k 


ox(), Fe 


Dimostriamo che: 


XIV. Se u(P) é una funzione definita nell’involucro dell’intervallo I di S la 
quale sia 1vt [VL] ym, , estste una funzione u equivalente ad u in I, la quale gode 
della seguente proprieta: per ogni (y,)Cl,xI, ¢ VL rispetto a t im I; la sua 
derivata ©" & misurabile in I e riesce: BS | : 

OX, OX, OX, lal, k 

Per h=1 il teorema é evidente. Supponiamo 1</A<yr. Siano ae bd ghi 
estremi di J,. Denotiamo con {J.} la famiglia degli intervalli di S contenuti 
in J;, per ognuno dei quali si ha: J,=J,x J, ove j, € lVintervallo di S™ che ha 
estremo inferiore a ed estremo superiore ¢, ed ove J, € un intervallo di Sie 
contenuto in J,. Consideriamo la funzione w, equivalente ad u in I e dedotta 


dalla « con il procedimento esposto nella dimostrazione del teor. VIII. Deno- 
tiamo con Po (Je, 1); Gon Je. 1)» Vesa Se») rispettivamente le variazioni 
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positiva, negativa e totale di Oy (Je, 4) = fw(E,)%do in J,, per ogni 4 €1,. 


Poniamo, per 7€1,, Je€ {Je}: Pa 
aerial pees si 
Va) = 6 aan Jes) alia 
ae ee 
an) 6 qa (Jz 1) ee 


Sia I, , lintersezione di J, con l’intervallo quadrato che ha come estremo inferiore 
il punto y di I, ed il lato di lunghezza /. La funzione di P: ®(, ;, t, 4) € misura- 
bile in J, come osservato nel corso della dimostrazione del teor. XI. Proviamo 
ora la misurabilita di V(J, ;, ¢, 7). Innanzitutto é facile vedere che se a<t<b, 


indicato con s(¢) il piu grande intero non negativo, tale che: a+s(t) - bel 


riesce: 
s(t) 


(11) Vd,,047) =lim >) al a a,”), y’, 9) — was, y’,n)| dy’. 

n— CO S=1 y,1 
Infatti, rpetendo un ragionamento analogo ad uno fatto nel corso della dimo- 
strazione del teor. VIII, si vede che, indicato con {6,} Vinsieme delle decompo- 
sizioni di (0, ¢), si ha: 


Vig tn) = a oe a |w(a,,¥’,7) — w(a,4,¥',y)| ay’. 
s= yt 
Poiché w(t, y, 7) & funzione di ¢ continua nei punti a e b e continua a sinistra 
in ogni punto interno ad J,, é facile dimostrare che 


Mm 


ety. soup 2 f |w( ay ea %—1,¥'.9)| ay’ 


s=1 Iy 1 


= im >‘ J | elas" 9,9) — wo, 9’, n)| dy’. 


N+ OO 5 lel 


Resta cosi provata la (41). La funzione 


walls 9m) = 3. f |o(al, y',9) — tw (als, 9 mi ay', 
s=1 Ty} 
per ogni fissato ¢ in J,, &€ misurabile rispetto ad (y, 7) e per ogni (y, nel Ris 
e aieion di ¢ costante a tratti (precisamente é costante in ognuno degli intervalli 
”) <t<o”). Pertanto Pn(t, y,) € misurabile in J. Dalla (11) segue che tale 
Cail anche la V(I, ;,¢,). Si ha per ogni nel, te Ty, eau ae 


P(J,,4n) =F(OU,,tn) +VF,,6n)] 


(12) 
q(Jy.t) =3(V(,,4) — GJ, t,n)]. 
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Per ogni y €I,, ¢€1,, le funzioni di Js D(J,, t,y) e V(J,, t, m) sono assolutamente 
continue in {J,} come segue dalla definizione stessa di © e dal teor. IV. Dalle 
(12) si trae allora che, per ogniy €/,,, t€J;, le funzioni di J,: p(J,, ¢, 7), g(t, 7) 
sono assolutamente continue in {J,}. 


Fissato 7 in J,,, consideriamo la decomposizione di LEBESGUE di Pea (Jz,): 
(13) Poy Jes )= Pm Jes eae if Peon cS n) ag 


ove pou ( Jz,4) & la componente singolare di pg Hal Js). Posto, persyiGie-sc 
Ieee}: 


P* (It na? Pogo Je») : is ; <b 
la (43) puo scriversi: ; 
(14) P(Iy 7) = P* (Sy tn) =n) dy I Poslt’, 9,0) atl. 


Da questa relazione si trae che, per 7 CI, e FONE la p*(J,, ¢, 4) & funzione asso- 
lutamente continua in {J,}. Dato che, per ogni 7€J,, la Pom (Je> 7) & funzione 
di J, additiva e non negativain {f,}, pery€I,, J,c{J,} la funzione dit: p*(J,, t, n) 
é monotona non decrescente in J,. Inoltre la funzione p*(L,,1,¢,4) € misurabile 
in I. Infatti dalle (12) e dalla misurabilita delle funzioni O(/, ;,7,7) e V(Z,,;, ¢,%) 
si trae quella della funzione p(L, ;, t, 7) € fpaes quella della funzione Poy (t See (i 


E cosi provata la misurabilita di f dy’ I Boag , yn) at’ e quindi, per la (14), 
Ty, a 
anche quella di p*(J, ;, ¢, 7). Pertanto V’insieme E dei punti di J nei quali non 
esiste finito max lim 5 hab y,t,¢,4) € misurabile ed inoltre la funzione cosi 
definita in T: 3 
= max lim n’ TEP Tot i) etl aol, 


if n—>0o 
eam) in E 


é misurabile in J. 
La funzione g(t, y, 4), per ogni (y, 7) CL, XJ, € funzione di ¢ monotona non 
decrescente in J,; infatti, se t’/>¢, si ha: p*(Iy,2, t, 7) SO lb ek t’,y), donde: 


max lim n'~* p*(Iy,1 ,t, 9) = max lim n* p* (1,2, t', 7). 
n— CO 


n—> © 


Pertanto la derivata é non negativa. Inoltre m(f, y, 7) € non negativa dato 


che tale é Pome. Hye Loiche, per yell, rer, la p*(J,; t,) & funzione di f, 
assolutamente continua in { I la p(t, y, y) & funzione di y misurabile e somma- 
bile in J, e si ha: 7 

1] e ii. , t € I, 
*(J,,t,n) =J plu y.n) dy | 
b* yt) rka 7 Rely 


Riesce pertanto: ~ 
eeaaraee wary yav|a ee tel, 
me ee eaathag-adnin | tnotlal, (jh 


y 
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Facciamo ora vedere che Z olt, y,m) & quasi ovunque nulla in J. Infatti, 
se cosi non fosse, dato che Be é sempre non negativa, esisterebbero un punto 7 


contenuto in J,, un numero positivo Qo (y) ed un insieme E(y)CJ/;, di misura 
h-dimensionale positiva, tali che: 


é Het, 
= p(t, ¥,n) = eolt 
ap PH YM) 09 (M) oe 


Poniamo, per 7 €1,,: 
=Oo(n) se (6, y) €E(n) 
Ree 
=o se (t,y) Cl, — E(y). 
Per pel vel, t€I,, si ha allora: 
t t z 
fol.xaa's [sz 9l.y0) dt Sobyn). 


Dalla (15) si trae ha 
Poywn(Le, n) = n) 2 eI [Po WA’ ny n) Ts Q (t’, y; n) | dy at’ 


e cio é assurdo. 


Ripetendo per 74, le considerazioni ora svolte per po, ,,, Sl perviene a di- 
mostrare l’esistenza di una funzione y(t, y, 7) definita per (y, 7)c1, xJ, e per 
t€I,, la quale é non negativa, monotona non decrescente rispetto a ¢ in Tasha: 
derivata rispetto a ¢ quasi ovunque nulla, é misurabile e sommabile rispetto ad 
y in I, ed € tale che: 


(16) g(t) =) vit y, 0) +f Jou) (t's Vs) ar’ dy ; 
Perel, j,e{ je} st ha: 

Pam Se.) = Poym Jer) — Geum Je”) = 8 (F401) — 9,40). 
Dalle (14) e (16) segue allora, per 7 €J,,, cel. IRS 


J [wty¥.n) — w(a, y,n)] dy 
ty (3 
=e Ie 9,0) — vl yn) +S [Pog 9) — Poult, 90) dt} dy. 


Da questa relazione si trae l’esistenza di una funzione # definita in J, equi- 


valente a w (e quindi ad ) in J, tale che, per t€J,, (y, Miedo, satlesces 


"? 


u(t, y,n) — (a, ¥,n) = p(t, ¥,n) — pl, ¥ 7) +S [Pagal V1) — Goyay(t’s ¥, 0) | at’. 


Ne viene che “(¢, y, 7) ¢ VL rispetto a ¢ in I, per ogni (y, 7) €I yo Lye che quash 
ou 
ovunque in J, riesce: Toa xm, donde l’asserto, per 1</hh<r. Nel caso h=r 


la dimostrazione é analoga a quella ora esposta. 
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The Mechanics of Non-Linear Materials with Memory 
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A.E. GREEN, R.S. RIvLin & A. J. M. SPENCER 


1. Introduction 


In a previous paper (GREEN & RIVLIN 1957), subsequently referred to as 
Part I, the form of the constitutive equations governing the deformation of a 
class of materials possessing memory was discussed. It was assumed that the 
stress in an element of the material depends not only on the deformation gra- 
dients in the element at the instant of time considered, but also on those at 
previous instants of time. The limitations imposed on the constitutive equation 
by the fact that it must be form-invariant under a rotation of the physical system 
considered (consisting of the body and applied forces) were examined. 

This was done by first considering that the stress depends on the deformation 
gradients at a number of discrete times up to the instant of measurement. Then, 
the number of instants of time was considered to increase indefinitely, so that 
the expression for the stress became a functional of the deformation gradients. 
In this analysis it was found that the form-invariance of the constitutive equation 
under a rotation of the physical system leads naturally to a particular form of 
dependence of the stress on the deformation gradients at the instant of meas- 
urement. It was assumed that, apart from this, the expression for the stress as 
a functional of the deformation gradients at times up to and including the instant 
of measurement is continuous. 

In the present paper, we do not make this assumption, but allow that the 
stress may have arbitrary polynomial dependence on the deformation gradients 
at the instant of measurement, while its functional dependence on the deforma- 
tion gradients at times preceding the instant of measurement is continuous. 
Under these conditions, the limitations imposed by isotropy of the material in 
its undeformed state on the form of the constitutive equation is considered. 


2. Special Dependence of the Stress on the Displacement Gradients 
We consider a three-dimensional body to undergo deformation described in 
a fixed rectangular Cartesian coordinate system x by 
1(t) =%(X;,2) (0), 
x(t) = X; (t <0) 


v 


(2.1) 


where X; and x; are the coordinates in the system x of a generic particle of the 
body at zero time and time 7 respectively. We assume that x(t) are single- 
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valued functions of the arguments, possessing continuous spatial derivatives up 
to any required order except possibly at singular points, lines and surfaces. 
If the deformation is to be possible in a real material we must have 


|ax (2)/2X,| > 0. 


We assume the stress components a; ;(¢) at time ¢ in the system x to be poly- 
nomial functions of the deformation gradients 0x;(t,)/OX; («=0, 1, 2,..., N) at 
the NV --1 instants of time'ty)'t2)'% .., ty, °t)(=2) between +=Ocand r= 7%) It 
has been shown (GREEN & RIvLINn 1957), that o;; may then be expressed in the 
form 


1 Ox; Ox; 
Onn == Vg 1 6;; OX, ax, ie ’ (2.2) 
where the notation 
a Ox; _ 0%;(2) 
6;; =9;;(), ax, ax, ) 
Ol GAl CEA Gs 
Boalt) a Te Boy = pall (23) 


EA) epi 1 — 8), 


is used, 0;; denotes the Kronecker delta and f and /,, are polynomial functions 
Oine, (vq) andi) 2 (z,)(« 0,172; «2.5, N). 


We may re-write (2.2) in a somewhat more succinct form by observing that 


2¢ (t) Dy Bare Cymn suv 8mu Sno» (2.4) 
r s 
so that 
Vp) 2¢! ox, OX, 7S? (2.5) 
where 


LDR Cp Cte mae ne fo Ire: (2.6) 


Hence F,, is a symmetric polynomial function of g,,(t,) and Ve(t) («=0, 4, 
2 eeu 


3. Isotropic Materials 
In Part I the restrictions required for a material which is isotropic when 
™<O were obtained using (2.2) and the corresponding results for (2.5) can then 
be deduced*. Here, however, we proceed directly from (2.5). Let % be a fixed 
rectangular Cartesian coordinate system related to x by 
pe hel) (3.1) 


Ce ie Ee. 


where A;; are constants satisfying the orthogonality conditions 
A,, A;,=4,;,4,; = 6:3, |A;,| =r ie (3.2) 


x A slight gap in the argument in Part I was completed later (GREEN & RIVLIN 
1958). 


6* 
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If %;(7) denotes the coordinates of x;(7) in the system % and X; =%;(0), we have, 
from (3.4), 


8; ;(T) ea ae ax is Poles BaP) (3.3) 
& (t) = |8i;(t)| =8(7)- 
Denoting the stress in the system % by o;; we obtain 
t= Ni; “Gn Span (3.4) 
Hence, from (2.5), (3.1) and (3.4), 


5 1 0%; 0X; 
aa Ox ve Ape een nae (3.5) 


If the material is isotropic at t =O then 


= . 4 OX; OX; rs a 
Oi; az ax, ax, [Fe Fn glitevs Pe Wma) \a (3.6) 


so that, using (3.5) and (3.6), we have 


Ox; O07; 
Oxy ON. 


{F [Epo (tals VE (twice Ave vgtinn tenet tn) leah) Owen) 


Multiplying this equation by the non-zero expressions 


OX, @X, 
Ou pee 
we obtain pa, 
ie [opal a) Veg (T.) | == re aA 1 Oo ag (tz); Vg (T.) | (3.8) 
Using the notation 
g(t) =|lg:;()], g(t) =(1,, (2) |, (3.9) 


it follows that F, are the components of a symmetric matrix polynomial in the 
matrices g (t,) in which the coefficients are scalar polynomials in /g(t,) (« =0, 4, 
2, ..., N) and traces of products formed from the matrices g (t,) (« =0, 1, 2,...,N). 
Since [g(r)]! is a continuous single-valued function of 8pq(t) and 41/g? has no 
singularities, we can omit the factor 4/(2g3) and the arguments [g(z,)]? in (2.5) 
and write, to any desired degree of approximation, 


c= cre, (3.10) 
where 


o=|lo,,||, ¢ =||ax,/aXx,|). (3.14) 
In (3.10) €’ is the transpose of ¢ and F is a symmetric matrix polynomial in the 


matrices g(t,) in which the coefficients are scalar polynomial functions in the 
traces of products formed from the matrices Fi, O0 eee 
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The expression (3.10), valid for materials which are isotropic at time tXO, 
can be written in a convenient alternative form. Let* 


Np q(T) == Ax Ox, Pa ) Ox» 0%, 729 6 12) 
R(t) = ||Apq(x) |, 
and 
Ox; Ox. 

C=|C;,|, C= Ox, aX, (3.13) 

Then 
C— 66,19 (0) =@ hit), hit) =e) 7 g(r)-c2: (3.14) 
The trace of products of matrices formed from g(t,) («=0,1, 2,..., N) 
can be expressed as the trace of products of matrices formed from C and h(r,) 
(x=0,1,2,...,.N). Also a symmetric matrix polynomial formed from the 


matrices g (t,,) can be expressed in the form e’ Le where L is a symmetric matrix 
polynomial formed from the matrices C and h(r,). It follows that we may write 
(3.10) in the alternative form 

o=9/C,h(z,)], (3.15) 


where o is a symmetric matrix polynomial in the matrices C,h(t,) («=1, 
2,...,.N), with coefficients which are polynomial functions of the traces of 
products formed from the matrices C, h(t,) (x =1,2,..., N). 


Conversely, we can show-that any constitutive equation of the form (3.15) 
can be expressed with any desired accuracy in the form (3.10), where F is a 
symmetric matrix polynomial in the matrices g(t,) with coefficients which are 
scalar polynomial functions in the traces of products formed from these matrices. 


4. Isotropic Tensor Functional 
When the material is isotropic initially we may make the transition from 
tensor functions to tensor functionals using either (3.10) or (3.15) as a starting 
point. In order to keep in line with the work of Part I we use (3.10) and suppose 
that 
Oj7=so- ae Gs) (4.1) 


where we assume that F, , is a symmetric matrix polynomial in g,, with coefficients 
which are functionals of g,,(t) in the range 0StSi. We assume that these 
functionals are continuous functionals of g,,(t) over the compact aggregate of 
functions which are continuous in the range 0<tSX?¢. We denote the Fourier 
half-range cosine coefficients of g,,(t) by Gy where 


t 
Gp) = Ff epqle) cos“ F2ax (a0), 
0 


* From (2.1) %,,(t) may be expressed as a function of %,,(¢), t, ¢ by eliminating 
X;, and then @*,,(t)/@%» can be evaluated. Alternatively, we can use the first expres- 
sion in (3.12) to find hy,(z). 
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Then each continuous functional in the expression for F., may be approximated 
by a polynomial P (G)?), Gi), ..., GY) which tends uniformly to the functional 
as g—> oo. It follows that F, can be expressed with any desired approximation 
by a polynomial in g,, and Gi), Gi), ..., Gil) with a sufficiently large choice of 
N and we write 
HE EBs Oma (N) 
Yi CIE ax, tes 8pqr Cpa Gpar---» Goal: (4.3) 
where F,, is a polynomial. 
We next consider a change of rectangular Cartesian axes of the form (3.1) 
and (3.2). If 6,; denotes the stress components at time ¢ in the system *, then 
G;; is given by 3. 4) and pitas Also, if the material is isotropic for tO, then 


Ea Epps Ghpnepy et labels (4.4) 


where 


oy (4.5) 


and g,,(t) is given by (3.3). From (3.3), (4.2) and (4.5), we have 


&rqg =Ap Move Smn> (4 6) 
Cte imbue ental ) 
Also, from (3.1), (3.3), (3.4), (4.3) and (4.4), we obtain 
On; OX; val 
t = («) a 
SB Ge Fs lBoa GP) — Arm Asn Ennl8pas GST} = 0 
and hence 
LONG ge y f Cy Ay mas wlan [ge , ey S 4.7 
Using the notation shee 7 ae oe 
G,=IGfl (48) 
it follows from (4. 7) and (4.6) that F is a symmetric matrix polynomial in the 
symmetric matrices g, G, («=0,1,...,.N) and equation (4.3) may be written 
in the matrix form 
o=cF ec’, (4.9) 
where F is a symmetric matrix polynomial in the matrices g, G, (a =0,.4, 2005 Ve 
Since g(t) =e'h(t)e, it follows that G,=c'H,,¢, where 
#1, = lV 
H®) — 2 
=F to cos aee dt (a > 0), (4.10) 


H®) = ne dt 
0 
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Hence, by an argument similar to that used in § 3, we may reduce (4.9) to the 
alternative form 
o =e(C, Hy, Hy 4. 4g (4.11) 


where o is a symmetric matrix polynomial in the arguments stated and C is 
given by (3.13). 


5. Further Development of the Equations for Isotropic Materials 
Since F in (4.9) is a symmetric matrix polynomial in the matrices g, G, 
(x0, 1,..., N) it is apparent that it may be expressed as the sum of a number 
of terms of the form 


te t 
/ (2) Shere Sea NL een A ALE TW? reece SO 
fife 8% on 09% co 8 ri 
x (tril, trIl,...trI,) 14 I’) diydt,...dtp, 


where II, , II,,..., 11, and II are matrix products formed from g, g(t), g(t2), 

--,9(tp) and are such that IL, II, ... MII is linear in each of the matrices g (1), 
9 (T:), ..., g(Tp), and II’ denotes the transpose of II. Let II, be a matrix product 
formed from g(t), 9(T2), ---, J (T,,) and possibly g; let II, be a matrix product 
formed from g (Tg, +1), J (T,+2)) --+» Y(T,,) and possibly g and so on, so that IL, 
is a matrix product formed from g (T,,_, +1), J (Tay-y+2)> +++» J (Tp,,) and possibly g, 
and let IT be a matrix product formed from g (T,,,+1), 9 (Ts, +2), +++» 9 (Tp) and 
possibly g. Then, the term (5.1) may be re-written as 


Coe t 

Ap IT, 

f() f f---f cos “1% *1 cos BATE... cos Str dr...d%, 
0 0 0 


t ¢ t 
. 4 aT, ag, TB, 
f feof cos 278 Att .. .cos — — trIE 25,41 -.. 4%, 
0 0 0 


. ° . 5 . . ; (5.2) 
t ¢ t 
t 
iy Ufrcedp cos Bact) 7 Touma i ++» COS Mn Ob te TL As ee tn hap 
0 0 0 
eae : Bu +1 % TBy+1 ap Tp / 
f f---f cos EE cos “PE (WW) dtp, 41 --- 4 tp. 
Os.0 0 
We now consider the factor 
bat t 
Op, UT, 
f f---f cos “1271 cos “BET ---cos =! tr dt, ...d%,. (5.3) 
0 0 0 


It has been shown (SPENCER & RIVLIN 1959) that the trace of a matrix 
product formed from the R symmetric 3 x3 matrices @p (P =1, 2, ..., R) may be 
expressed as a polynomial in the traces of matrix products formed from those 
listed in Table 1 by replacing K,, K,,..., Kz by all permutations of 1, 2,...,R 
seven at a time. 
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be expressible as a polynomial in the invariants listed in Table 2 and the invariants 
formed from these by replacing 1,, T2,.-., 7 by all possible permutations of 
Ty, T2,--+, Tg, Seven ata time. 
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Since II, is linear in each of the matrices g(t), 9 (T2), +++; g (tz,), tr TI, must 


i 


Table 1 
trax,, tr ak,, tr ak; 
tr ax, 4x, tr ax, ak,, tr ak, ak,; 
trax, ax, Ax, tr ax, ax, Ak, tr ax, ak, ak, ; 
trax, 4x,4x,4x,, tt ax, 4x, 4x, ax, tr ax, ax, 4k, 4k,, 
tr ax, ax, 4x, 4K, 4K, 
tr ax, ax, x, UK, UK, tr ax, 4x, 4x, IK, aK,, 
trax, 4x, 4K, 4K,4K,9%K,, tt ax, 4x, 4K, 4K, 4K, ak,, 
tr ax, ax, 4K, 4K,4K,%x,, tt ax, 4x, ax, ax, ak, 
tr ax, 4x, 4x, 9K, 4K,4K,, tT 4x, 4x, 4x, 4K, AK, ak, 


tr ax, Ax, Ix, UK, 4K, 4K, IK,- 


Denoting the invariants in Table 2 which involve none of the matrices 


g(%), J (Ts), ---, 9(T7) by tr@®® («=1, 2,3), those which involve g(t,) only 
or g and g(z,) only by tr & («=1, 2, 3), those which involve g(t,) and g(t.) 


Table 2 
CEG tL Oem ath ges 
trgg(t%), trg’?g(r); 
trg g(t) g(t), trg?g(%) g(t); 
tr g 9(%) 9(T2) g(t), te g? G(T) Y(T.2) G(T), 
tr g(%) 9(T2) 9 G(T) g?; 
tr g 9(%) 9(T2) 9 (Ts) G(T), tr g? 9 (7%) 9 (T2) G(T) G(T), 
tr g(%) 9 G(T) 9(Ts) G(T) 9%, = tr G(T) G(T2) 9 G(Ts) G(T) 9, 
tr g(t) 9 (72) 9(t3) 9 G(T) G?; 
tr g 9(%) G(T) G(T) G (Ts) (Ts), tr g? g(t) 9(T2) G(T) Y(Ts) G (Ts); 
) 


9 \%3) g 
tr g 9 (T,) 9(T2) G(T3) F(T) G(T) G(T.) ; 


trg(%), trg(%) g(t), tr g(t) g(t2) g(t), 

tr g(%) G(T) G(%) G(T), te G(T) G(T.) G(Ts) G(Ts) G (TS), 
tr g(%) G(T) (Ts) G(T) G(T) G(T), 

tr g(t) 9(T2) 9 (Ts) G(T) g (T5) G(T) 9 (7)- 


or J, g(t,) and g(t.) only by tr &?) (x =1, 2, 3), and so on, we see that the term 
(5.3) may be expressed as a polynomial in expressions of the forms 


t aa 2 
tr &.°), iu p(t) traMdr, ff fo(%, t.) @O?) dt, dt, 
00 


(5.4) 


t t 
era @ (Ts Vay ce-4 Tq) UROL @ bed Gamera tee 


Where the functions ” (ti), P(%, T.), --+, P(T%, Te, ---, T) are analytic functions 
of their arguments. Each of the other factors in (5.2) except the last may be 


similarly expressed. 
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We now consider the last factor in (5.2), 7.e 
t ¢t t 
J ff cos Fett Set cos SPATE (TT HID’) dtg,41-..dtp. (5.5) 
0 oO 


Since II+II’ is a symmetric isotropic matrix polynomial in the matrices 
9 (Tpyu+1)> J (Ten+2)) ++» 9 (Tp) and g, linear in each of them except g, it follows 
(SPENCER & RIVLIN 1958) that it may be expressed as an isotropic matrix poly- 


nomial of the form ; 
+ I’ => 9,(%+%), (5.6) 


where x, denote the matrix products formed from those listed in Table 3 by 
replacing 7, T2,..-,% by all possible permutations of 1,11, Tg, +9) +--+) Tp SiX 
at a time, while gw, are polynomials in the invariants obtained from those listed 
in Table 2, by replacing 11, t2,..-,t, by all possible permutations of T,,4,, 
Teu+2»+++, Tp seven at a time. x’ denotes the transpose of x. Since II + II’ is 
linear in each of the matrices g (t,,+1), 9 (Tay+2)> +++» J (Te), P, and X,+ x, cannot 
involve any of these matrices in common. We may therefore express the factor 
(5.5) as an isotropic matrix polynomial in which the matrix terms are of the forms 


10, fooler P+ xO) a 

te FE 

el Gps ta) a (5.7) 
00 


t ¢t 


IE rece cavclas (x + M0) dt dt,... dt, 
00 0 


where x) are the matrix products listed in Table 3 which do not involve any 
of the matrices g(t), g(T2), ---,9(T,.), x are those which involve only g(t) 


Table 3 
q; 
Gins 9; 
g(u), 9? 9(t%); 
9(%) 9(%), g(ugg(t%), Gg? g9(%) g(t), g(t) 9? G(T), 
g ae 9? g(t), 9 9(%) G(T.) 9; 
9 9(%) G(T) J (Ts), 9(%) 9 G(T2) J (Ts), 
9° G(%) G(T) G(T), 9 (%) 9? G(T) G(T); 
9 9(7,) G(T2) J (Ts) g?, 9(%) 9 9(T2) G(T) g?, 
9(%) G(T) 9 G(t3) 9, G(T) FI(T2) 9? G(T); 
9 G(T) G(T) G(T) F(T), 9 (7) 99(T) g (t3) 9 (Ta) 
9(%) G(T) GGT) G(T), — 9? G(T) F(T) F(Ts) Y (Ta), 
9 (7%) 9? G(T) G(T) G(T), G(T) F(T2) 9? G(T) Y (Ta); 


) 
9 9 (7) G(T) G(T) G(T) G(Ts), F(T) Y Y(T2) F(T) G(T) YTS)» 
9(%) 9 (72) G G(s) G(T) G(T); 
9(%), 9(%) G(T), G(%) G(T) G(T), G(T) G(T) Y (Ts) G(T), 
9 (%) 9 (Tz) G(T) G(T 4) 9 (Ts), G(T) G(T) Y(T3) g (Ta) 9 (Ts) G(T) - 


or g and g(t,), x? are those which involve only g(t) and g(t) or g, g(t) and 
g(t) and so on, and the functions y are analytic functions of their arguments. 
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It follows that (5.2) and hence F may be expressed as an isotropic matrix 
polynomial in which the matrix terms take the forms (5.7) and the coefficients 
are polynomials in the invariants (5.4) and functions of ¢. In general, the matrix 
polynomial will contain more than one term of each of the forms (5.7), except 
the first. Alternatively, we can bring the coefficients under the integration signs 
in the terms of the forms (5.7) to derive the result that F may be expressed 
in the form 


t 
Bis OTs! OG. Yr ek 
v vy 0 


tt 
+ DS LA, (er, te) P+ XP’) dy dty+ ++ + (5.8) 


0 

t t 
J: F f 3, (t%, To, eon te) (x + x’) Gt at, 5 4U,; 
0 0 


where the #’s are functions of ¢ and of their indicated arguments and polynomials 
in the invariants (5.4). 
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Mathematische Begriindung der Scheibentheorie 
(xweidimensionale Elastizitatstheorie) 


DIETRICH MORGENSTERN 


Vorgelegt von ELI STERNBERG 


I. Problemstellung 


In der Begriindung der Elastizitatstheorie gewinnt man aus dem Hookeschen 
Gesetz unter gewissen, bekannten Vernachlassigungen die Fundamentalgleichun- 
gen der linearen dreidimensionalen Elastizitatstheorie: Fiir das Spannungsfeld o; , 
und das Verschiebungsfeld , gilt mit dem daraus abgeleiteten Verzerrungsfeld! 


(1) Vir =% (Win + Ups) » 

rg oh suilereen in 210i 
(2 a) Vik 2G OrR 0; k Em 
bzw. die Umkehrung 

G 

(2b) Cin = 26 Yin + 6:4 (21711) 
und 
(3) Cri 
Elimination ergibt dann 
(4) G {eth He oe Hiren} = 0. 


Dabei gilt zwischen den drei positiven Konstanten E (Elastizitatsmodul) G 
(Gleit- oder Schubmodul) und m (Querkontraktionszahl 22) die Relation 
2G(m+1)=Em. Zur Bestimmung des Spannungs- und Verschiebungszustandes 
gehoren zu den Differentialgleichungen (1) bis (4), die Innern eines Bereiches 
gelten, noch Randbedingungen: 


a) ”;0;,—=/, an den Teilen der Oberflache (mit der Normalen n;), wo die 
(5) Oberflachenkraftdichte /, gegeben ist. 


b) vorgeschriebene Werte fiir ~;, wo die Verschiebung gegeben ist. 


1 Indizes bei uv, o, y durchlaufen 1, 2, 3; bei den spateren Buchstaben 1, T, ¢ 
durchlaufen die Indizes 1, 2. Summationskonvention und Ableitungsbezeichnung 


w wird verwendet. 


aa: Ox; 
L 
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Scheibenformige Kérper, das sind Zylinder (Achsenrichtung sei %,) mit geringer 
Hohe h, die nur an dem Mantel des Zylinders (dem Rand der Scheibe) durch 
Krafte und Verschiebungen in der %,, ¥,-Richtung beansprucht werden (/; =0) — 
kurz Scheibe2 genannt — sind fiir praktische Anwendungen wichtig; hier gelten 
also auBer den Differentialgleichungen (1) bis (4) in dem Bereich 


{x,, %EG, |x| Soh 
die Randbedingungen 
a) o,=0 fir %=+ 3h; {%, %fE@ 
(5%) 1; 6;, = Ty fit 4x, tae, [%e\ee eg mut fy = 0- 


b). wget fir Ln, X76 Sine Sale 


Dabei bedeuten %,, Jt, eine Zerlegung des Randes des zweidimensionalen Be- 
reiches ©. Der Einfachheit halber sollen dabei die in den Randbedingungen 
auftretenden Funktionen /, und u;* als unabhaingig von x, angenommen 
werden. 


Fiir technische Zwecke gewinnt man nun zweidimensionale Gleichungen zur 
Lésung dieser Randwertaufgabe; dazu setzt man o;, = 0, laBt die Gleichungen (3) 
mit 7=3 weg und nimmt alle Funktionen als nur von x, und x, abhangig an. 
Es entsteht dann das Gleichungssystem 


(6) v; = Verschiebung, Verzerrung ¢;, = $ (Uj, + Ugh), 

4 
(7a) Cie = ZG Tin — Oih — 
bzw. die Umkehrung 

G 
(7b) Tp = 2603, + obra (Db) 
und 
(8) Tip =O; 
was durch Elimination auf 
1 

(9) G {in ae ms vaya} =0 


fiihrt. Diese Differentialgleichungen (9), die Gleichungen der ,,zweidimensionalen 
Elastizitatstheorie“ gelten im Innern des zweidimensionalen Bereiches @; dazu 
kommen noch Randbedingungen: 


a) n;T;,—=f, an dem Teil R, der Randkurve von ® wo die Randkrafte t, 
as gegeben sind (k =1, 2), 
10 
b) v; =u an dem komplementaren Teil Jt, des Randes, wo die Verschie- 
bungen gegeben sind. 


Iho Gegensatz zur technischen Platte, die auch — oder nur — seitliche Krafte 
(in der ¥,-Richtung) aufnehmen soll. 
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Es ist das Ziel dieser Untersuchung, diese durch Plausibilitatsbetrachtung (die 
insbesondere die Wahl der vernachlassigten Gleichung betrifft) gefundenen Glei- 
chungen als wesentlich fiir diinne Scheiben nachzuweisen. Es wird dazu gezeigt 
werden, daB bei h->0 die tiber die Scheibendicke gemittelten Werte der Lésungen 
des Randwertproblems (4) (5) gegen die Lésung von (9) (10) streben. 


II. Formulierung als Extremwertaufgabe 


Fir die dreidimensionale Elastizitatstheorie ist als Prinzip vom Minimum 
der Formanderungsarbeit bekannt, daB die Verschiebung durch die Extremums- 
eigenschaft 


(11) As( =6f[ [| Dy Sh, rw)'| dada dy — [fu f,d0 = Minimum 
i,k 
Ri 
(do = Oberflachenelement) 


charakterisiert werden kann; dabei werden alle hinreichend oft differenzierbaren 
u-Felder, die den geometrischen Bindungen (5*) geniigen, zugelassen. Der Mini- 
malwert bei der Scheibendicke h sei a, (h). 


Analog gehért zu den Gleichungen (9) (10) das Extremwertproblem 


(12) A,{v} =6/{[2ce+ 4 (Uo) da day — [fas = Minimum 


m—1 
Ri 


mit dem Minimalwert a,. 


Diese Charakterisierung der Losungen durch Extremwerteigenschaften soll 
nun nach einem friiheren Vorbild®, welches beispielsweise die erste Randwert- 
aufgabe bei hAAu—Au=0 betraf, fiir den Beweis der behaupteten Limes- 
eigenschaft benutzt werden ?. 


III. Aufstellung zweier Ungleichungen 


In A, werde 
4h 
(13) v= > | ud Xs 


th 


als Vergleichsfunktion eingesetzt; dann folgt 


3 MORGENSTERN, D.: Singulare Stérungstheorie partieller Differentialgleichungen. 
J. Rat. Mechanics and Analysis 5, 203—216 (1956). 

4 Die Existenz von Lésungen hat auf Grund der zugehérigen Extremwerteigen- 
schaften K. FRIEDRICHS nachgewiesen: On the boundary value problems of the theory 
of elasticity and Korwn’s inequality. Annals of Mathematics 48, 441—471 (1947). 
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und aus der Schwarzschen Ungleichung erhalt man 


1h th 


(14) Ga jf rat vinid. (5 Gay)? << (> 1711)? 4 Xs, 


ere 


und es ergibt sich fiir 


(15) A,{ys + ef ff [27 si wor (3m) Ax, AXA Xs — ffm j,do} . 


Nun gilt offensichtlich 
» (16) 


1 2< ye 1 2 
Se Sige ee 


und es bleibt also die folgende Ungleichung 


A,{o} S > effi lz” 


a \2 
5 (2% + 33) x 
x dx, dx, dx, — [ jrusao} = 5 As {u}. 


Umgekehrt soll nun aus jedem Funktionenpaar v,, v, ein v;-Funktionensystem 
erzeugt und in den Ausdruck fiir A, eingesetzt werden: 


(18) u 
dann gilt 


=r (¢ =A5 2) Uy a Oe 


Vin =Sin(s 2 =1,2), Yig=3%Mji(0=1,2), Yys=w, 


und man erkennt, daB man nach der im Sinne der Ungleichung (16) optimalen 
Substitution 


(19) 


2 
1 
 m—4 2b 


da die Integration in der x,-Richtung leicht durchfiihrbar ist, die Ungleichung 


patacef]|$ erat (helt $ tentay 


ee Dee 


x Ox, 4%, — ee Ag{v} +62 Lf (grad w)?dx,d x5 
R, 
erhalt. 


IV. AbschluB des Beweises 


Tecunare erkennt man aus der Ungleichung (17), da8 fiir die aus der 1 Osis 
) fiir die Scheibe der Dicke h durch (13) gewonnenen Funktionen v™ gilt 


U 1 L 
(21) A, {ol} <— Ay {ul} = + a; (h). 
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Andrerseits gilt wegen (20) fiir das aus der Lésung des zweidimensionalen Pro- 
blems v gewonnenen u 


4a; (h) = poe 7 As{u} S 4a{o} +65 [ (grad w)?dx, dx, 
(22) 


= a, + CF | [ (grad wax dy, 


wobel w = — 


2 
7 >» ¢,, einzusetzen ist. 
1=1 


Wenn also dies letzte Integral konvergiert, erkennt man 


1 
(23) lim —- as (i) S a, 


_und wegen A,{v™} => a, und (21) auch 


(24) jim A, {ul} = ay. 


Falls fiir die Losung des zweidimensionalen Problems das Integral /f (grad w)?dx, dx, 
divergiert, kann man statt dessen eine Naherungsfunktion verwenden, fiir die das 
Integral endlich bleibt, aber a, bis auf eine beliebig kleine Differenz erreicht 
wird. Es folgt dann dasselbe Ergebnis. 

In dem Hilbertschen Raum der Funktionenpaare v,,v, mit der Dirichlet- 
Norm 


(25) || 2|| = (@, v) =G/f)xe+ 4 Oe |dudx 


1aBt sich A, so ausdriicken 


(26) A, {v} = (u,v) — 2(v, g) = ||» — g|? — lel, 


wobei g zu der Loésung des zweidimensionalen Problems gehért. Fiir die Ab- 
stande von diesem Element g gilt dann, wenn Mt die lineare Menge der erlaubten 
Funktionenpaare bezeichnet, die die Randbedingungen (10) erfiillen, 


(27 inf ls — gl? =a, + [lell?=lle — 2 
und 
(28) |e — gl? 4e.+|[e]F oem), 


woraus nach dem Lemma vom eindeutig bestimmten Lot auf konvexe Mengen 
die Konvergenz vv im Sinne unserer Norm (25) folgt. 

Damit ist also bewiesen, daB die nach (13) gemittelten Werte der Lésungen 
des dreidimensionalen Problems bei h-—>0 im Mittel gegen die Lésung der Glei- 
chungen (9) (10) konvergieren®. Aus dieser Mittelkonvergenz im Sinne der Norm 
(25) folgt, wie leicht gezeigt werden kann, auch die Mittelkonvergenz im iiblichen 
Quadratmittel, wenn man 2, v2 nur bis auf Transformationen v,>v,+4%,+ 4, 
Ug—>V,—a4x,+C, die infinitesimalen Bewegungen entsprechen, als bestimmt an- 
sieht. 


5 Diese zweidimensionalen Gleichungen sind nicht zu verwechseln mit den Gln. (4) 
fiir x,-unabhangige Spannungszustande (unendlich dicke Scheiben), wo alle GroBen 
von *, unabhangig sind. 


vy. SchluBbemerkung 


5 4 AuBer der sich hier anschlieBenden Frage nach der Re Contes Koirae 
_ die im Anschlu8 an dieses Ergebnis zu losen méglich sein diirfte, gibt es eine 
-_ Reihe anderer begrifflicher Liicken beim Aufbau der Mechanik, insbesondere der 
in der ,,technischen Mechanik‘‘ verwendeten N aherungstheorien, die durch ahn- 
_ liche mathematische Grenzwertsatze geschlossen werden sollten; so zum Beispiel 
fehlt eine exakte Begriindung fiir die technische Naherungstheorie der Balkens 
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